Planar stability criteria for viscous shock 
waves of systems with real viscosity 



Kevin Zumbrun* 
February 1, 2008 



Abstract 

We present a streamlined account of recent developments in the 
stability theory for planar viscous shock waves, with an emphasis on 
applications to physical models with "real," or partial viscosity. The 
main result is the establishment of necessary, or "weak" , and sufficient, 
or "strong" , conditions for nonlinear stability analogous to those estab- 
lished by Majda [Ma. 1-3] in the inviscid case but (generically) sepa- 
rated by a codimension-one set in parameter space rather than an open 
set as in the inviscid case. The importance of codimension one is that 
transition between nonlinear stability and instability is thereby deter- 
mined, lying on the boundary set between the open regions of strong 
stability and strong instability (the latter defined as failure of weak 
stability). Strong stability holds always for small-amplitude shocks of 
classical "Lax" type [PZ.1-2, FreS]; for large-amplitude shocks, how- 
ever, strong instability may occur [ZS, Z.3]. 
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1 Introduction: structure of physical equations 



Many equations of physics take the form of hyperbolic conservation laws 

(1.1) u t + J2 Fj (u) Xj =o, 

j 

with associated viscous conservation laws 

(1.2) D i ^i/i^ir)/:,,;,. 

incorporating neglected transport effects of viscosity, heat conduction, etc. 
(more generally, hyperbolic and viscous balance laws 1 including also zero- 
order derivative terms C(U), as especially in relaxation and combustion 
equations). Here, U, F J G W\ B^ k G M nxn , x G R d , and t G K. 

Examples 1.1. Euler and Navier-Stokes equations, respectively, of gas- 
or magnetohydro dynamics (MHD); see Q1.44J1 below. 

A fundamental feature of is the appearance of shock waves 

(1.3) U(x,t) = U(x - st) 



U- x\ < st, 
U+ x\ > st, 



discontinuous weak, or distributional solutions of (|1.1|) determined by the 
Rankine-Hugoniot conditions 

(RH) s[U] = [F(U)}; 

see, e.g., [La. 1-2, Sm]. Here and elsewhere, [/i(Z/)] := h(U+) — h{U ). Solu- 
tion ()1.3() may be uniquely identified by the "shock triple" (£/_, U + , s). 

Such waves in fact occur in applications, and do well-approximate ex- 
perimentally observed behavior. However, in general they occur in only one 
direction, despite the apparent symmetry (x, t, s) — > (— x, —t, s) in equations 
(|1.1|) . That is, only one of the shock triples (£/_, U + , s) (U+, U-, s) is typi- 
cally observed, though they are indistinguishable from the point of view of 
(RH). The question of when and why a particular shock triple is physically 
realizable, known as the shock admissibility problem, is one of the oldest and 

1 Outside the scope of these lectures, but accessible to the same techniques; see, e.g., 
[God, Z.3, MaZ.l, MaZ.5, Ly, LyZ.1-2, JLy.1-2]. 
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most central problems in the theory of shock waves. For an interesting dis- 
cussion of this issue from a general and surprisingly modern point of view, 
see the 1944 roundtable discussion of [vN]. 
Two basic approaches to admissibility are: 

1. Hyperbolic stability in the Hadamard sense, i.e., short-time bounded 
stability, or well-posedness of (|1.3|) as a solution of 1)1.1(1 . also known as dy- 
namical stability [BE]: that is, internal consistency of the hyperbolic model. 
Here, there exists a well-developed theory; see, e.g., [Ma. 1-3, Me. 1-4, FMe] 
and references therein. 

2. Consistency with viscous or other regularization, in this case the 
viscous conservation law (|1.2|) . (a) A simple version is the "viscous profile 
condition" , requiring existence of an associated family of nearby traveling- 
wave solutions 



of the viscous conservation law p.2jl : see, e.g., [Ra, Ge, CF, Be, Gi, MP, Pe, 
MeP], and references therein. This is the planar version of the prepared- 
data "vanishing-" or "small-viscosity" problem (SV) treated for general, 
curved shocks in the parallel article of Mark Williams in this volume [W]. 
The viscous profile condition, augmented with the requirement that U be 
a transverse connection with respect to the associated traveling-wave ODE, 
is sometimes known as structural stability [BE, ZS, Z.3-4]. (b) A more 
stringent version is the "stable viscous profile condition" , requiring stability 
under perturbation of individual profiles H1.4j) with viscosity coefficient v 
held fixed. We denote by (IT) the associated problem of determining long- 
time viscous stability. 

Definition 1.2. Long-time viscous stability is defined as the property that, 
for some appropriately chosen norms \ ■ \x and \ ■ \y , for initial data Uq 
suffiently close to profile U in \ ■ \x, the viscous problem (|1.2I) . v = 1, has a 
(unique) global solution U(-,t) that converges to U as t — > oo in \ ■ |y. We 
refer to the latter property as asymptotic \ ■ \x — ► | • |y viscous stability. 

Remark 1.3. One may also consider the question whether solutions of 
(|1.2|) converge on a bounded time interval for fixed initial data as v — > to 
a solution of 1 . 1|) . that is, the unprepared-data (SV) problem. This was 
considered for small-amplitude shock waves in one dimension by Yu [Yu], 
and, more recently, for general small-variation solutions in one dimension in 
the fundamental work of Bianchini and Bressan [BB] ; in multiple dimensions 
the problem remains completely open. This is a more stringent requirement 



(1.4) 




lim U(z) = U± 
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than either of 2(a) or 2(b); indeed, it appears to be overly restrictive as an 
admissibility condition. In particular, as discussed in [Fre.1-2, PreL, L.4, 
Z.3-4], there arise in (MHD) certain nonclassical "overcompressive" shocks 
that are both stable for fixed v and play an important role in solution 
structure, yet which do not persist as v — > 0. The (LT) and (SV) problems 
are related by the scaling 

(1.5) (x,t,u)^(x/T,t/T,u/T), 

with T — > oo, < t < T, the difference lying in the prescription of initial 
data. 

Conditions 1 and 2(a) may be formally derived by matched asymptotic 
expansion using the rescaling Q1.5|) and taking the zero- viscosity limit, as 
described, e.g., in [Z.3], Section 1.3. They have the advantage of simplicity, 
and for this reason have received the bulk of the attention in the classical 
mathematical physics literature; see, e.g., the excellent surveys [BE] and 
[MeP]. However, rigor (and also rectitude; see Section 1.4, [Z.3]) of the theory 
demands the study of the more complicated, but physically correct condition 
2(b), motivating the study of the long-time viscous stability problem (LT). 
It is this problem that we shall consider here. 

In contrast to the hyperbolic stability theory, progress in the multidimen- 
sional viscous stability theory has come only quite recently. The purpose of 
this article is to present an account of these recent developments, with an 
emphasis on (i) connections with, and refinement of the hyperbolic stabil- 
ity theory, and (ii) applications to situations of physical interest, i.e., real 
viscosity, large amplitude, and real (e.g., van der Waals-type) gas equation 
of state. Our modest goal is to present sharp and (at least numerically) 
computable planar stability criteria analogous to the Lopatinski condition 
obtained by Majda [K, Ma. 1-3] in the hyperbolic case. 

This is only the first step toward a complete theory; in particular, eval- 
uation of the stability criteria/classification of stability remain important 
open problems. Preliminary results in this direction include stability of gen- 
eral small-amplitude shock profiles [Go. 1-2, HuZ, PZ.1-2, FreS]; geometric 
conditions for stability, yielding instability of certain large-amplitude shock 
profiles [GZ, BSZ, FreZ, God, Z.2-4, Ly, LyZ.1-2]; and the development 
of efficient algorithms for numerical testing of stability [Br. 1-2, BrZ, BDG, 
KL]. On the other hand, the techniques we use here are completely general, 
applying also to relaxation, combustion, etc.; see, e.g., [Z.3] and references 
therein. 
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We begin in this section with some background discussion of a mainly 
historical nature concerning the common structural properties relevant to 
our investigations of various equations arising in mathematical physics, at 
the same time introducing some basic energy estimates of which we shall 
later make important use. 

1.1. Symmetry and normal forms. We may write (|l.lj) and (|1.2|) in 

quasilinear form as 

(1.6) U t + Y, Aiu xi=° 



and 

(1-7) U t ■ =uJ2(B jk U Xk ) Xj , 

j j,k 

where A 7 := dF^{U) and B^ k := B 3k (U). We assume the further structure 
u 7 \ ( A\, A\\ „„- t /0 



(w > "-•-»>• &)' "-wis 



u 1 £ W 1 r , u 11 G W typical in physical applications, identifying a distin- 
guished, "inviscid" variable u , with 

(1.9) ^J2^u>0\^, 

9 > 0, for all £ E M d . Here and below, aM denotes spectrum of a matrix or 
linear operator M. In the case of an unbounded operator, we use the simplest 
definition of spectrum as the complement of the resolvent set p(M), defined 
as the set of A € C for which A — M possesses a bounded inverse with respect 
to a specified norm and function space; see, e.g., [Kat], or Appendix A. 

1.1.1. Inviscid equations. Local stability of constant solutions (well- 
posedness) requires hyperbolicity of (|1.6j) . defined as the property that 
:= Ylj have real, semisimple eigenvalues for all £ € R d . As pointed 
out by Godunov and Friedrichs [G,Fr], this may be guaranteed by symmetriz- 
ability, defined as existence of a "symmetrizer" ^4° such that A is symmetric 
positive definite and A 3 ' := A A 3 are symmetric. Left-multiplication by ^4° 
converts (|1.6j) to (quasilinear) symmetric hyperbolic form 

(1.10) A°U t + Y,&U Xj =Q. 



6 



A nonlinear version of this procedure is an invertible coordinate change 
U — > W such that considered as an equation in W takes the form 

(1.11) u(w) t + J2 Fj ( u (w)) Xj =0, 



where ^4° := dU /dW is symmetric positive definite and A 3 := dF 3 {dU / dW) 
are symmetric. This is more restrictive, but has the advantage of preserving 
divergence form; see Remarks 11.121 1-2 below. 

Symmetric form yields hyperbolic properties directly through elementary 
energy estimates/integration by parts, using the Friedrichs symmetrizer re- 
lation 

(1.12) &(U,SU Xj ) = -\(U,Si j V) 

for self-adjoint operators S £ C nxn , U € C n (exercise). Here and below, (•, •) 
denotes the standard, (complex) L 2 inner product with respect to variable 
x. We shall require also the following elementary bounds. 

Lemma 1.4 (Strong Sobolev embedding principle). For s > d/2, 

(1-13) l/lL°°(R d ) - I/Il 1 ^) - C\f\H*(M. d )> 

where f denotes Fourier transform of f . 

Proof. The first inequality follows by Hausdorff- Young's inequality, the sec- 
ond by 

\f(0\^ = 1/(0(1 + m(i + l£ITV 

<i(i+iciTW(0(i+ier)ii*<ci/itf- 

□ 

Lemma 1.5 (Weak Moser inequality). For s = ^ \ctj\ and k > d/2, 

r 

l(a ai wi)---(a a ^ r )| L 2 (Rd) <J2\vi\ H s iRd) (J[\vj\ L 

(1-15) 

<C^|i; i | fr . (R dj(jJ|^| £r * 
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Proof. By repeated application of the Hausdorff- Young inequality \f*g\]j> < 
|/|i2 \g\u>, where * denotes convolution, we obtain the first inequality, 



k^o • • • {d^v r )\ L , = Kr^i) * • • • * (r^r)i 



L 2 



(1.16) 



< \^ s vi | * — * \v r 



L 2 



+ ■•• + 



* • • • * \£, s v r 



L 2 



i=l 



The second follows by (f 1 . 13|> . 



□ 



Proposition 1.6 ([Fr, G]). Symmetric form (|1.1U() implies local well- 
posedness in H s , and bounded local stability \U(£)\h s < C\Uq\h s , provided 
Ai(-) e C s and s > [d/2] +2. 



Proof. Integration by parts together with (|1.1U|) yields the basic I? estimate 



±(U,A U) t = ~(U,A o t U) + (U,A°U t ) 

= 1 -(U,A° t U)-(U,Y J tiU Xj ) 
(1.17) i 

= l(U,(A t +^2Ai ] )U) 

3 

< C\U\^2 |t7|jyl,oo < C|C/|^2 |f^|_H" s ■ 

Here, we have used (|TT2l in equating (U,^ A^U Xj ) = -f<E7,£j A J Xj U), 
original equation Ut = —J2jA J U x . in estimating A® < C\Ut\ < C2 1 C^c I i n 
the second-to-last inequality, and Ijl.lMjl in the final inequality. 



S 



A similar, higher-derivative calculation yields the H s estimate 
(1-18) s 

\ct{dl,A°dlU)) t = ±J2(d r x U,A°d r x U) +Y j {dlXJ,A a dlU t ) 

r=0 r=0 r=0 

s s 

= ^{%u,J!S%v) - ^2(d r x u,^2A d r x A^u x .) 

r=0 r=0 j 

r=0 j 
s r 

r=0 £=1 

<C|C/|^(|c%. + (|[/[^)> 

where the final inequality follows by (|1.15|) : see Exercise 1 1 . 81 b elow . 
Defining 

1 i 

(1.19) C(t):=-(Y,(O r x ,A°d r x U)), 

r=0 

we have, therefore, the Ricatti-type inequality 
(1-20) (t<C{C + ( s ), 

yielding ((t) < C£(0) for small t, provided £(0) is sufficiently small. Ob- 
serving that C 1/2 is a norm equivalent to |?7[#£, we obtain bounded local 
stability, provided a solution exists. 

Essentially the same a priori estimate can be used to show existence and 
uniqueness of solutions. Define the standard nonlinear iteration scheme (see, 
e.g., [Fr, Ma.3]) 

(1.21) A°(U n )U? +1 + ^ A- j {U n )TJ^ x = 0, U(0) = U . 

j 

For U n G H s , an H s solution TU n+1 of (fOTj) may be obtained by lin- 
ear theory; see Remark 1.3.81 Section 3.2. By the estimate already obtained, 
we find, that T takes the ball B := {U : [t^|L°°([o,r];_ff a (x)) — %\Uo\h s (x)} 
to itself, for r > sufficiently small. A similar energy estimate on the 
variation e := T(Ui) — T(U2), for Uj G B yields that T is contractive in 
L°°([0, t]; L 2 {x)) on the invariant set B, and stable in L°°[0, r]; H s (x), yield- 
ing existence of a fixed-point solution U G L°°([0, T]; H s (x)) (Exercise II. 9|1 . 
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Likewise, uniqueness of solutions may be obtained by a stability estimate 
on the nonlinear variation e := U\ — U2, where U\ and U2 denote solutions 

of unj). □ 



Remarks 1.7. 1. Clearly, we do not obtain global well-posedness by this 
argument, since solutions of Ricatti-type equations in general blow up in 
finite time. Indeed, it is well-known that shock-type discontinuities may 
form in finite time even for arbitrarily smooth initial data, corresponding to 
blow-up in H ] see, e.g., [La. 1-2, J, KIM, Si]. 
2. Using the strong Moser inequality 

r 

(1.22) |(9 Qi Vl )---(a Q ^ r )| i2(IRd) <c^kl^ ( ^)(IIl^l^(K d )) 

i=l j^i 

for s = \ctj\ (proved using Gagliardo-Nirenberg inequalities [T]), the same 
argument may be used to show that smooth continuation of the solution is 
possible so long as It/l^/i.oo remains bounded; see [Ma.3], Chapter 2. 

Exercise 1.8. Verify the final inequality in (|1.18j) by showing that 

(1.23) \d r ~ s A(U)U x \<C ni<3<r+i-s\dx j \. 

£[aj|<r+l-s 

Exercise 1.9. If u n £ H s {x) are uniformly bounded in H s , and convergent 
in -£/^(x), show that lim n — >C>Q u n 6 H s , using the definition of H s as the set of 
v € L 2 such that, for all 1 < r < s, (Vjd^tp) ^ C\<P\l 2 f or all test functions 
(j) G Cq°, together with the fact that limits and distributional derivatives 
commute. 

Symmetrizability is at first sight a rather restrictive requirement in more 
than one spatial dimension. However, it turns out to be satisfied in many 
physically interesting situations, in particular for gas dynamics and MHD. 
Indeed, a fundamental observation of Godunov is that symmetrizability is 
closely related with existence of an associated convex entropy. 

Definition 1.10. A hyperbolic entropy, entropy flux ensemble is a set of 
scalar functions [r]^ 3 ) such that 

(1.24) dndF j = dq j , 
or equivalently 

(1-25) n(U) t + J2qi(U) Xj =0 

j 

for any smooth solution U of (jl.lj) . 
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Proposition 1.11 ([God, Mo, B, KSh]). For U lying in a convex set 
U, existence of a convex entropy n is equivalent to symmetrizability of (jl.lj) 
by an invertible coordinate change U — > W := dr/ (known as an "entropy 
variable"), i.e., writing 

(1.26) U(W) t + J2 Fj (W) Xj =0, 

j 

we have A := (dll/dW) = (cPr?)" 1 symmetric positive definite and A? : = 
(dFi /dW) = A> A symmetric. 

Proof. (Exercise) (=^») Differentiate (|1.24|) and use symmetry of d 2 r\. {<=■) 
Reverse the calculation to obtain (|1.24|) with drj := W, then note that dr\ is 
exact, due to symmetry of dW := A . □ 

Remarks 1.12. 1. Symmetrizability by coordinate change implies sym- 
metrizability in the usual quasilinear sense (exercise), but not the converse. 
In particular, (nonlinear) symmetrization by coordinate change preserves 
divergence form, whereas (quasilinear) symmetrization by a left-multiplier 
A does not, cf. iPfll and (Ol) . 

2. Symmetrizability by coordinate change implies also (|1.25p . which 
yields the additional information that J r/dx, without loss of generality equiv- 
alent to the L 2 norm of U, is conserved for smooth solutions. Thus, we find 
in the discussion of Remark 11.71 that blowup occurs in a derivative of U and 
not in U itself. 

3. For gas dynamics and MHD, there exists a convex entropy in the 
neighborhood of any thermodynamically stable state, namely the negative 
of the thermodynamical entropy s; see, e.g. [Kaw, MaZ.4, Z.4]. In particular, 
for an ideal gas, there exists a global convex entropy. 

4. For the hyperbolic shock stability problem, hypotheses of hyperbol- 
icity, symmetrizability, etc., are relevant only in neighborhoods of U± and 
not between. Thus, a shock may be stable even if U± are entirely separated 
by unstable constant states, as, e.g., for phase-transitional shocks in van der 
Waals gas dynamics; see [Pre. 3, B-G.2-3]. 

1.1.2. Viscous equations. Analogous to (jl.l()|) in the setting of the 
viscous equations (jl.2j) is the symmetric hyperbolic-parabolic form 

(1.27) A°W t + &W Xj = J2(& k W Xk ) Xj + 

3 3,k 




11 



A symmetric positive definite, A 3 U symmetric, B^ k 
with 



block-diag {0, b jk 



} 



(i.28) 5>^' fc >#iei 2 , e>o, 



for all ^eR d , and 



(1.29) 



G 



(. 



~g{d x W) 







) 



with g = 0(\W X \ ), to be achieved by an invertible coordinate change U — > 
14 7 combined with left-multiplication by an invertible lower block-triangular 
matrix S(U); see [Kaw, KSh] and references therein, or Appendix Al, [Z.4]. 

Remark 1.13. As pointed out in the references, ^4° may be taken without 
loss of generality to be block-diagonal, thus identifying "hyperbolic" and 
"parabolic" variables w 1 and w 11 . Indeed, w 1 may be taken without loss of 
generality as u 1 and w 11 as any variable satisfying the (clearly necessary) 
integrability condition B^ k U Xk = /# fc (l7)wg for all j, k [GMWZ.4, Z.4]. 

Similarly as in the hyperbolic case, we may deduce local well-posedness 
of (|1.27|) directly from the structure of the equations, using energy esti- 
mates/integration by parts. Here, we shall require also a standard but es- 
sential tool for multidimensional parabolic systems, the Garding inequality 



for Lipshitz W G R nxn satisfying uniform ellipticity condition (|1.28|) . and 
< 6 < 9, with C = C(6, \d x V k \ L oo) < C 2 \d x & k \ L o a /\0-6\ for some uniform 
C 2 > 0. 

Exercise 1.14. (i) Prove (jl.3U|) with C = in the case ty k = constant, 
using the Fourier transform and ParsevaVs identity, (ii) Prove (|1,.30|) with 
C = in the case that W k varies by less that \6 — 9\/d 2 from some constant 
value, i.e, oscillation V < 2\8 — 9\/d 2 . (in) Prove the general case using a 
partition of unity {xr} such that oscillation ti k < 2\9 — 9\/d 2 on the support 
of each \r, and the estimates 



Y,^J^ k d x J) =Y,12( d *iXrffi k d x Xrf) + 0{\d x f\ L *\(d x Xr)f\L*) 



(1.30) 



Y,{d Xj f:V k d x J) > ~9\d x f\\ 2 - C\f\\ 2 



(1.31) 



r j,k 



12 



and 

(1.32) \8 x f\ L 2 = \Xrd x f\ L 2 = (\9xXrf\v + \(9xXr)f\v>: 

r r 

Remark 1.15. The Garding inequality (jl.30f) is an elementary example of 
a pseudodifferential estimate. Pseudodifferential techniques play a funda- 
mental role in the analysis of the curved shock problem; see [Ma. 1-3, Me. 4], 
[GMWZ.l, GMWZ.3-4, W], and references therein. 

Proposition 1.16 ([Kaw]). Symmetric form (jl.27|) implies local well- 
posedness in H s , and bounded local stability \U(t)\ji s < C\Uq\h s , provided 
A>(-), & k {-) G C s and s > [d/2] + 2. 

Proof. For simplicity, take g = 0; the general case is similar. Similarly as in 
11.171 we have 

^{W,A°W) t = ^{W,A° t W) + {W,A°W t ) 



(1.33) 



\{w,a*w) - nr.^.i'ir,. - J2(& k w Xk ) Xj ) 
hwj^ + iw 1 ,^^.™ 1 } + 0(\d x w n \\W\ 



J 



< C\W\ 2 L2 \W 

< c\w\ 2 L2 \w\ Hld/2]+2 , 

where in the second-to-last inequality we have used (|1.3()|) with 9 = 9/2, 
together with Young's inequality \d x w n \\W\ < (l/2C)\d x w n \ 2 + (C/2)\W\ 2 
with C > sufficiently large, and used the original equation to bound \ Wt\ < 
C(|Wa;| + \d x wII \)> an d in the final inequality we have used the Sobolev 
inequality l|1.13j) . 

Likewise, we obtain by a similar calculation 



(1.34) - } 



l(J2{d r x ,A°d r x W)) t < C\W\ 2 H s(\W\ HW2]+2 + \W\ s Hld/2]+2 ) 



r=0 



<C(\W\ 3 HS + \W\ S + 1 ), 



by the Moser inequality ()1.15j) and s > [d/2] +2, from which the result 
follows by the same argument as in the proof of Proposition 11.61 □ 
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The following beautiful results of Kawashima et al (see [Kaw, KSh] and 
references therein) generalize to the viscous case the observations of Go- 
dunov et al regarding entropy and the structure of the inviscid equations. 

Definition 1.17 ([Kaw, Sm]). A viscosity-compatible convex entropy, en- 
tropy flux ensemble (rj, q 3 ) for (|1.7|) is a convex hyperbolic entropy, entropy 
flux ensemble such that 

(1.35) d 2 V J2 &k &k>0 
and 

(1.36) rank $ld 2 r) ^ = rank ^ = r 

j,k j,k 

for all nonzero £ G 

Exercise 1.18. Using (jl,8|) - ()1.9|) . and the fact that 

(1.37) (d^^B^id 2 ^- 1 / 2 = {d 2 nr 1 ' 2 {d 2 r ] B^){d 2 r 1 )- 1 l 2 

is similar to , show that d 2 r\B^ k symmetric is sufficient for (|1.35|) - ()1.36|) . 

Proposition 1.19 ([KSh, Yo] 2 ). A necessary and sufficient condition 
that a system ([1.2ft , (|1.9jl can be put in symmetric hyperbolic-parabolic form 
()1.27|) with G = and A> symmetric (A not necessarily block-diagonal) by 
a change of coordinates U — > W{U) for U in a convex set U, is existence 
of a convex viscosity-compatible entropy, entropy flux ensemble r\, q> defined 
on U, with W = dr](U). 

Proof. Identical with that of Proposition II . 1 ll as concerns ^4° and A J . Re- 
garding & k = B^ k (d 2 r])~ 1 = (d 2 r])~ 1 (d 2 r]B : ' k )(d 2 T])~ 1 , we have that condi- 
tions Yljk^ ' £j£k > anci rank^J-B- 3 = r f° r an nonzero £ G R d are 
equivalent both to (OSl - lfOn]! and to (TT~28|l and B jk = block-diag {0, b jk } 
(exercise). The latter assertion depends on the observation that ^ ■ k B^^j^ > 

together with structure (jl.Hjl implies that 3? k & £j£k is block-diagonal 
and vanishing in the first diagonal block. This completes the argument. □ 

Established for d 2 r)B jk symmetric in [KSh]. The observation that 11. 351 - 111. 3 (ill suffice 
seems to be due to Yong [Yo]. 
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Definition 1.20 ([Kaw]). Augmenting the local conditions (jl.27JI - (jl.28|) . 

we identify the time- asymptotic stability conditions 3 of "first- order symme- 
try, " A 3 symmetric, and "genuine coupling: " 

(GC) No eigenvector of ^ ^dF 3 lies in ker ^ £j£ k B jk , for ^ G R d . 

j 

Proposition 1.21 ([Kaw] 4 ). Symmetric hyperbolic-parabolic form (|1,27[) 
together with the time- asymptotic stability conditions of first- order symme- 
try, A 3 symmetric, and genuine coupling, (GC), implies L 1 n H s — > H s 
time- asymptotic stability of constant solutions, U = constant, provided that 
A?, & k G C s and s > [d/2] + 2, with rate of decay 

(1.38) \U -U\ H s(t) <C(l + t)-i\U -U\ H s nLl (0) 

for \U — U\H a nL 1 (0) sufficiently small, equal to that of a d-dimensional heat 
kernel. 

Proposition 1.22 ([Kaw] 5 ). Existence of a convex viscosity- compatible 
entropy, together with genuine coupling, (GC), implies global well-posedness 
of (|1.7|) - (jl.9j) in H s , and bounded stability of constant solutions, 

(1.39) \U-U\ H s(t) < C\U-U\ H s(Q), 

provided A 3 , B jk G C s and s > [d/2] + 2. 

Propositions 11.211 and 11.221 depend on a circle of ideas associated with 
the phenomenon of "hyperbolic-parabolic smoothing", as indicated by the 
existence of parabolic-type energy estimates 
(1-40) 

\W(t)\ 2 H s + J Q (idxWll^ + ld^llsys < c(\W(0)\ 2 H s+ \W{s)\ 2 L2 ds) 

for s sufficiently large. We defer discussion of these important concepts to 
the detailed treatment of Sections 3 and 4. 

Exercise 1.23. Result (|1.39|) was established for gas dynamics in the sem- 
inal work of Matsumura and Nishida [MNi]. The key point is existence of 
an L 2 energy estimate 

(i.4i) \w(t)\h+ fid^n^ciwmh, 

Jo 

3 This refers to stability of constant solutions; see Proposition 1 1 . 2 ll below. 
4 Established in [Kaw] for d > 3, or d > 1 in the case of a convex entropy. 
5 Established in [Kaw] for "physical" convex entropies satisfying p. 43^ below. 
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which, coupled with the general machinery used to obtain (|1,4U|) . yields an 
improved version of (|1,4U|) in which the J \U(s)\ 2 ^ 2 ds term on the righthand 
side does not appear. This in turn implies (jl.39|) . Prove (f 1 .411) in the general 
case using the viscous version 

(1.42) V (U) t + Q j (U) Xj = J2(dvB jk U Xk ) Xj - £ U^nB^ 

j j,k j,k 

of entropy equation (|1.25|) and Exercise \l.lJ^[ ii) , under the assumption that 
\W\h s (hence \W\l°° ) remains sufficiently small, for s sufficiently large. 

The Gdrding inequality is not required for gas dynamics or MHD, for 
which there hold the strengthened ellipticity condition 

(1-43) Y,^P d2r l B3k fx k )> e \ dx f\ 2 - 

Under this assumption, show using (jl.42|) that (|1,39|) holds for s = 0, inde- 
pendent of the size of \W(0)\i2 . 

Remark 1.24. Similarly as in the inviscid case, the negative of the ther- 
modynamical entropy serves as a convex viscosity-compatible entropy for 
gas and MHD in the neighborhood of any thermodynamically stable state; 
see [Kaw, MaZ.4, Z.4]. In particular, for an ideal gas, there exists a global 
convex viscosity-compatible entropy. 

Example 1.25. The Navier-Stokes equations of compressible gas dynam- 
ics, may be written as 

p t + div(pu) =0, 

divr 

(1.44) 

{pu)t + div(pu <g> u) + Vp = pAu + (A + ^)Vdivu, 
(p(e + 7^ 2 ))t + div(p(e + ^u 2 )u + pu) = div(r • u) + kAT, 

where p > denotes density, u € M. d fluid velocity, T > tempera- 
ture, e = e(p,T) internal energy, and p = p(p,T) pressure. Here, r := 
Adiv(u)/ + 2pDu, where Dujk = ^(u J Xj + ui k ) is the deformation tensor, 
X(p,T) > and p(p,T) > are viscosity coeffients, and k(p,T) > is the 
coefficient of thermal conductivity. The thermodynamic entropy s is defined 
implicitly by the underlying thermodynamic relations e = e(v,s), T = e s , 
p = —e v , where v = p^ 1 is specific volume. (These may or may not be inte- 
grable for arbitrary choices of p(-, •), e(-, •); on the other hand, each choice 
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of e satisfying e s = T > uniquely determines functions p and e.) Thermo- 
dynamic stability is the condition that e be a convex function of (v,s), or 
equivalently (exercise) ex, Pp > 0. The Euler equations of compressible gas 
dynamics are (|1.44j) with the righthand side set to zero. 

Evidently, (|1.44|) is of form Q1.2|) . (|1.8|) . with "conservative variables" 
U = (/?, pu, p(e+^u 2 )) . By Remark |l.l31 we have that symmetric hyperbolic- 
parabolic form (|1.27j) . if it exists, may be expressed in the "natural variables" 
W = (p,u,T). Indeed, this can be done [KSh], with 

fPp/p \ 
(1.45) 1° = pl d , 

\ pe T /Tj 



({p P )u-i p P t \ 

(1.46) j2 A ^j=[ v P t l p(u-t)h p T e , 

3 \ p T i (pe T /T)u ■ £J 

and 

/0 \ 

(1.47) &k tek = p\i\ 2 I d + {p + 

3,k \o o r- 1 ^! 2 / 

whenever p„, > (thermodynamic stability), in which case we also have 
the time-asymptotic stability conditions of symmetry of A 3 and genuine 
coupling (GC) (by inspection, equivalent to p p ^ 0). If only ex > 0, as for 
example for a van der Waals-type equation of state, we may still achieve sym- 
metric hyperbolic-parabolic form, but without symmetry of A 3 or genuine 
coupling, by dividing through p p from the first equation (~ first row of each 
matrix 1°, A> , & k ). That is, we obtain in this case local well-posedness, 
but not asymptotic stability of constant states. Similar considerations hold 
in the case of MHD; see [Kaw, KSh]. 
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2 Description of results 



We now turn to the long-time stability of viscous shock waves. Fixing the 
viscosity coefficient u, consider a viscous shock solution (|1.4j) of a system of 
viscous conservation laws (jl.2|) of form (jl.9|) . Changing coordinates 

if necessary to a rest frame moving with the speed of the shock, we may 
without loss of generality arrange that shock speed s vanish, so that (|1.4[) 
becomes a stationary, or standing-wave solution. Hereafter, we take v = 1 
and s = 0, and suppress the parameters v and s. 

2.1. Assumptions. The classical Propositions 11.211 and 11.221 con- 
cern long-time stability of a single, thermodynamically stable equilibrium, 
whereas a viscous shock solution typically consists of two different ther- 
modynamically stable equilibria connected by a profile that in general may 
pass through regions of thermodynamical instability; see Remark ll.121 4. and 
examples, [MaZ.4, Z.4]. Accordingly, we make the following structural as- 
sumptions, imposing stability at the endstates, but only local well-posedness 
along the profile. 



Assumptions 2.1. 

(Al) For U in a neighborhood of profile U(-), there is an invertible 
change of coordinates U — > W such that ([1.2ft may be expressed in W coor- 
dinates in symmetric hyperbolic-parabolic form, i. e, in form ()1.27|) with A 
symmetric positive definite and (without loss of generality) block- diagonal 
and symmetric, & k = block-diag {0,^} with V satisfying uniform 

ellipticity condition (jl,28j) . and G = ft\ with g = 0(\dxW\ 2 ). 

(A2) At endstates U±, the coefficients A> , W k defined in (Al) satisfy 
the constant- coefficient stability conditions of first- order symmetry, A± sym- 
metric, and genuine coupling, (GC). 

Remark 2.2. Conditions (Al)-(A2) are satisfied for gas- and magnetohy- 
drodynamical shock profiles connecting thermodynamically stable endstates, 
under the mild assumption er > on the equation of state e = e(p, T) re- 
lating internal energy e to density p and temperature T: in particular for 
both ideal and van der Waals-type equation of state. (Exercise: verify this 
assertion for gas dynamics, starting from form (jl.45j) - ([1.47j) and modifying 
all coeffients M by the transformation 

(2.1) M -» X M + (1 - X )block-diag {p~\ I d , 1}M, 
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where x = x{W) 1S a smooth cutoff function supported on the region of 
thermodynamic stability p p > 0.) 

To Assumptions l2.il we add the following technical hypotheses. 

Assumptions 2.3. 

(HO) F 3 , B 3k , W(-), A G C q+1 , q > q(d) := [d/2] + 2. 

(HI) A^iJJ) is uniformly definite, without loss of generality A\^ > 8 > 0. 
(For necessary conditions, detyl} 1 ^ is sufficient.) 

(H2) det^4± / 0, or equivalently detA]_ ^ 0. 

(H3) Local to U, the solutions of (|1.4I) form a smooth manifold {U S }, 
5 £ U C R*, with U° = U. 

Condition (HO) gives the regularity needed for our analysis. Condition 
(H2) is the standard inviscid requirement that the hyperbolic shock triple 
(U-,U+, 0) be noncharacteristic. Condition (HI) states that convection in 
the reduced, hyperbolic part of (|1.27|) governing coordinate w 1 is either up- 
or downwind, in particular uniformly noncharacteristic, everywhere along 
the profile; as discussed in [MaZ.4, Z.4], this is satisfied for all gas-dynamical 
shocks and at least generically for magnetohydro dynamical shocks. Condi- 
tion (H3) is a weak form of (implied by but not implying) transversality of 
the connection U as a solution of the associated traveling-wave ODE. 

These hypotheses suffice for the investigation of necessary conditions for 
stability. In our investigation of sufficient conditions, we shall require two 
further hypotheses at the level of the inviscid stability problem. 

Assumption 2.4. 

(H4) The eigenvalues o/^ • have constant multiplicity with respect 
to £ e R d , 0. 

Denote by a^(£), j = 1, . . . ,n the eigenvalues of ^ ■ -A±£j> necessarily 
real by (A2), indexed by increasing order. These are positive homogeneous 
degree one and, by (H4), locally analytic on £ G M. d \ {0}. 6 Here, and 
elsewhere, real homogeneity refers to homogeneity with respect to positive 
reals. Let 

(2.2) P±(Z,t) :=ir + J2iCjAi 

3 

6 In fact, they are globally analytic, since they maintain fixed order; however, this is 
unimportant in the analysis. 
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denote the frozen-coefficient symbols associated with (|1.6(l at U = U±. Then, 
detP±(£,r) = has n locally analytic, positive homogeneous degree one 
roots 

(2.3) ir = -iaf(£), r = l,...,n, 

describing the dispersion relations for the frozen-coefficient initial- value prob- 
lem (IVP). The corresponding objects for the initial-boundary- value prob- 
lem (IBVP) are relations 

(2.4) «£l = Mr(f,i"), r = l,...,n, 

£ = (£i, ■ ■ ■ ,£d) ='■ (£i>Oi describing roots of det Q±(r) = 0, where 

d 

(2.5) Q±(£,t) :=i£ 1 + (^)- 1 (ir + ^^ ± ) 

denote the frozen-coefficient initial-boundary-value symbols associated with 
ljl.6|) at J7 = Z7±. Evidently, graphs (|2.3j) and (j2.4|) describe the same sets, 
since det A 1 det Q-t (r) = det P-j- and det ^43j_ 7^ 0; the roots ir describe char- 
acteristic rates of temporal decay, whereas \x = i£i describe characteristic 
rates of spatial decay in the x\ direction. 

Definition 2.5. Setting £ = (£i,...,£<i) =: we define the glancing 

sets G(P±) as the set of all (£, r) smc/i that, for some real £1 and 1 < r < n, 
t = — a^(£i,£) and (da^/d£i)(£i,£) = 0: i/taf is, the projection onto (£,t) 
of the set of real roots (£,t) of det P± = at which (|2.3|) is not analytically 
invertible as a function (|2.4|) . T/ie roote (£, r) are called glancing points. 

Assumption 2.6. 

(H5) Each glancing set G(P±) is the (possibly intersecting) union of 
finitely many smooth curves r = rfr (£), on which the root £1 o/ir+a^(-,£) = 
aas constant multiplicity s q (by definition >2). 

Among other useful properties (see, e.g., [Me. 3, Z.3, MeZ.2]), condition 
(H4) implies the standard block structure condition [Ma. 1-3, Me. 2] of the 
inviscid stability theory [Me. 3]. Condition (H5), introduced in [Z.3], imposes 
a further, laminar structure on the glancing set that is convenient for the 
viscous stability analysis. In either context, glancing is the fundamental 
obstacle in obtaining resolvent estimates; see, e.g., [K, Ma. 1-3, Z.3, MeZ.l, 
GMWZ.1-4]. The term "glancing" derives from the fact that, at glancing 
points (£,t), null bicharacteristics of detP± lie parallel to the shock front 
xx = 0. 
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Remarks 2.7. 1. Condition (H5) holds automatically in dimensions one 
(vacuous) and two (trivial), and also for the case that all characteristics 
a r (£) are either linear or convex/concave in £1. (Exercise [GMWZ.2]: prove 
the latter assertion using the Implicit Function Theorem.) In particular, 
(H5) holds for both gas dynamics and MHD. 

2. Condition (H4) holds always for gas dynamics in any dimension, 
and generically for MHD in one dimension, but fails always for MHD in 
dimension greater than or equal to two. For a refined treatment applying 
also to multidimensional MHD, see [MeZ.3]. 

Assumption 2.8. For definiteness, consider a classical, "pure" Lax p- 
shock, satisfying 

(L) > s = > a+, 

where := a- (1,0, . . . , 0) denote the eigenvalues of A± := dF 1 (U±), and 
I = 1 in (H3). See, e.g., [ZS, Z.S~4] for discussions of the general case, 
including the interesting situation of nonclassical over- or undercompressive 
shocks. 

2.2. Classical stability conditions. Recall the classical admissibility 
conditions described in the introduction of: 

1. Structural stability, defined as existence of a transverse viscous profile. 
For gas dynamics, this is equivalent to the (hyperbolic) Liu-Oleinik admis- 
sibility condition; see [L.2, Gi, MeP]. In general, it may be a complicated 
ODE problem involving both and the specific form of the viscosity B^ k . 

2. Dynamical stability, defined as local hyperbolic well-posedness, or 
bounded, bounded-time stability of ideal shock ()1.3|) . Following Majda 
[Ma. 1-3], we define weak dynamical stability (Majda's Lopatinski condition) 
as the absence of unstable spectrum 3?A > for the linearized operator (ap- 
propriately defined) about the shock, an evident necessary condition for sta- 
bility. (Recall, equations (jl.lj) are positive homogeneous degree one, hence 
spectrum lies on rays through the origin and instabilities, should they occur, 
occur to all orders.) 

In the present context of a Lax p-shock, under assumptions (A1)-(A2) 
and (H2), the necessary condition of weak dynamical stability may be ex- 
pressed in terms of the Lopatinski determinant 

(2.6) A(|, A) := det (r^, ■ ■ ■ , rj,r+. 1 , • • • ,r+, \[U) + i[F*]) 
as 

(2.7) A(£, A) ± for £ £ R^ 1 and 5?A > 0, 



21 



where = (£, A) denote bases for the unstable (resp. stable) subspace 
of the coefficient matrix 

(2.8) A±:={A 1 )- 1 {\ + iA^)± 

arising m the IBVP symbol Q± with A := ir, F% := V ;/ , /•''£,. and 
A^ := Ylj^i-^ij- Functions Tj , and thus A, are well-defined on 5RA > 
and may be chosen analytically in (£,A), by a standard lemma of Hersch 
[H] asserting that A± has no center subspace on this domain; see exer- 
cise I2.1UI below. Zeroes of A correspond to "normal modes" , or solutions 

{U, X) = e xt e l ^' x (U(xi), X) of the linearized perturbation equations, where 
xi = X(x, t) denotes location of the shock surface. We define strong dynam- 
ical instability as failure of (|2.7|) . 

Under assumptions (A1)-(A2), (H2), and (H4), r- , and thus A, may 
be extended continuously to the boundary KA = 0; see [Me. 3, CP]. Strong 
dynamical stability (Majda's uniform Lopatinski condition) is then defined 
as 

(2.9) A(£,A)/0 for £ 6 R d_1 , !RA > 0, and (£, A) ^ (0,0). 

Strong dynamical stability under our assumptions is sufficient for dynamical 
stability, by the celebrated result of Majda [Ma. 1-3] together with the result 
of Metivier [Me. 3] that (H4) implies Majda's block structure condition. 

Remarks 2.9. 1. In the one-dimensional case, £ = 0, both (|2.7|) and H2.9|) 

reduce to nonvanishing of the Liu-Majda determinant 

(2.10) * := det (rjy ■ ■ ,r~,r+ +1 ,--- ,r+ [[/]), 

where denote eigenvectors of A\. associated with characteristic modes 
outgoing from the shock. As a single condition, this is generically satisfied, 
whence we find that shocks of Lax type (L) are generically stable in one 
dimension from the hyperbolic point of view. 

2. As pointed out by Majda [Ma.1-3], for gas dynamics the region of 
strong dynamical stability is typically separated from the region of strong 
dynamical instability by an open set in parameter space a := (U-,U+, s) 7 
of indeterminate stability, on which A(£, A) has roots A lying precisely on 
the imaginary axis, but none with KA > 0; thus, the point of transition from 

7 More precisely, the subset of a corresponding to Lax p-shocks on which A was defined, 
or an open neighborhood with A extended by 12.611 . 
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stability to instability is not determined in the hyperbolic stability theory. A 
similar situation holds in MHD [Bl, BT.1-4, BTM.1-2]. The reason for this 
at first puzzling phenomenon, as described in [BRSZ, Z.4], is that, under 
appropriate normalization, A(£, •) takes imaginary A = ir to imaginary 
A for \t\ sufficiently large relative to |£|: more precisely, for (£,t) lying 
in the "hyperbolic regions" TL± bounded by glancing sets Q(P±). Thus, 
fixing £, and explicitly noting the dependence on parameters a, we find that 
imaginary zeroes A = vr of A Q (£, ■) for which (£, r) £ 7i±, such as occur for 
gas dynamics and MHD, generically persist under perturbation in a, rather 
than moving into the stable or unstable complex half-planes 5RA < and 
KA > as would otherwise be expected; see exercise 12 . 1 1 1 below . 

Exercise 2.10. Under (Al)-(A2), (H2), show that det (A±{£, A)-i£i) = 
implies detQ^ i£jA 3 +\) = 0, A± defined as in (|2.8|) . violating hyperbolicity 
of if 3ftA 0, = 0. Thus, A± has no center subspace on the 

domain 3ft A > 0. It follows by standard matrix perturbation theory [Kat] 
that the stable and unstable subspaces of A± have constant dimension on 
9ftA > 0, and the associated A±-invariant projections are analytic in (£,A). 
This implies the existence of analytic bases r^, by a lemma of Kato asserting 
that analytic projections induce analytic bases on simply connected domains; 
see [Kat], pp. 99-102. 

Exercise 2.11. Let /(a, A) be continuous in a and A £ C, with the addi- 
tional property that /(a, •) : R — ► R. (a) Show that in the vicinity of any root 
(0, Ao), Ao real, such that /a(0, Aq) exists and is nonzero, there is a continu- 
ous family of roots (a,A(a)) with A(a) real (Intermediate Value Theorem). 
More generally, odd topological degree of /(0, •) at A = Ao is sufficient; for 
even degrees f(a, A) = (a + A 2 ) m , a € R, Ao = is a counterexample, (b) 
If f is analytic in \, show that /(a, A) = /(a, A), where z denotes complex 
conjugate of z, hence nonreal roots occur in conjugate pairs. 

2.3. Viscous stability conditions. We'll both augment and refine 
the classical stability conditions through a viscous stability analysis, at the 
same time providing their rigorous justification. 

2.3.1. Spectral stability. We begin by adding to the conditions of 
structural and dynamical stability a third condition of (viscous) spectral 
stability. Let L denote the linearized operator about the wave, i.e., the 
generator of the linear equations Ut = LU obtained by linearizing (|1.2|) about 
the profile U, and the family of operators obtained from L by Fourier 

transform in the directions x parallel to the front, indexed by frequency £. 
Explicit representations of L and are given in Section 3.2. 



23 



Definition 2.12. Similarly as in the hyperbolic stability theory, we define 
weak spectral stability (clearly necessary for viscous stability) as the absence 
of unstable L 2 spectrum 9?A > for the linearized operator L about the 
wave, or equivalently 

(2.11) A G" a(L~) for i G R^ 1 and 5RA > 0. 

We define strong spectral stability (neither necessary nor sufficient for viscous 
stability) as 

(2.12) A G" a(Lg) for f G R d_1 , 5?A > 0, and (£, A) ^ (0,0). 
We define strong spectral instability as failure of Q2.11JI . 

Conditions ()2.11|) and 1)2.12(1 may equivalently be expressed in terms 
of the Evans function D(^,X) (defined Section 5), a spectral determinant 
analogous to the Lopatinski determinant of the hyperbolic C£IS6 ? ELS 

(2.13) D(i, A) + for £ G R^ 1 and 3?A > 
and 

(2.14) D(i, A) / for f G M d_1 , 5?A > 0, and (£, A) ^ (0, 0), 

respectively, where zeroes of D correspond to normal modes, or solutions 

U = e^ t e l ^' x U(xi) of the linearized perturbation equations: equivalently, 
eigenvalues A of the operator obtained from L by Fourier transform in 
the directions x parallel to the front. Under assumptions (Al)-(A2), D may 
be chosen analytically in (£, A) on {f, A : KA > 0} \ {(0, 0)}; see [ZS, Z.3-4] 
or Section 5 below. 

Remarks 2.13. 1. It is readily verified that LU' = LqU' = 0, a conse- 
quence of translational invariance of the original equations (|1,2|) , from which 
we obtain G fj e ss (L) and in particular D(0,0) = 0. The exclusion of the 
origin in definition ([2.14(1 is therefore necessary for the application to stabil- 
ity of viscous shock profiles. If there held the stronger condition of uniform 
spectral stability, 

(2.15) 9ftr(L) < 0, 

we could conclude exponential stability of the linearized solution operator 
e LT by the generalized Hille-Yosida theorem (Proposition IA.111 Appendix 
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A) together with the high-frequency resolvent bounds we obtain later. How- 
ever, in the absence of a spectral gap between and Kct(L), even bounded 
linearized stability is a delicate question. Moreover, we here discuss non- 
linear asymptotic stability, for which we require a rate of linearized decay 
sufficient to carry out a nonlinear iteration. This is connected with the rate 
at which the spectrum of of moves into the stable complex half-plane 

5RA < as £ is varied about the origin, information that is encoded in the 
conditions of structural and refined dynamical stability defined just below. 

2. In the one-dimensional case £ = 0, strong spectral stability is roughly 
equivalent to linearized stability with respect to zero-mass initial perturba- 
tions; see [ZH, HuZ, MaZ.2-4]. Zero-mass stability is exactly the property 
that was used by Goodman and Xin [GoX] to rigorously justify the small- 
viscosity matched asymptotic expansion about small-amplitude shock waves 
in one dimension; indeed, their argument implies validity of matched asymp- 
totic expansion in one dimension for any structurally and dynamically stable 
wave about which the associated profiles exhibit zero-mass stability. 

2.3.2. Refined dynamical stability conditions. A rigorous connec- 
tion between viscous stability and the formal conditions of structural and 
dynamical stability is given by the following fundamental result (proved in 
the Lax case in Section 5). 



Proposition 2.14 ([ZS, Z.3, MeZ.2] 8 ). Given (A1)-(A2) and (H0)- 
(H4), and under appropriate normalizations of D and A, 



uniformly on £ £ R d_1 , > 0, for p := |£, A sufficiently small, where 7 
is a constant measuring transversality of connection U as a solution of the 
associated traveling-wave ODE, A is an appropriate Lopatinski condition, 
and I is as in (H3). In the present, Lax case, A is as in (J2.6|) and I = 1. 

Remark 2.15. In the absence of (H4), (|2.16|) still holds along individual 
rays through the origin, provided that they lie in directions of analyticity 
for A, in particular for IRA > 0, by the original argument of [ZS]. This 

8 Proved in [ZS] in the strictly hyperbolic case, along rays through the origin lying 
outside the glancing set; full version proved in [Z.3] under the additional hypothesis (H5), 
and in [MeZ.2] for the general case. 



o(\lM l ) 



(2.16) 
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is sufficient for the necessary stability conditions derived below; however, 
uniformity is essential for our sufficient stability conditions. 9 

Definition 2.16. We define refined weak dynamical stability as (|2.7j) aug- 
mented with the second-order condition 3i/?(£,ir) > for any real £, r 
such that A and D^' X (p) := D(p£,pX) are analytic at (£,ir) and (£,ir, 0), 
respectively, with A(£, ir) = and A,\(£, ir) ^ 0, where 
(2.17) 

0(i iT) . = (d/d P Y +1 D(plpir)\ p=0 = {d/dpY +1 D{plp\) 

(d/d\)A(lir) " (d/dX)(d/dpYD(plp\) ]p=ox=iT - 

Under (A1)-(A2), (H2), (H4), analyticity of A and D^ x (p) are equiv- 
alent, and hold for all real (£,r) lying outside the glancing sets Q{P±) [ZS, 
Z.3-4]. In this case, we define strong refined dynamical stability as (|2.7|) 
augmented with the conditions that: (i) 9ft/3(£, ir) > 0, A\(£,ir) 7^ 0, and 
{r^, . . . , r p _ l7 rp +1 , . . . , r+}(£ , ir) are independent, where are defined as 

in 1)2. 6|) (automatic for "extreme" shocks p = 1 or n), for any real (£, 
r) ^(-P±) such that A(£, ir) = and A is analytic, and (ii) the zero- 
level set in (£, r) of A(£, ir) = intersects the glancing sets Q± transversely 
in R^ -1 x IR (in particular, their intersections are trivial in dimension d = 2). 

Remark 2.17. The coefficient (3 may be recognized as a viscous correction 
measuring departure from homogeneity: note that (3 would vanish if D were 
homogeneous degree I. It has a physical interpretation as an "effective vis- 
cosity" governing transverse propagation of deformations in the front [Go. 2, 
GM, ZS, Z.3-4, HoZ.1-2]. 

2.3.3. Main results. With these definitions, we can now state our 
main results, to be established throughout the remainder of the article. 

Theorem 2.18 ([ZS, Z.3]). Given (A1)-(A2), (HO)-(HS), structural sta- 
bility 7^0, and one- dimensional inviscid stability A(0, 1) 7^ 0, weak spec- 
tral stability and weak refined dynamical stability are necessary for bounded 
Cq° — > L p , viscous stability in all dimensions d > 2, for any 1 < p < 00, 
and for profiles of any type. 

Theorem 2.19 ([Z.3-4]). Given (A1)-(A2), (H0)-(H5), strong spectral 
stability plus structural and strong refined dynamical stability are sufficient 

9 This point was not clearly stated in [Z.3]; indeed, l|2.1f)0 is hidden in the detailed 
estimates of the technical Section 4 of that reference, which depend on the additional 
hypothesis (H5). For a more satisfactory treatment, see [MeZ.2, GMWZ.3-4, Z.4]. 
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for asymptotic L C)H S — > H s nL p viscous stability of pure, Lax-type profiles 
in dimensions d > 3, for q(d) < s < q, q and q(d) as defined in (HO), and 
any 2 < p < oo, or linearized viscous stability in dimensions d > 2, for 
1 < s < q and 2 < p < oo, with rate of decay (in either case) 

(2.18) \U(t) - U\ H s < C(l + ^-^-^^[[/(O) - U\ L i nH . 

approximately equal to that of a (d — 1)- dimensional heat kernel, for any 
fixed e > and \U(0) — ^liinH 8 sufficiently small. Strong spectral plus 
structural and strong (inviscid, not refined) dynamical stability are sufficient 
for nonlinear stability in all dimensions d > 2, with rate of decay (|2.18j) . 
e = 0, exactly equal to that of a (d — 1)- dimensional heat kernel. Similar 
results hold for profiles of nonclassical, over- or undercompressive type, as 
described in [Z.3-4]. 

Remarks 2.20. 1. Rate 1)2.18(1 with e = is sharp for weakly inviscid stable 
shocks, as shown by the scalar case [Go. 3, GM, HoZ.1-2]. (As pointed out by 
Majda, scalar shocks are always weakly inviscid stable [Ma. 1-3, ZS, Z.3].) 

2. Given (A1)-(A2), (H0)-(H3), tA(0,1) ^ is necessary for one- 
dimensional viscous stability [ZH, MaZ.3] with respect to Cq° perturbations. 
(Note: C^( Xl ) £ C^(x).) 

3. Under (A1)-(A2), (H0)-(H3) plus the mild additional assumptions 
that the principal characteristic speed %>(£) be simple, genuinely nonlinear, 
and strictly convex (resp. concave) with respect to £ at £ = (£i,£) = 
(l,0rf_i), the sufficient stability conditions of Theorem 12 . 191 are satisfied for 
sufficiently small-amplitude shock profiles [FreS, PZ]. On the other hand, 
the necessary conditions of Theorem 12 . 1 81 are known to fail for certain large- 
amplitude shock profiles under (A1)-(A2), (H0)-(H3) [GZ, FreZ, ZS, Z.3- 
4]- 

4. The boundary between the necessary conditions of Theorem 12 . 1 81 and 
the sufficient conditions of Theorem 12.191 is generically of codimension one 
(exercise; see [Z.4], Remark 1.22.2). Thus, transition from viscous stability 
to instability is in principle determined, in contrast to the situation of the 
inviscid case (Remark 12. 9|) . 

5. It is an important open problem which of the conditions of refined 
dynamical and structural stability in practice determines the transition from 
viscous stability to instability as shock strength is varied; see [MaZ.4, Z.3-4] 
for further discussion. 
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3 Analytical preliminaries 



The rest of this article is devoted to the proof of Theorems 12.181 and 12.191 
We begin in this section by assembling some needed background results. 

3.1. Profile facts. Consider the standing wave ODE 
(F 1 ) / ([/) / = 0, 



U 1 , (F 1 ) 1 G M n " r , U n , (F 1 ) 11 G W, where '"" denotes d/dx and super- 
scripts I and II refer respectively to first- and second-block rows. Integrat- 
ing from — oo to x gives an implicit first-order system 

(3.2) (ByV = (F 1 ) // ([/) - (F 1 ) // (J7_) 
on the r-dimensional level set 

(3.3) (F 1 ) 1 (U) = (F 1 ) / (C/_). 

Lemma 3.1 ([MaZ.3]). Given (Al), the weak version detA^ ^ of 
(HI) is equivalent to the property that (j3.3|) determines a nondegenerate r- 
dimensional manifold on which (|3,2() determines a nondegenerate first-order 
ODE. Moreover, assuming (A1)-(A2) and det^^ 7^ 0, (H2) is equivalent 
to hyperbolicity of rest states U± . 

Remark 3.2. This result motivated the introduction of the weak form 
detA^ ^ of (HI) in [Z.3, MaZ.3]; the importance of the strong form 
A\ x > (< 0) was first pointed out in [MaZ.4]. 

Corollary 3.3 ([MaZ.3]). Given (A1)-(A2), and (H0)-(H2) 

(3.4) \d k x (U-U±)\<Ce-W, x^O, 
for < k < q + 2, where q is as defined in (HO). 

Proof of Lemma \3.1\ This result was proved in somewhat greater generality 
in [MaZ.3]. Here, we give a simpler proof making use of structure (Al)- 
(A2). From (Al), it can be shown (exercise; see [Z.4], appendix Al) that 
dll/dW is lower block-triangular and (B 1 ) 11 (dU/dW) = block-diag {0,6}, 
b nonsingular, with d(F 1 ) 1 /dw 7^ if and only if A^ 7^ 0, with A 1 as 
in (Al). Considering (|3,3|) - (|3,2|) with respect to coordinate W, we readily 
obtain the first assertion, with resulting ODE 

(3.5) b(w n )' = (F 1 ) II (U(w I (w II ),w 11 ) - (F 1 )^.), 
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where w I {w 11 ) is determined by the Implicit Function Theorem using rela- 
tion (F 1 ) 1 {U (W)) = (F 1 ) 7 ([/_). 

Linearizing (|3.5[) about rest point W = W± and using the Implicit Func- 
tion Theorem to calculate dw 1 /dw 11 , we obtain 

(3.6) (w 11 )' = M±w n := (S)- 1 (-a2iau«i2 + "2 2 )±w 7/ , 

where a± := (Si* 1 /e?PF)±, and b± are evaluated at W±. Noting that 
det(— ei2\OL^oi\2 + 0122) = det a/ detail, with an / as a consequence 
of detAj-L 7^ 0, we find that the coefficient matrix M± of the linearized 
ODE has zero eigenvalues if and only if deta±, or equivalently det A\. van- 
ishes. On the other hand, existence of nonzero pure imaginary eigenvalues 
i£i would imply, going back to the original equation (jl.27|) linearized about 
W±, that det(i^i,4 1 - (fB n )± = for some real £1 ^ 0. But, this is pre- 
cluded by (A1)-(A2) (nontrivial exercise; see (K3), Lemma 13.181 below). 
Thus, the coefficient matrix M has a center subspace if and only if (H2) 
fails, verifying the second assertion. □ 

Proof of Corollarv \3.t\ Standard ODE estimates. □ 

For general interest, and to show that our assumptions on the profile 
are not vacuous, we include without proof the following generalization of a 
result of [MP] in the strictly parabolic case, relating structure of a viscous 
profile to its hyperbolic type. 

Lemma 3.4 ([MaZ.3]). Let denote the dimension of the stable subspace 
of M + and d- the dimension of the unstable subspace of M_ , M± defined as 
in (13. Then, existence of a connecting profile U together with (Al)-(A2) 
and (H1)-(H2), implies that 

(3.7) I := d + + d- — r = i + + i_ — n, 

where i+ denotes the dimension of the stable subspace of and i- denotes 
the dimension of the unstable subspace of A l _; in the Lax case, 1 = 1. 

Remark 3.5. I = t in the case of a transverse profile, I as in (H3), whereas 
+ i_ — n = 1 may be recognized as the Lax characteristic condition (L). 
Assumption 12.81 could be rephrased in this case as t = 1 = i + + i_ — n. 

Proof. See [MaZ.3], Appendix Al, or [Z.4], Appendix A2. □ 
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3.2. Spectral resolution formulae. Consider the linearized equations 



(3 _ 8) U t = LU:=- Y^(A j U) Xj + Y,(B*U Xh ) Xj , U(0) = U , 

3 jk 

where 
(3.9) 

B jk ._ A j v .- dF j {U(x x ))v - (dB jl (U(xi))v)U'(xx). 

Since the coefficients depend only on x%, it is natural to Fourier-transform 
in x, x = (xi, . . . , Xd) ='■ (xi, x), to reduce to a family of partial differential 
equations (PDE) 

, " V 

U t = LfJ := (B U U')' - {A 1 ^)' -i Ai ZjU + * Yl Bil ^fi' 

(3.10) " j^i j^i 

+ i Y,(B lk tkU)' - Y, B3 %ZkU, U(0) = U 

in (x\,t) indexed by frequency £ € where '"" denotes d/dx\ and 

U = U(xi,£,t) denotes the Fourier transform of U = U(x,t). 

Finally, taking advantage of autonomy of the equations, we may take the 
Laplace transform in t to reduce, formally, to the resolvent equation 

(3.11) (A - L~)U = U , 

where U(xi,^, A) denotes the Laplace-Fourier transform of U = U(x,t) and 
Uq(xx) denotes the Fourier transform of initial data Uq(x). A system of ODE 
in x\ indexed by (£, A), (|3.11|) may be estimated sharply using either explicit 
representation formulae/variation of constants for systems of ODE or Kreiss 
symmetrizer techniques; see, e.g., [Z.3-4] and [GMWZ.2], respectively. The 
following proposition makes rigorous sense of this procedure. 

Proposition 3.6 ([MaZ.3, Z.4]). Given (Al), L generates a C° semi- 
group \e Lt \ < Ce 70 * on L 2 with domain V{L) := {U : U, LU £ L 2 } (Lu 
defined in the distributional sense), satisfying the generalized spectral reso- 
lution (inverse Laplace-Fourier transform) formula 

(3.12) e Lt f(x) = VN. / / e**+ A '(A-L^V>i,0<W 

J-y-ioo Jffi d -! 
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for 7 > 70, t > 0, and f G D{L), where f denotes Fourier transform of f . 
Likewise, 
(3.13) 



ft f"y-\-too r _ ^ 

e L ^f(s)ds = P.V. / e^ +A *(A-L £ ~)-V T (^i,<e,A)ci|dA 

for 7 > 7o, < t < T, and f G L 1 ([0, T]; P(L)), where g denotes Laplace- 
Fourier transform of g and 



(3.14) f T (x,s):-- 




for < s < T 
otherwise. 



Remark 3.7. Bound ()3.13|) concerns the inhomogeneous, zero- initial-data 
problem Ut — LU = /, U (0) = 0. Provided that L generates a C° semigroup 
e LT , the "mild", or semigroup solution of this equation is defined (see, e.g., 
[Pa], pp. 105-110) by Duhamel formula 



(3.15) U(t) = [ 

Jo 



t 

L(t-s) 



f(s)ds. 



This corresponds to a solution in the usual weak, or distributional sense, with 

somewhat stronger regularity in t. Formally, V := (A — L)^ 1 f T describes the 
Laplace-Fourier transform of the solution of the truncated-source problem 
V t — LV = f T , V(0) = 0. Evidently, U(s) = V{s) for < s < T. 

Proof. Let u n — > u and Lu n — ► /, u n G L)(L), u, f G L 2 , "— ►" denoting 
convergence in L 2 norm. Then, Lu n —r Lu by u n — > u, where Lu is defined 
in the distributional sense, but also Lu n —r f by Lu n — > /, whence Lu = f 
by uniqueness of weak limits, and therefore u G L)(L). Moreover, T>(L) is 
dense in L 2 , since Cq° C L)(L). Thus, L is a closed, densely defined operator 
on L 2 with domain T>(L); see Definition IA.1| Appendix A. 

By standard semigroup theory, therefore (Proposition IA.111 Appendix 
A), L generates a C°-semigroup \e LT \ < Ce 70 * if and only if it satisfies 
resolvent estimate 

(3-16) | (A _ L) -fc| i2 <^_ 

|A - 7oT 

for some uniform C > 0, for all real A > 70 and all k > 1, or equivalently 

(^7) \TT\ < C \( X - L ) kU \v 

(3 - 17) lUlL2 ~ |A- 7 o| fc 
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for all k > 1 and 

(3.18) range(A-L) = L 2 . 

In the present case, the latter condition is superfluous, since it is easily 
verified by results of Henry [He, BSZ] that A — L is Fredholm for sufficiently 
large real A. More generally, ()3.18[) may be verified by a corresponding 
estimate 

(3.19) \u\„ < WzE]El» 

I A — -yo I 

on the adjoint operator L* , equivalent to (|3.18j) given (|3.17|) (or, given (|3.18|) . 
to (|3.17l) . k = 1) and often available (as here) by the same techniques; see 
Corollary IA.17| Appendix A. 

A sufficient condition for ()3.17|) . and the standard means by which it is 
proved, is 

(3.20) \W\ L , < ^- L)S~y L , 

|A-7o| 

for all real A > 70, for some uniformly invertible transformation S: in this 
case the change of coordinates S : U — > W defined by 

(3.21) W := (&)V 2 (dW/dU){U(xi))U. 

Once (|3.17|) - (|3.18|) have been established, we obtain by standard properties 
of C° semigroups the bounds [Se^S -1 ] < e 70 * and \e Lt \ < Ce 70 * and also 
the inverse Laplace transform formulae 

(3.22) e Lt f(x) = P.V. / e x \\-L)~ l fd\ 

J 7— ioo 

for / G V(L) and 

(3.23) / e L ( T ~ t )/(x) = P.V. / e A *(A-L)" 1 / T (A)dA 

JO J-y—ioo 

for / G L 1 ([0, T]; P(L)), for any 7 > 70, where / denotes Laplace transform 
of /; see Propositions lA.lll IA.24| and I A. 2 5l Appendix A. Formulae (|3.12|) - 
(|3.13j) then follow from ()3.22|) - l|3.23|) by Fourier transform/distribution the- 
ory. Thus, it remains only to verify (|3.20|) in order to complete the proof. 
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To establish (|3.20j) . rewrite resolvent equation (A — L)U = f in W- 
coordinates, as 

(3.24) XA°W + &W Xj + Y,(B jk W Xk ) x . +CW + Y, = A°f, 

3 jk j 

where / := (dW/dU)(U(xi))f, C is uniformly bounded, and coefficients 
A 3 , B jk satisfy (Al) with (5 = 0. (Exercise: check that properties (Al) 
are preserved up to lower-order terms CW + ^ w x-i using commutation 
of coordinate-change and linearization.) Using block-diagonal form of ^4°, 
we arrange by coordinate change W — > (A^^^W if necessary that A = I 
without loss of generality; note that lower-order commutator terms arising 
through change of coordinates are again of the form CW + D^w^. 1 ., hence 
do not change the structure asserted in (|3.24|) . 

Taking the real part of the complex inner product of W against (|3.24|) . 
and using (|1.28[) . (|1.30|) together with symmetry of A\ x terms and Young's 
inequality, similarly as in the proof of Proposition II . 16] we obtain 

(3.25) ^X\W\ L 2 + 6\w n \ 2 m <7o|W|| 2 + |/| L2 |Ty| L2 

for 70 > sufficiently large. Dropping the favorable w 11 term and dividing 
both sides by |W|/,2, we obtain 

(3.26) (m - l0 )\W\ L 2 < \f\ L 2 

Noting that S(X — L)S~ 1 W = /, we are done. □ 

Remark 3.8. The properties and techniques used to establish Proposition 
13.61 are linearized versions of the ones used to establish Proposition 11.161 
For a still more direct link, see Remark IA. 181 and Exercise IA.221 Appendix 
A, concerning the Lumer-Phillips Theorem, which assert that, in the favor- 
able VF-coordinates, existence of a C° contraction semigroup \e Lt \ < e 70 * is 
equivalent to the a priori estimate 

(3.27) (d/dt)(l/2)\W\ 2 L2 = K(u,Lu) < 7o|^| 2 , 

the autonomous version of ()1.17j) . That is, Proposition 13.61 essentially con- 
cerns local linearized well-posedness, with no global stability properties ei- 
ther asserted or required so far. 

Besides validating the spectral resolution formulae, the semigroup frame- 
work gives also a convenient means to generate solutions of the nonau- 
tonomous linear equations arising in the nonlinear iteration schemes de- 
scribed in Section 1, based only on the already- verified instantaneous ver- 
sion of (|3.27|) : specifically, Proposition IA.301 Appendix A, assuming (|3.27|) . 
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guarantees a solution satisfying |n(t)|^2 < e 7o *|u(0)|^2. A standard way of 
constructing the autonomous semigroup is by semidiscrete approximation 
using the first-order implicit Euler scheme (see Remark IA.12|) . the nonau- 
tonomous solution operator then being approximated by a concatenation of 
frozen-coefficient semigroup solutions on each mesh block; see [Pa]. It is 
interesting to compare this approach to Friedrichs' original construction of 
solutions by finite difference approximation [Fr]. 

3.3. Asymptotic ODE theory: the gap and conjugation lemmas. 

Consider a general family of first-order ODE 

(3.28) W'-A(x,A)W = F 

indexed by a spectral parameter A G Q C C m , where W € C^, x S E and 
'"" denotes d/dx. 

Examples 3.9. 1. Eigenvalue equation (L^ — \)U = 0, written in phase 

coordinates W := (W,{w n )'), W := {dW/dU){U(x))U , with A := (£, A) 
and F := 0. 

2. Resolvent equation (L? — X)U = f, written in phase coordinates 
W := (W, (w 11 )'), W := {dW/dU){U(x))U, and |F| = 0(\f\). 

Assumption 3.10. 

(hO) Coefficient A(-, A), considered as a function from f2 into C°(x) is 
analytic in A. Moreover, A(-,A) approaches exponentially to limits A± as 
x — > ±oo ; with uniform exponential decay estimates 

(3.29) \(d/dx) k (A - A±)| < C ie ~ 9 ^ /C2 , for x ^ 0, < k < K, 
Cj, 6 > 0, on compact subsets o/O. 

Lemma 3.11 (The gap lemma [KS, GZ, ZH]). Consider the homo- 
geneous version F = of (|3.28|) . under assumption (hO). IfV~(A) is an 
eigenvector of A- with eigenvalue /j,(A), both analytic in A, then there exists 
a solution of (|3.28|) of form 

(3.30) W(A, x) = V(x, A)e^ A)x , 

where V is C 1 in x and locally analytic in A and, for any fixed 6 < 9, 
satisfies 

(3.31) V{x,A) = V-(A) + 0(e- § W\V~(A)\), x<0. 
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Proof. Setting W(x) = e' lx V(x), we may rewrite W' = AW as 

(3.32) V' = (A_ - fil) V + 6V, 9 := (A- A_) = 0(e~ e \ x \), 

and seek a solution V(x,A) — > V - (a;) as x — > oo. Choose 9 < 9\ < 6 such 
that there is a spectral gap |K(<tA_ — (/i + #i))| > between <rA_ and 
+ 9\. Then, fixing a base point Ao, we can define on some neighborhood 
of Ao to the complementary A_-invariant projections -P(A) and Q(A) where 
P projects onto the direct sum of all eigenspaces of A_ with eigenvalues 
p, satisfying 9i(/2) < 5ft(/u) + 9\, and Q projects onto the direct sum of the 
remaining eigenspaces, with eigenvalues satisfying 5?(/i) > + 6\. By 

basic matrix perturbation theory (eg. [Kat]) it follows that P and Q are 
analytic in a neighborhood of Ao, with 
(3.33) 

e (A—ni)x p < C {e 01x ), x>0, e iA —» I)x Q <C{e 8lx ), x<0. 
It follows that, for M > sufficiently large, the map T defined by 
TV(x) = V~ + f X e i ~ A —^ x -^P9{y)V{y)dy 



(3.34) 



-oo 
-M 

e^-^-v)Q0(y)V(y)dy 



is a contraction on L°°(— oo, —M]. For, applying (|3.33j) . we have 
(3.35) 

\TVi -TV 2 \ {x) <C\V 1 -V 2 \oo(J X e^-y^e^dy + J e^^e^dy) 

< c 2 \Vi - v 2 \ooe~° M < hy x - y 2 u. 

By iteration, we thus obtain a solution V £ L°°(—oo, —M] of V = TV 
with V < C3IV - " |; since T clearly preserves analyticity V(A,x) is analytic 
in A as the uniform limit of analytic iterates (starting with Vb = 0). Differ- 
entiation shows that V is a bounded solution of V = TV if and only if it is 
a bounded solution of (|3.32j) (exercise). Further, taking V% = V, V 2 = in 
(|3.35j) . we obtain from the second to last inequality that 

(3.36) \V-V~\ = \T(V)-T(0)\ < C 2 e 5x \V\ < cJ x \V~\, 

giving (|3.31|) . Analyticity, and the bounds 1)3.31(1 . extend to x < by 
standard analytic dependence for the initial value problem at x = — M. □ 
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Remark 3.12. The title "gap lemma" alludes to the fact that we do not 
make the usual assumption of a spectral gap between //(A) and the remaining 
eigenvalues of A_, as in standard results on asymptotic behavior of ODE 
[Co]; that is, the lemma asserts that exponential decay of A can substitute 
for a spectral gap. Note also that we require only analyticity of fi and not 
its associated eigenprojection ILj, allowing crossing eigenvalues of arbitrary 
type (recall, 11^ is analytic only if \x is semisimple; indeed, 11^ blows up 
at a nontrivial Jordan block [Kat]). This is important in the following 
application; see Exercise 13, 161 below. 

Corollary 3.13 (The conjugation lemma [MeZ.l]). Given (hO), there 
exist locally to any given Ao € £1 invertible linear transformations P+[x, A) = 
I + 0+(x, A) and P—(x, A) = I + G_(x, A) defined on x > and x < 0, 
respectively, <f>± analytic in A as functions from Q, to C°[0,±oo), such that: 

(i) For any fixed < 9 < 9 and < k < K + 1, j > 0, 

(3.37) \(d/dA) j (d/dx) k Q± \ < C(j)CiC 2 e- e|l|/C2 for x ^ 0. 

(ii) The change of coordinates W =: P±L, F =: P±G reduces (|3.28f) to 

(3.38) .Z' - A ± Z = G for x ^ 0. 
Equivalently, solutions of (|3.28f) may be factored as 

(3.39) w=(i+e±)z ± , 

where Z± satisfy the limiting, constant-coefficient equations (|3.38j) and Q± 
satisfy bounds (|3.37|) . 

Proof Substituting W = P-Z into equating to JHIHHJ, and rearrang- 

ing, we obtain the defining equation 

(3.40) PL = A-P- - P-A, P -» I as x -> -oo. 

Viewed as a vector equation, this has the form P'_ = AP-, where A ap- 
proaches exponentially as x — > — oo to its limit A-, defined by 

(3.41) A-P := A_P-PA_. 

The limiting operator A- evidently has analytic eigenvalue, eigenvector pair 
[i = 0, P- = I, whence the result follows by Lemma l3.11l for j = k = 0. The 
x-derivative bounds < k < K + 1 then follow from the ODE and its first 
K derivatives, and the A-derivative bounds from standard interior estimates 
for analytic functions. A symmetric argument gives the result for P + . □ 
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Remark 3.14. Equation (|3,39|) gives an explicit connection to the invis- 
cid, bi-constant-coefficient case, for bounded frequencies |(£, A)| < R (low 
frequency ~ inviscid regime). For high-frequencies, the proper analogy is 
rather to the frozen-coefficient case of local existence theory. 

Exercise 3.15 ([B]). Use Duhamel's formula to show that 

(3.42) PV{x) = V- + ! e (A— /*/)(*-!/) pe(y)V(y)dy 

J — oo 

and 

(3.43) QV(x) = QV(-M) + f e^-^ x -^Q6(y)V(y)dy, 

J-M 

hence V(x) = TV(x) for the unique solution V of (|3.32j) determined by 
conditions V(-oo) = V_ = PV- and QV(-M) = 0. 

Exercise 3.16. (i) If {rj} and {Ik} are dual bases of (possibly generalized) 
right and left eigenvectors of A_, with associated eigenvalues fj,j and 
show that {rjl* k } is a basis of (possibly generalized) right eigenvectors of the 
operator A- defined in (|3.41[) . with associated eigenvalues fij — [i k . (ii) Show 
that n = is a semisimple eigenvalue of A~ if and only if each eigenvalue 
o/A_ is semisimple. 

3.4. Hyperbolic parabolic smoothing. We next introduce the circle 
of ideas associated with hyperbolic-parabolic smoothing and estimate ()1.40(l . 

Lemma 3.17 ([Hu, MaZ.5]). Let A = block-diag {ajl mj } be diagonal, 
with real entries dj appearing with prescribed multiplicities rrij in order of 
increasing size, and let B be arbitrary. Then, there exists a smooth skew- 
symmetric matrix-valued function K(A, B) such that 

(3.44) Re (B - KA) = 5?block-diag B, 

where block-diag B denotes the block- diagonal part of B, with blocks of di- 
mension mj equal to the multiplicity of the corresponding eigenvalues of A. 

Proof. It is straightforward to check that the symmetric matrix $IKA = 
(1/2) (KA — AtR) may be prescribed arbitrarily on off-diagonal blocks, by 
setting Kij := [ai — aj) _1 Mjj, where My is the desired block, i ^ j. Choos- 
ing M = OlB, we obtain R(B - KA) = !ftblock-diag (B) as claimed. □ 

Lemma 3.18 ([KSh]). Let A , A, and B denote real-valued matrices such 
that A is symmetric positive definite and A := A A and B := A°B are 
symmetric, B > 0. Then, the following are equivalent: 



37 



(KO) (Genuine coupling) No eigenvector of A lies inker B (equivalently, 
in ker B ). 

(Kl) block-diag LBR = R*BR > 0, where L := 0*(I°) 1/2 and R := 
are matrices of left and right eigenvectors of A block- diagonalizing 
LAR, with O orthonormal. Here, as in Lemma \3.11\ block-diag M denotes 
the matrix formed from the diagonal blocks of M, with blocks of dimension 
equal to the multiplicity of corresponding eigenvalues of LAR. 

(K2) (hyperbolic compensation) There exists a smooth skew- symmetric 
matrix-valued function K {A, B, A ) such that 

(3.45) Re [B - KA^j > 0. 

(K3) (Strict dissipativity) 

(3.46) ^a{-iiA-\i\ 2 B) < -#j£| 2 /(l + |£| 2 ), 6 > 0. 

Proof. (KO) 44> (Kl) by the property of symmetric nonnegative matrices M 
that v l Mv = if and only if Mv = (exercise), so that block-diag {R^MBRj} 
has a kernel if and only if ^Rj^RBRja = 0, if and only if WtBRja = 0, where 
Rj denotes a block of eigenvectors with common eigenvalue. 

(Kl) => (K2) follows readily from Lemma l3.181 by first converting to the 
case of symmetric A and B < by the transformations A — > (A ) 1 / 2 AA )^ 1 / 2 , 
B — > (A ) 1 / 2 BA9)~ 1 I 2 , from which the original result follows by the fact that 
M > 44> (A°) 1 ^ 2 M(A ) 1 ^ 2 > 0, then converting by an orthonormal change 
of coordinates to the case that A is diagonal and B < 0. Variable multiplic- 
ity eigenvalues may be handled by partition of unity/interpolation, noting 
that $t(B — KA) < persists under perturbation. 

(K2) => (K3) follows upon rearrangement of energy estimate 

= 5R<((C + |£| 2 )i° + i£K)w, (A + i£A + \£\ 2 B)w) 

(3.47) = ft\(w, ((C + |£| 2 )i° + i£K\w) + \£\ 2 {w, »(B - KA)w) 

+ U(w, -i\^K{A Q )- x Bw) + C\£\ 4 (w,Bw), 

which yields 

$t\{w,((l + C\Z\ 2 )A° + i£K^jw) 

( 3 - 48 ) < -(o\£\ 2 \w\ 2 + M\£\ 3 (\w\\Bw\ + (C/M) \£\ 4 \Bw\ 2 } 

< -6\£\ 2 \w\ 2 
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and thereby 

(3.49) »A(1 + \Z\ 2 )M 2 < -e^lw] 2 , 

for M, C > sufficiently large and 9, 9\ > sufficiently small, by positivity 
of i° and ^(S-iOl). 

Finally, (K3) =r- (Kl) follows by by first-order Taylor expansion at £ = 
of the spectrum of LAR — i^LBR (well-defined, by symmetry of LAR), 
together with symmetry of R t BR. □ 

Corollary 3.19. Under (Al)-(A2), there holds the uniform dissipativity 
condition 

(3.50) 3ME^ " J2^B jk )± ^ ~^| 2 /(1 + l^l 2 )- 

Moreover, there exist smooth skew- symmetric "compensating matrices" K±(^), 
homogeneous degree one in such that 

(3.51) ft^fc&B'* ~ KiOiA )- 1 ^^)^ >0>0 

j,k k 

for all£e R d \ {0}. 

Proof. By the block-diagonal structure of B jk (GC) holds also for A± and 

Qjk ._ (A°)- l $lBi k , since 

(3.52) 

ker J2 C]Ck& k = ker £ = ker ]T = ker £ 

Applying Lemma 13.181 to 
(3.53) 

A°:=Al, A:=((A^^A k ), B : ((.*"; 1 ^ 



we thus obtain 1)3. 5 lj) and 



(3.54) ^[(I )-^ »> V ' ~ E^^ ifc )] ± ^ - il^| 2 /(l + Iff 
01 > 0, from which we readily obtain 

(3.55) ( y,'^ aj - 2>&# fc )± ^ -« 2 /u + iei 2 ) 

3 j,k 
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and thus l)3.5U|) (exercise, using M > 9i (i°) ± 1/2 M(i°) ± 1/2 > and 
a(A°)± 1/2 M(A°)± 1/2 > 6 o a^^M > 0, together with S > 9 & aS > 6 
for S symmetric). Because all terms other than K in the lefthand side of 
1)3.51)1 are homogeneous, it is evident that we may choose K (•) homogeneous 
as well (restrict to the unit sphere, then take homogeneous extension). □ 

3.4.1. Basic estimate. Energy estimate (|3.47|) - (|3.49)l may be recog- 
nized as the Laplace-Fourier transformed version of a corresponding time- 
evolutionary estimate 
(3.56) 

(d/dt)(l/2)((A°W,W) + (A°W X ,W X ) + (Kd x W,W)) 

= -(W X MB-KA)W X ) - (W^KiA^BW^-CiW^BW^), 
< -e(\W x \ 2 L2 + \BW x \ 2 m ) 

for the one-dimensional, linear constant-coefficient equation 

(3.57) A°W t + AW x = BW xx , 

which in the block-diagonal case B = block-diag {0,6}, b > 0, yields 1)1.40)1 
for s = 1 by the observation that 

(3.58) £{W) := (1/2) ({A°W, W) + {Kd x W,W)\ + C(A°W X ,W X ) 

for C > sufficiently large determines a norm £ 1//2 (-) equivalent to | • 
This readily generalizes to 

(3.59) (d/dt)S(t) < -9{\d x W\ 2 H .-i + \Bd x W\ 2 Ha ), 



£{W) := (l/2)( (A°W,W) + (Kd x W,W)} 

, , +c((A°W x ,W x ) + (Kd x W x ,W x ))+--- + 

(3.60) V ') 

+ c a - 1 (iA Q d?- 1 w,£r 1 w) + (Kd^-^d^w)) 

+ C s {A d s x W,d s x W), 

where £ 1 ^ 2 (•) is a norm equivalent to \ -\h s i yielding 1)1. 40|) for s > 1. We refer 
to energy estimates of the general type (K3), ()3.56j) - ()3.60)) as "Kawashima- 
type" estimates. 
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Proposition 3.20. Assuming (A1)-(A2), (HO) for W- = 0, energy esti- 
mate Q1.4UJI is valid for all q(d) < s < q, q(d) and q as defined in (HO), so 
long as \W\h s remains sufficiently small. 

Proof. In the linear, constant-coefficient case, Q3.51JI together with a calcu- 
lation analogous to that of (|3.56j) and (|3.59j) yields 

(3.61) (d/dt)£(t) < -0{\d x W\ 2 H „^ + \d x w n \ 2 H s) 

for 

£{W) := (1/2) ((1° W,W) + (K-{d x )W,W) 

+ c((A°_d x w, d x w) + (K. (d x )d x w, d x w) )+••• + 



(3.62) 

+ C s - 1 [{A a _d s ~ 1 W,d x - 1 W) + (K_{d x )d'- x W,&'- x W) 



+ c s {A°_d s x w,d s x w), 

where operator K_(d x ) is defined by 

(3.63) K^f){0--=iK-(0f(0 

with K_ (■) as in (|3.51[) . where g denotes Fourier transform of g, for all s > 1. 

In the general (nonlinear, variable-coefficient) case, we obtain by a sim- 
ilar calculation 

(3.64) (d/dt)£(t) < -6{\d x W\ 2 Hs -, + \d x w u \ 2 H s) + C\W\\ 2 

for C > sufficiently large, for 
(3.65) 

S(W) := (1/2) ((A°(W)W,W) + (K_(d x )W,W) 



+ C( K (A°(W)d x W, 8 X W) + (K-(d x )d x W, 8 x W)j +■■■ + 

+ C s - 1 ((A°(W)d s x - 1 W,d s x - 1 W) + (K-{d x )%T x W,dZr x W) 
+ C s (A°(W)d x W,d s x W), 

provided s > q(d) := [d/2] + 2 and the coefficient functions possess sufficient 
regularity C s , i.e., s < q. 

Here, we estimate time-derivatives of A terms in straightforward fash- 
ion, using integration by parts as in the proof of Proposition 11.61 We esti- 
mate terms involving the nondifferential operator K~(d x ) in the frequency 
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domain as 
(3.66) 

(d/dt)(l/2){K(d x )dlW,ff x W) = (d/dt)(l/2){iK(0(iO r W,m r W) 

= {iK(0(iO r W,(iO r W t ) 

= ((^ r w,-^(0(i°_r 1 (E^ i -)(^) r ^) 

3 

using Plancherel's identity together with the equation, written in the fre- 
quency domain as 

(3.67) W t = -J2 %M°_ Y l A j _W + H, 

j 

where 



(3.68) 



j 

+ Y J {^)-\& k W Xk ) X] +(A°)- 1 G{W X) W x \ 



By a calculation similar to those in the proofs of Propositions ll.6| and ll.lb1 
(exercise, using smallness of \A°_ — A 3 (W) \ ~ \W\ and the Moser inequality 
Q1.15|l ). we obtain 

(3.69) \d r x H\ L 2 < C\d r x +2 w n \ L 2 +C\W\ Hr+ i(\W\ Hld/2]+ 2 + \W\ r Hld/2]+2 ). 

Thus, using homogeneity, |-ftT(£)| < C|£|, together with the Cauchy-Schwartz 
inequality and Plancherel's identity, we may estimate the final term on the 
righthand side of ()3.66|) as 

((iKmYWMYH) < C\dl +1 W\ L2 \\d r x H\ L i 

< C\d r x +1 W\ L2 \d r x +2 w n \ L 2 +e\W\ 2 Hr+1 , 

any e > 0, for |VF|^[,i/2]+2 sufficiently small: that is, a term of the same form 
arising in the constant-coefficient case plus an absorbable error. 

Combining A - and X(^)_-term estimates, and using (|3.51|) together 
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with 
(3.71) 

J2(d r x W Xj ,& k (WW x W Xk ) = (m r W,J2tek& k (iO r W) 

3,k j,k 

+ o(\d r x +1 w\ 2 L2 \w\ Lac ) 

> ((iO r W,J2^kBL\iO r W)-e\d x +1 W\ 2 L2 

for any e > 0, for |V^|loo sufficiently small, we obtain the result, similarly 
as in the constant-coefficient case. □ 

Exercise 3.21. (Alternative proof) Show using (fl~34l together with (l3~66l) - 
(I3~7B that 

(3.72) (d/dt)£ < C\W\ 2 Hs ., - OiQdZWfa + \d x +1 w n \ L 2) 
for C > sufficiently large and 9± > sufficiently small, for 

s 

(3.73) £{W) := (1/2) + (l/2C)(i$T_ (^J^W, a*" 1 ^)- 

r=0 

Using the H s interpolation formula 
(3.74) 

|/|^* ^/^C^I/l^+^-^C- 1 /^)!/!^, / 9=( S2 - S ,)/(«2-«i), 

valid for s\ < s* < S2, taking s\ = 0, S2 = s, and C sufficiently large, show 
that lj3~72|l zmp/ies (j31>5j> . 

3.4.2. Alternative formulation. A basic but apparently new obser- 
vation is that energy estimate (|3.65|) implies not only (|1.4Uj) but also the 
following considerably stronger estimate. 

Proposition 3.22. Under the assumptions of Provosition Iff 2 (A 
(3.75) 

\W(t)\ 2 H s < C{e- e ^\W(0)\ 2 Hs + f e^ {t - s) \W{s)\ 2 L2 ds) , 9 1 > 0, 

J o 

so long as \W\h 3 remains sufficiently small. 

Proof. Rewriting (|3.64|) using equivalence of £ 1 ^ 2 and H s as 

(3.76) (d/dt)£(t) < -OtSit) + C\W(t)\ 2 L2 , 

we obtain by Gronwall's inequality e eit £\ < C J* * e eiS |W(s)|£ 2 ds and thereby 
the result (exercise). □ 
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Proof of Proposition \1.22L Exercise, using (|3.75|) together with the fact that 
\U{t)\ L 2 < C\U(0)\ L 2 for \W\lo* sufficiently small (Exercise EH)) . □ 

In Appendix B, we sketch also a simple proof of Proposition 11,211 using 
(|3.75j) together with linearized stability estimates. This may be helpful for 
the reader in motivating the nonlinear stability argument carried out in 
Section 4 for the more complicated variable-coefficient situation of a viscous 
shock profile. 
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4 Reduction to low frequency 



We now begin our main analysis, carrying out high- and mid- frequency esti- 
mates reducing the nonlinear long-time stability problem (LT) to the study 
of the linear resolvent equation (A — L^)U = f in the low- frequency regime 

|(£, A) | small described in Proposition 12.141 The novelty of the present treat- 
ment lies in the simplified argument structure based entirely on energy es- 
timates. More detailed linearized estimates were obtained in [Z.4] using 
energy estimates together with the pointwise machinery of [MaZ.3]. 

4.1. Nonlinear estimate. Define the nonlinear perturbation 

(4.1) U:=U-U, 

where U denotes a solution of (|1.2|) with initial data Uq close to U. Our first 
step is to establish the following large- variation version of Proposition 13.221 

Proposition 4.1 ([MaZ.4, Z.4]). Given (A1)-(A2), (HO)-(HS), and As- 
sumvtion \2.8l 

(4.2) \U(t)\ 2 HS <C{e- e2t \U(0)\ 2 HS + f e- 9 ^\U(s)\l 2 ds), 2 > 0, 

J o 

so long as \U\h" remains sufficiently small, for all q(d) < s < q, q(d) and q 
as defined in (HO). 

4.1.1. "Goodman-type" weighted energy estimate. The obvi- 
ous difficulty in proving Proposition 14.11 is that (A2) is assumed only at 
endpoints W±, hence we cannot hope to establish smoothing in the key 
hyperbolic modes w 1 through the circle of ideas discussed in Section 3.4, 
except in the far field x\ — > ±oo. The complementary idea needed to treat 
the "near field" or "interior" region x\ £ [— M, M] is that propagation in the 
hyperbolic modes thanks to assumption (HI) is uniformly transverse to the 
shock profile, so that signals essentially spend only finite time in the near 
field, at all other times experiencing smoothing properties of the far field. 

The above observation may be conveniently quantified using a type of 
weighted-norm energy estimate introduced by Goodman [Go. 1-2] in the 
study of one-dimensional stability of small- amplitude shock waves. Con- 
sider a linear hyperbolic equation A°Wt + Ylj=i A*W Xj = with large- but 
finite- variation coefficients depending only on x±, 

(4.3) |4J<©, [&(y)dy<oo, j = 0,l,...d, 
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A? symmetric, ^4° > 0, under the "upwind" assumption 
(4.4) A 1 > 9 > 0. 

Define the scalar weight a{x\) := e^o 1 -( 2C ®/ e )(y) d v ^ positive and bounded 
above and below by (|4.3j) . where C > is sufficiently large. 

Then, the zero-order Goodman's estimate, in this simple setting, is just 

(d/dt)(l/2)(W,aA°W) = {W,aA°W t ) 

= -(W,a^A j W Xj ) 

3 

( ' 5) =(l/2){W,(aA 1 ) Xl W) 

= (l/2){W,a(-(2CQ/e)A 1 + A l Xi )W) 
< -(C/2)(W,aOW), 

yielding time-exponential decay in | W\^2 on any finite interval x\ € [— M, M] . 
Exercise 4.2. Verify the corresponding s-order estimate 

s s 

(4.6) (d/dt)(l/2)Y,(d r x W,aA°d r x W) < -(C/2) £>£W, aQd^W) 

r=0 r=0 

for C > sufficiently large, yielding time- exponential decay in \W\h s on 
any finite interval x\ G [— M, M]. 

Remark 4.3. Proper accounting of the favorable effects of transverse prop- 
agation is a recurring theme in the analysis of stability of viscous waves; see, 
e.g., [Go.1-2, L.4, LZ.1-2, LZe.l, SzZ, L.l, ZH, BiB, GrR, MeZ.l, MGWZ.b 
4]- 

4.1.2. Large-amplitude hyperbolic parabolic smoothing. Com- 
bining the weighted energy estimate just described above with the techniques 
of Section 3.4, we are now ready to establish Proposition 14.11 

Proof of Proposition \4-l\ Equivalently, we establish 

(4.7) \W(t)\ 2 H s < C(e~ e ^\W(0)\ 2 HS + f ' e- e ^-^\W(s)\l 2 ds), 9 1 > 0, 

J o 

where 

(4.8) W(x, t) := W(x, t) - W{xi) := W(U{x, t)) - W{U{ Xl )). 
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(Exercise: Show |£/|#r- ~ |W|ip- for all < r < s, in particular and s, 
using s > [d/2] +2 and same estimates used to close the energy estimate in 
the proof of Proposition 11.161 ) 

Noting that both W and W satisfy (|1.27j) . we obtain by straightforward 
calculation [MaZ.4, Z.4] the nonlinear perturbation equation 

(4.9) 3 j ' k 

+ M l W X1 + ( M 2 W X1 ) Xj + M l W *i + G ( x > d - W ) > 

3 3 

where 

(4.10) Ai(x,t) :=Ai(W(x,t)), B jk {x,t) := B jk (W(x,t)), 



(4.11) M 1 =M 1 (W,W):=A 1 (W)-A 1 (W) = (J dA\W + 9W)dB\w 



(4.12) Ml = M 3 JW, W) := B jl - B j = i , rl n , , . 



(4.13) M{ = Mi( Xl ) := {dG/dW X3 )\ dxW = ( °, °, 



and 



(4.14) 



G 



9 J ' \g(d x w, d x w) - g(d x W, d x W) - dg(W)d x W 


o(\d x w\ 2 



(In the case that there exists a global convex viscosity-compatible entropy, 
we may further arrange that G = 0.) 

To clarify the argument, we drop terms M and G in (|4.9[) . The reader 
may verify that these generate harmless, absorbable error terms in our en- 
ergy estimates (exercise; see [MaZ.4, Z.4]). With this simplification, we 
reduce to the familiar situation 

(4.15) A°W t + Y #W Xj = Y,(& k W Xk ) Xj , 

3 j,k 
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with the difference that A* = A> (W + W) and B jk = & k (W + W) are no 
longer approximately constant, but for small |W|#s approach the possibly 
rapidly varying values A*(W) and W (W) along the background profile. 
By Corollary 13.31 and Assumption (HI), respectively, we have 

(4.16) \(d/d Xl ) k AHW)\, \{d/d Xl ) k & k (W)\ < 8(xi) := Ce- 9 ^ 
and 

(4.17) A 1 > 6 > 0, 
for some C, 6 > 0. Define 

(4.18) a{xi) := e-> 

where C* > is a sufficiently large constant to be determined later. 

Decompose W into the sum of near- and far-field parts Wn and Wp, 

w = w N + wp~ + Wp 

:= (l - X + (xi) ~ X~(xi))W + X + (xi)W + X ~(xi)W, 

where x^ 1 are C°° cutoff functions supported on x\ > M and x\ < —M and 
one on x\ > 2M and x\ < — 2M, respectively, for M > sufficiently large, 
with 

(4.20) Kx^i)! <emin{| X ± |,|l-X ± |}, < r < s + 1, 

where e > is a sufficiently small constant to be determined later. Up to 
lower-derivative error terms that are harmless for e sufficiently small, both 
Wn and Wp satisfy the same equation (|4.15j) satisfied by W. Again, we'll 
drop these terms, leaving it to the reader to verify that they can indeed 
be absorbed in the estimates. (The main point is to notice that first-order 
contributions dxB lk W Xk have zero first row, hence do not affect the 
structural assumption A^ symmetric.) 

On the respective supports x\ > M, x\ < —M of Wp, coefficients A> 
and & k remain arbitrarily close in H s to A J ± and B± provided |W|#s is 
taken sufficiently small and M sufficiently large. Thus, using the argument 
of Proposition 13.201 we obtain 

(4.21) (d/dt)£±(W±(t)) < -9{\d x W±\ 2 Hs ^ + IdxWp 1 ^*) + C\W±\ 2 L2 
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for C > sufficiently large, for 

£$(W) := {1/2)({A°W,W) + (K±(d x )W,W)) +C{A°d x W,d x W) 

( 4 - 22 ) + C(K±{d x )d x W,d x W)) + C 2 (A°d 2 W,d 2 W) + ... 

+ C s (A°d s x W,d s x W), 

where operator K±(d x ) is defined by 

(4.23) K^m):=iK±(Of(0 

with K±(-) as in Q3.51JI . where g denotes Fourier transform of g. 

To estimate Wn, we use a different technique, substituting Goodman- for 
Kawashima-type estimates in controlling wj^ terms. Combining calculations 
(fOS)) and (JJU), we obtain 
(4.24) 

^(W N ,aA°W N ) t = ^(W N ,aA° t W N ) + (W N ,aA°(W N ) t ) 
= ±{W N ,aA° t W N ) - (W N ,Y,<xA j (W N ) Xj + J2 a (B jk (W N ) Xk ) Xj ) 

3 j,k 

= {^(w T N , (aA l u (W)) X1 - ^<«% , a& k (wj^) Xk } } 

+ l(W N ,aJi$W N ) + K,^ (a^ - Ai^W)))^) 

j 

j i 
+ i(<,^aIi 2 (<), j }-^(<,a, lJ B lfc (<) Xfe )} 

< aeu;#) +9{d x w I J,ad x wjj)^ 

+ C(|Wi| 0O + |W x | o0 ) y a|W| 2 + C y al^w^HWl 

+ cc,y eaK J |(|^| + |^< 7 l)> 

for some c > 0, C > independent of C*. Choosing C* sufficiently large 
that C* >> (7(1 + and C* >> C 2 /6Q on the support [-2M, 2M] of 

Wat, and estimating the unbracketed terms in the final line of (|4,24|) using 
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Young's inequality, we find that 

1 ~ C 9 

(4 25) 2^'°^°^* - ~i[( wI N' aQvjI N) ~ ^(dxWN,ad x w^) 

+ C(C*)(w I r j,aw I J). 

More generally, defining 

s 

(4.26) £ N (W) := (1/2) £ c~ r (d r x W, aA^W), 

r=0 

for c = c(C*) > AC*/ 9 > sufficiently large, we obtain by telescoping sum 

that 

(4.27) 

{d/dt)£ N {W N (t)) <J2 c ' r (- Y^ d >N^ ed >N)) 

r=0 

+ E c ~ r ( " li^N^dl+Wjj) + C(a)(5>if,a^O 

r=0 



01 > and Ci > depending on C*, where in the second inequality of Q4.27B 
we have used the fact that a and are bounded above and below on the 
support [-2M,2M] of W N . 

Summing Q4.21JI and (|4.27|) . and using the fact that (£f) 1/2 and (£ N ) 1/2 
are each equivalent to H s , we obtain 

(4.28) (d/dt)£(W(t)) < -0 2 £(W(t)) + C 2 \W(t)\ 2 L2 , 
where 

(4.29) £{W) := £p{Wp) + £ N (W N ) + £p(Wp) 

with £ l / 2 (W) equivalent to |W|#s. Applying Gronwall's inequality, similarly 
as in the proof of Proposition 13^21 we obtain e e *£\% < C f* e 02S \W {s)\\ 2 ds 
and thereby the result. □ 

Remark 4.4. Note that Q4.2|) gives pointwise- in-time control on the H s (x) 
norm of W in terms of a time- weighted L 2 (t) average of the L 2 (x) norm. We 
will take advantage of this fact in the linearized analysis of Section 4.3, where 
we establish L 2 -time-averaged high-frequency estimates by a simple Parseval 
argument, thus avoiding the complicated pointwise bounds of [MaZ.4, Z.4]. 
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4.2. Linearized estimate. We next establish the following estimate 
on the linearized inhomogeneous problem 

(4.30) U t -LU = h + d x f 2 , U(0) = U , 

where L is defined as in (|3.8[) . assuming the low- frequency bounds of [Z.3]. 

Proposition 4.5. Given (A1)-(A2), (H0)-(H5), assumption (|2.8|). and 

strong spectral, structural, and strong refined dynamical stability, the so- 
lution U{t) := Jle L{t - s ^f{s)ds of satisfies 
(4.31) 

[ T e- 8e ^ T -^\U(s)\ 2 L2 ds < Cil + Ty^l^U^ 
Jo 

+ Ce- 2 ^ T {\U \ 2 Hl + \LU \ 2 m ) 

+ C f\~ 2 ^\\fM\ 2 m + \dM 2 m)ds 

+ c( j\l + T- s)- (d - 1)/4+e |/i( S )| L i ds)" 

+ c[f\i +t- s) -cd-D/4+ 6 -i/ 2 | /2(s)Ui ds y 

/or any fixed e > 0, /or some C = C(e) > sufficiently large, for any Uq, 
fli f2 such that the righthand side is well-defined. In the case of strong 
(inviscid, not refined) dynamical stability, we may take e = and C > 
fixed. 

Remark 4.6. Somewhat stronger linearized bounds were obtained in [Z.4] 
by a more detailed analysis (pointwise in time, with no loss in regularity, 
and for the initial-data rather than the zero-data inhomogeneous problem). 
However, these will suffice for our argument and are much easier to obtain. 

4.2.1. High- and mid-frequency resolvent bounds. Our first step 
is to estimate solutions of resolvent equation (J3.2U|) . 

Proposition 4.7. (High-frequency bound) Given (A1)-(A2) and (HO)-(Hl), 

(4.32) KA-^rV^) ^ C for |(f,A)| >fland»A>-0, 

for some R, C > sufficiently large and 9 > sufficiently small, where 
| • |^i denotes H 1 operator norm, \f\fji ■= |(1 + \d xi \ + |^|)/|l2. 
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Proof. A Laplace-Fourier transformed version (now with respect to (t, x)) 
of nonlinear estimate (|4.27|) . s = 1, carried out on the linearized equations 
(l3~%l) written in IV-coordinates, W := (dW/dU)(U)U, yields 
(4.33) 

ka((i + \l?)\w\ 2 +\d xi w\ 2 ) < -0i(i£i 2 |wf + l^^i 2 ) 

+ c x (\w\ 2 + (i + \£\ 2 )\w\\f\ + \d xi w\\d x j\) 

for some C\ > sufficiently big and 9\ > sufficiently small, where | • | 
denotes | • \l 2 ( Xi ) an d 

(4.34) (A - L~)U = f. 

We omit the proof, which is just the translation into frequency domain of the 
proof of Proposition 14. 1( carried out in the much simpler, linearized setting. 
For related calculations, see the proof of Proposition l3~o1 or the proof of (K3) 
in Lemma 13.181 

Rearranging, and dividing by factor |(1 + |£| + |9 :ri )|P^|, we have 

(4.35) + + ||| + \d xl \)W\ < Ci|((l + \i\ + \d xi \)f\ + d\W\, 

which yields (IP2l with 9 := 9 1 /2 (i.e, 9?A + 6>i > 9 x /2 > 0) and C := 2Ci/0, 
for |(£,A)| > 2Ci (i.e., sufficiently large that C\W may be absorbed in the 
lefthand side of flQ5J) ). □ 

Proposition 4.8. (Mid-frequency bound) Given (A1)-(A2), (HO)-(Hl), 
and strong spectral stability (|2.14|) . 

(4.36) \{\- L^)- 1 \ fll{xi) <C for FT 1 < |(|,A)| < R and 5?A > -9, 

for any R > 0, and C = C(R) > sufficiently large and 9 = 9{R) > 
sufficiently small, where | ■ |^i is as defined in Proposition \4- 1\ 

Proof. This follows by compactness of the set of frequencies under consid- 
eration together with the fact that the resolvent (A — Lg) is analytic with 

respect to H 1 in (£, A) for A in the resolvent set p^ x (L) of with respect to 

H 1 once we establish that A lies in the /9^i(L|) whenever R~ l < |(£, A)| < R 

and 5RA > —9. Here, the .fP-resolvent set p^i(L^) is defined as the set 

of A such that (A — L|) has a bounded inverse with respect to the H 1 - 
operator norm [Kat, Yo, Pa], and similarly for p L 2(L^), to which we refer 
elsewhere simply as p(Li); see Definition IA.3| Appendix A. Analyticity in 
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A on the resolvent set of (A — L)^ 1 holds for general operators L (Lemma 
IA.51 Appendix A), while analyticity in £ on the resolvent set of (A — Lg) 
follows using standard properties of asymptotically constant-coefficient ordi- 
nary differential operators [He, Co, ZH] from analyticity in £ of the limiting, 
constant-coefficient symbol together with convergence at integrable rate of 
the coefficients of to their limits as x\ —* ±00. Alternatively, we may 
establish this directly for 9 sufficiently small, using the smoothing properties 
induced by (A1)-(A2), (H0)-(H2); see Exercises EU 1-2, Appendix A. 
Under assumptions (Al)-(A2), (H0)-(H2), p^ 1 (L^) and p L i{L£) agree 

on the set R~ l < |(£, A)| < R, 5?A > —9, provided 9 is chosen sufficiently 
small relative to R~ . For, standard considerations yield that, on the "do- 
main of consistent splitting" comprised of the intersection of the rightmost 
components of the resolvent sets of the limiting, constant-coefficient opera- 
tors as x\ — > ±00, the spectra of the asymptotically constant-coefficient 

ordinary differential operators with respect to H s , for any < s < q(d), 
q(d) as defined in (HO) (indeed, with respect to any reasonable norm), con- 
sist entirely of isolated L 2 eigenvalues, i.e., A — is Fredholm in each H s ; for 
further discussion, see [He, GZ, ZH, ZS, Z.3-4] or Section 5, below. By (K3), 
the domain of consistent splitting includes the domain under consideration 
for 9 sufficiently small, giving pgi^Lg) = piji{L^) as claimed. Alternatively, 
we may establish this by the direct but more special argument that (|4.35[) 
together with |(A - L^)~ l \ L 2^ < C\ implies |(A - -i^) -1 ^^) < C 2 . 

With this observation, the result follows by the strong spectral stability 
assumption (|2.12|) . since analyticity of the resolvent in (£, A) implies that 
{(£, A) : A G p L 2(L^)} is open, and therefore contains an open neighborhood 

of {(I, A) : KA > 0} \ {(0,0)}. (Recall, a L i(L) := p L 2{L) c .) □ 

Remarks 4.9. 1. As the argument of Proposition 14.81 suggests, uniform 
mid-frequency bounds are essentially automatic in the context of asymp- 
totically constant-coefficient ordinary differential operators, following from 
compactness/continuity of the resolvent alone. 

2. Resolvent bounds (|4.32l) for A with negative real part, correspond 
roughly to time-exponential decay in high-frequency modes, similarly as in 
(l4~2l . Indeed, a uniform bound |(A - ^-)~ 1 |^i < C for all f, KA > -9 < 0, 
or equivalently |(A — L)~ 1 \ H i < C for all 3?A < —9 < would imply time- 
exponential linearized stability |e L *|^i < C{9)e~ et for any 9 < 9, by a 
general Hilbert-space theorem of Priiss [Pr]; see Remark IA.23| Appendix 
A. Contrasting with estimate 1)3.20(1 corresponding to local well-posedness, 
what we have done here is to use the additional structure (A2), (HI) to push 
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7o into the negative half-plane. Note that we made no stability assumptions 
in the statement of Proposition 14.71 

3. High-frequency linearized resolvent bounds may alternatively be ob- 
tained by Kreiss symmetrizer techniques applying to a much more general 
class of systems; see [GMWZ.4], In particular, we may dispense with the 
symmetry assumptions in (A1)-(A2) connected with energy estimates using 
integration by parts, replacing them everywhere with sharp, spectral condi- 
tions of hyperbolicity/parabolicity. By corresponding pseudodifferential es- 
timates, it might be possible to recover also the nonlinear time-evolutionary 
estimate Q4.7|l for this more general class of equations. This would be a very 
interesting direction for future investigation. 

4.2.1. Splitting the solution operator. Letting L as usual denote 
the linearized operator defined in (|3.8jl . define a "low- frequency cutoff" 

Si(t)f(x) := 

(437) / / e^(X-L^f( Xl ,i)did\ 

of the solution operator e Lt described in (|3.12|) . similarly as in the constant- 
coefficient argument carried out in Appendix B, where 9±, 02 > are to 
be determined later. From Propositions I4.7H4.81 we obtain the following 
decompositions justifying this analogy. 

Corollary 4.10. Given (Al)-(A2), (HO)-(Hl), and strong spectral stabil- 
ity, let \e Lt \ < Ce 70 * denote the C° semigroup on L 2 generated by L, with 
domain V(L) := {U : U, LU € L 2 }. Then, for arbitrary 8 2 > 0, Q x > 
sufficiently small relative to 62, and all t > 0, 
(4.38) 
e Lt f(x) = S 1 (t)f(x) 

/-01 + too r 
/ I|||2 + |cj A |2|>0 +£ ) e^' 5+At (A - L|) _1 /(xi, |) d£d\ 
-61-iao JK. d - 1 111- 

for all f £ T>(L) n H 1 , where Ip denotes the indicator function for logical 
proposition P ( one for P true and zero for P false ) and f denotes Fourier 
transform of f . Likewise, for < t < T , 
(4.39) 

e L{t - s) f(s)ds = [ S 1 (t-s)f(s)ds 
Jo 

/ I||| 2+ | 3A | 2 | >ei+e2 e^- £+A '(A - L l )- 1 f'r{x 1 , 1 X) didX 

-di— too JM d_1 
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for all f € L 1 ([0, T]; 2?(L) ri-ff 1 (x)), where f denotes Laplace-Fourier trans- 
form of f and truncation f T is as in (|3. 14|) . 

Proof. Using analyticity on the resolvent set p(L) of the resolvent (A — L) -1 
( Lemma IA.5I Appendix A), the bound 

(4.40) | (A - L)- l f\ < C(|(A - L)- l \\Lf\ + \f\)\X\~ l 

coming from resolvent identity (|A.4|) (Lemma lA.7|) . and the results of Propo- 
sitions l4~TH4,81 we may deform the contour in (|3.22|) using Cauchy's Theorem 
to obtain 
(4.41) 

e Lt f(x) = I e xt (X - L)- 1 f{x) dX 

isRA=6»2-|QA| 2 ; |3A| 2 <6»i+6»2 
»— 0i+ioo 



I m ^\>B 1+ e 2 e Xt (X-L)- 1 f(x)dX 

-9i—ioo 

I f e*^ A *(A-L*)-V>i,£)<W 

JsRA=6» 2 -|3 , A| 2 ; |OA| 2 <9i+6»2 Jl&fi- 1 

^x?\>e 1+ e 2 / e^ +xt (X - L^) -1 f(xi,£) didX. 



Using exercise 14 . 1 1 1 b elow . we may exchange the order of integration to write 
the first term of the last line as 

(4.42) / / e^' £+A *(A — Lj!)~ x f{x\, (■) d£dX. 

J^eRd- 1 J5RA=e 2 -|SA| 2 ; |SA| 2 <6»i+e 2 

Deforming individual contours in the second integral, using analyticity of 
(A — Lj) on the resolvent set p(L^) together with the results of Propositions 
I4.7H4.81 and exchanging orders of integration using the result of Exercise 
I4.1H we may rewrite this further as 
(4.43) 

/ / e*~* +At (A-L^7>i,0<W 

J\i\ 2 <o 1 +8 2 Jm=e 2 -\z\\2~\i\2 s 

9 1 +iVeT+0 2 



J — 01 — i\J 01+02 



/ / S-^ x \X-L^f(xi,l)dldX 

J\i\ 2 <e x +e 2 Jka=02-|3A| 2 -||| 2 

r- 0i+V0i+02 r 

+ / / %A| 2 | + ||| 2 >0i + 02 e ^' +A< ( A - L |) _1 ^'^^ A - 



55 



Substituting back into Q4.41JI and recombining terms, we obtain (|4,38|) , A 
similar computation, together with the observation that 

(4.44) fT( Xl ,i,X):= [ T e- Xs f( Xl ,i,s)ds 

Jo 

is analytic in (£, A) as the absolutely convergent integral of an analytic in- 
tegrand, yields (|Q9|1 . □ 

Exercise 4.11. For g E H (x), show using Parseval's identity and the 
Fourier transform definition \g\H s (x) := |(1 + |£| 2 ) s ^ 2 5(£)Il 2 (§) that 

(4-45) f \9\h {xi) dt<C{l + L^)\g\l 1{x) 

J\(\>L 

converges uniformly to zero as L — > +oo ; where g(xi,£) denotes the Fourier 
transform of g with respect to x. 

Remark 4.12. Defining the "high-frequency cutoff" Si{t) := eLt ~ S\{t), 
we have jj e L ^ds = jj Si(t - s)ds + jj £ 2 (t - »)ds, where, by (lOfll) . 
(4.46) 

S 2 (t-s)f(s)ds = / S 1 (t-s)f(s)ds 



o 



/— f 1+100 /• _ ^ 

/ I^+Wi^+fl,^^^ - ^)~ 1 / T (^i, f, A) dfdA. 
-01— ioo JIR 1 * -1 

Note that Si(t) = — S 2 (t) ^ for £ < 0, in violation of causality. "Causality 
error" J. S\(t — s)f(s) ds is the price for fixing the cutoff T in the second 
term of (|4.46|) . 

4.2.3. Proof of the main estimate. Now, assume the following low- 
frequency estimate of [Z.3], to be discussed further in Section 5. Similar 
results hold for profiles of nonclassical, over- or undercompressive type, as 
described in [Z.3-4]. 

Proposition 4.13 ([Z.3] 10 ). Under the assumptions of Proposition 
\Si(t)(h + 9 x f2)\ L 2 ix) < C(l + t)- (d - 1)/4+ 1/i| L i (a; ) 

( ' } +c(i + t)-^y^-^\ f2 \ LHx) 

for all t > 0, for any fixed e > 0, for some C = C(e) > sufficiently large, 
for any f\ € L x (x) and f 2 G L x (x). In i/ie case of strong (inviscid, not 
refined) dynamical stability, we may take e = and C > fixed. 
10 For a more detailed exposition, see [Z.4]. 
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Corollary 4.14. Under the assumptions of Proposition 

(4.48) V(x,t):= [ Si{T -t)f(x,s)ds 

Jo 

satisfies 
(4.49) 

j\-se l( T-t )lv{t)l 2 L2{x) dt < C ^J\ X + T _ sr (^i)/4 +£ -i/ 2|/(s)|Li dfi 

/or ffl/ZT > 0, /or any f(x,t) G ^([0, T]; L 1 ^)). 
Proof. Using (|4.47|) and the inequality 

(4.50) e^^l + t - s )-(d-l)/4+e-l/2 < C(l + T - 8 )-{d-i)/4+*-l/2 

for < s < t < T, we obtain 
(4-51) 

< C £ e^-*) ( + t - S )- (d - 1)/4+e - 1/2 |/( S )| L1( ,)d S ) 2 



T 



<c 2 / (l + r-^w^^i/wi^cfa 



Lemma 4.15. Under the assumptions of Proposition 



□ 



(4-52) |5i(*)/)|i2(*) < Ce-^|/I 

for all t < 0, for any f € H (x). 

Proof. Immediate, from representation (|4.37|) and Proposition 14.81 □ 
Corollary 4.16. Under the assumptions of Proposition 

(4.53) W!(x,t):=J Si(T -t)f(x,s)ds 
satisfies 

(4.54) £ e-»W-V\WMl Hx) dt < C j\-^ T ^\f(t)f Hl{x) dt 
for allT> 0, for any f(x,t) e L 2 ([0, T],H 1 (x)). 
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Proof. Using the triangle inequality, and (|4.52|) . we obtain 
(4.55) 



-89i (T-t) 



S 1 (t-s)f(s)ds 



L 2 (x) 



(It 



< c 



< c 



< Co 



T 



T 



-46>i (T-t) 



- 2fll(T - f) |/( S )l| 1( ,)« 



- 2ei{T ~ t] \m\ 2 H Hx) d S . 



□ 



Lemma 4.17. Under the assumptions of Proposition \J78[ 
(4.56) 
W 2 (x,t) := 

9i+ioo 



P.V. 

satisfies 



! I l5l 2 +|3A| 2 |>e 1 +e 2 ' : 



(A-L f -)-V T (xi,e,A)d^A 



(4-57) j Q e-^ T -*\W 2 (t)\l Hx) dt<cj o e-^ T -*\f(t)f Hl(x) dt 

for allT> 0, for any f{x,t) G L 2 ([0, T],H 1 (x)). 
Proof. Using the relation 

(4.58) Ch(X) = F(e-^ xt h) (-3FA) 

between Fourier transform Th and Laplace transform Ch of a function /i(t), 
we have 



(4.59) Fe eit W 2 {t){-k) = LW{-Q X + ik) = (A - L)" 1 /:/(-0i + ifc). 
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Together with Parseval's identity and Propositions I4.7h4.8l this yields 
^1^)1^) dt = \e^W 2 (t)\l 2(x ^ T]) = 



T e M " l \W 2( 



o 



I |?| 2 + |SA| 2 |>6I 1 +6I 2 ( A ~~ L i) f T ~ d ^ 

-6\ ~ioc 



(4.60) 



L*( X1 ,$ 



i— too 

T 



c [ ^\f{t)^ Hl{x) dt 

yielding the claimed estimate upon multiplication by e~ 2dlT . □ 

Lemma 4.18. Under the assumptions of Proposition 
(4.61) 

W (x,t) : = 

/ %| 2+ | 3A | 2 |> ei+ , 2 e^ +A4 (A - L^Uoixu^didX 

satisfies 

(4-62) £ e- 2e ^\Wo(t)\h {x) dt < Ce-^ T (\U \ 2 H1{X) + \LU \ 2 m ) 
for all T>0, for U , LU G H\x, [0, T]). 

Proof. Without loss of generality taking ^ 0, we have by (j4,58|) . Parse- 
val's identity, Propositions 14, 7^4.81 and the resolvent identity 

(4.63) (A - L^Uo = A" 1 ((A - L^Lfo + U ) 

( (|A.4|) . Appendix A), that 

e 2 ^\Wo(t)\ 2 L2{x) dt = \e^W,(t)\l Hx ^ = 

-Oi+ioo 2 

d\ 



o 



(4.64) J-Ox-ioo 

< c 



I ICI 2 +|SA|2|>e 1 +e 2 ( A ~ L l) u c 

-9i —ioo 



L2 ( Xl ,0 



il — too 



\M- 2 dX{\Lf \ Al{xii) + \U \ LZ{xi)i) ) 

< C(\LU \h1(x) + \Uo\l*(x)), 
from which the claimed estimate follows upon multiplication by e~ 2SlT . □ 
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Proof of Proposition ^ -5\ Immediate, combining the estimates of Proposi- 
tion Corollaries 14.141 and 14.161 and Lemmas 14.181 and 14.171 and using 
representation 
(4.65) 

U(T) = e LT U + f e L ^ l \h + d x f 2 ) dt =: Sx(T)U + V + W + W 1 + W 2 
Jo 

coming from jOity - jOty . □ 

4.3. Nonlinear stability argument. We are now ready to establish 
the result of Theorem 12.191 assuming for the moment the bounds asserted 
in Proposition 14. 13l on the low- frequency solution operator S\(t). Define the 
nonlinear perturbation 

(4.66) U:=U-tJ, 

where U denotes a solution of (jl.2|) with initial data Uq close to U. Then, 
Taylor expanding about U, we may rewrite (|l.'2j) as 

(4.67) U t -LU = d x Q(U,d x U)=:d x f 2 , [7(0) = U , 
where 

\Q(u,d x u\ <c(\u\\d x u\ + \u\ 2 ), 

(4.68) \a x Q(u,a x u\ < c(\u\\d 2 x u\ + \a x u\ 2 + \u\\a x u\), 
\a 2 x Q(u,a x u\ < c(\u\\alu\ + \a x u\\a 2 u\ + \a x u\ 2 ), 

so long as U remains uniformly bounded. 

Lemma 4.19. So long as U remains uniformly bounded, 

\h{t)W {x) <c\u{t)\ 2 HHx) , 

\d x f2(t)\H 1 (x) < C\U(t)\ 2 H[d/2]+2(x y 

Proof. Immediate, by Q4.68JI and \d x U\ L oo^ < \U\ H [d/2]+2/ x y □ 

Proof of Theorem M.ltA Set 8#i := 9 2 , 9 2 > as in the statement of Propo- 
sition aim define 

(4.70) C(t):= sup |*7(r)| L 2(l + T)^- e , 

0<T<t 

where e > is as in the statement of Theorem 12. 191 By local well-posedness 
in H s , Proposition ll.lbl and the standard principle of continuation, there 
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exists a solution U G H s (x), U% G H s ~ 2 (x) C L 2 (x) on the open time- 
interval for which |i7|if« remains bounded. On this interval, £ is well-defined 
and continuous. 

Now, let [0, T) be the maximal interval on which \ U\fjs^ remains bounded 
by some fixed, sufficiently small constant 5 > 0. By Proposition 14.11 

\U{t)\%. <C\U(0)\ 2 Hs e- et + C f e- e ^\U(r)\l 2 dr 
(4.71) 70 

< C 2 {\U(0)\h +CW 2 )(1 + r)-( d - 1 )/ 2+2£ . 

Combining (|4.31f) with the bounds of Lemma 14.191 we thus obtain 
(4.72) 

fe- m - T) \U(T)\l 2 dr < C(l + t)-( <J - 1 )/ 4fe |L/bl!i 
Jo 

+ Ce- 2ei *(|C/o| L 2 + |LC/ |HO 

+ C f e- 2e ^- T) \d x h{r)\ 2 Hl dr 
Jo 

+ C(j\l + t - r)-^ 1 )/^-^!^^)!^ drf 

< C (i + t)-^/ 2 ^\u -u\ 2 LlnH3 

+ C((t) A ( J\l + t - r )-( d - 1 )/ 4 +- 1 /2 (1 + r) -(d-l)/2+2e dr 

< C 2 (|C7 - C7|i lnH3 + C C (t) 4 ) (1 + r)-^/ 2 ^ 



for d > 2 and e = 0, or d > 3 and e sufficiently small (exercise). 
Applying Proposition 14. II a second time, we obtain 

\U(t)\ 2 HS <Ce~ e ^\U -U\ 2 H s+C f e- e ^\U{r)\ 2 L2 dr 

(4.73) Jo 

< C 2 (\U - U\ 2 LlnHs + <7C(i) 4 )(l + rr (d - 1)/2+2e 

and therefore 

(4.74) ((t) < C(\U - U\ LlnHs + CC{t)'- 



so long as |J7|h s an d thus |t7| remains uniformly bounded. 

From ()4.74|) . it follows by continuous induction (exercise) that 



(4.75) C(t) <2C\U(0)\ 



L 1 nH s 
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for |f/(0)|/,i n #s sufficiently small, and thus also 

(4.76) \U(t)\ H s < 2C(l + t)^ d - 1 ^\U(0)\ LlnHs 

as claimed, on the maximal interval [0, T) for which \U\jj a < S. In particular, 
\U(t)\H« < 5/2 for \U(0)\ L i nH s sufficiently small, so that T = +oo, and we 
obtain both global existence and H s decay at the claimed rate. Applying 
1)1.13(1 . we obtain the same bound for \U(t)\i,°°, and thus, by LP interpolation 
(formula (lETIi Appendix B), for \U(t)\ LP , all 2 < p < oo. □ 

Together with the prior results of [Z.3], the results of this section com- 
plete the analysis. We'll review those prior results in the low-frequency 
analysis of the next section. 

Remarks 4.20. 1. Sharp rates of decay in L p , 2 < p < oo may be obtained 
by a further, bootstrap argument, as in the related analysis of [Z.4]. 

2. The approach followed here, yielding high-frequency estimates en- 
tirely from resolvent bounds, also greatly simplifies the analysis of the one- 
dimensional case, replacing the detailed pointwise high-frequency estimates 
of [MZ.3]. Pointwise low-frequency bounds for the moment are still required; 
however, a new method of shock-tracking introduced in [GMWZ.3] might 
yield an alternative approach based on resolvent bounds. This would be an 
interesting direction for further investigation. At the same time, the new 
approach yields slightly more general results: specifically, we may drop the 
assumption that the hyperbolic convection matrix A*{U(x)) have constant 
multiplicity with respect to x, and we may reduce the regularity requirement 
on the coefficients from C 5 to C 3 and the regularity of the perturbation from 
L 1 n H 3 to L 1 n H 2 . 11 

3. Kawashima-type energy estimates are available also in the "dual", re- 
laxation case, at least for small-amplitude shocks. Thus, we immediately ob- 
tain by the methods of this article, together with prior results of [Z.3], a cor- 
responding result of nonlinear multi-dimensional stability of small-amplitude 
relaxation fronts. In the one-dimensional case, an energy estimate has re- 
cently been established for large-amplitude profiles in [MZ.5]; in combination 
with the methods of this article, this recovers the large-amplitude nonlin- 
ear stability result of [MZ.5] without the restrictive hypothesis of constant 

11 More precisely, to C 4 and L 1 n H 3 by the arguments of this article alone. To obtain 
the weaker regularity stated, we must substitute for the uniform H 1 — » H 1 high- and mid- 
frequency resolvent bounds of Propositions l4.7l and l4.8l the corresponding L 2 — » L 2 bounds 
obtained in [MaZ.3] by pointwise methods or in [GMWZ.4] by symmetrizer estimates. 
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multiplicity of relaxation characteristics. Besides being mathematically ap- 
pealing, this improvement is important for applications to moment-closure 
models. 

4. The analysis so far, with the exception of the large-amplitude version 
of the "magic" energy estimate (|4.2|) . is all "soft". The novelty lies, rather, 
in the argument structure, which for the first time successfully integrates 
Parseval- and semigroup-type (i.e., direct, pointwise-in-time) bounds in the 
analysis of long-time viscous shock stability. See [KK] for an interesting ear- 
lier approach based entirely on Parseval's identity and Hausdorff- Young's 
inequality, which may be used in cases of sufficiently fast decay to establish 
Jo°° l^( s )li 2 ^ s < i n the shock-wave setting, this occurs for dimensions 
d = 3 and higher [GMWZ.2]. For further discussion/comparison, see Re- 
mark [231 below. More generally, our mid- and high-frequency analysis could 
be viewed as addressing the larger problem of obtaining time-algebraic de- 
cay estimates from spectral information in the absence of either a spectral 
gap or sectorial-type resolvent bounds; see Remark IA.23I 

Note in particular that our conclusions are independent of the particu- 
lar setting considered here, depending only on the availability of nonlinear 
smoothing-type energy estimates favorable in the highest derivative. Thus, 
our methods may be of use in the study of delicate stability phenomena 
arising in other equations, for example, stability of traveling pulse- or front- 
type solutions of nonlinear Schrodinger's equation. Moreover, they are truly 
multi-dimensional, in the sense that they do not intrinsically depend on pla- 
nar structure/decoupling of Fourier modes, so could be applied also to the 
case of a background wave u varying also in the x directions: for example 
to traveling waves in a channel (xi,x) Eixl,lc R ' -1 bounded. 

Indeed, there is no requirement that the resolvent equation even be 
posable as an (possibly infinite-dimensional) ODE. Thus our methods may 
be of use in interesting situations for which this standard assumption does 
not hold: in particular, stability of irrational- speed semi- discrete traveling- 
waves for non-upwind schemes (see discussion, treatment of the upwind 
case in [B, BHR]), or stability of shock profiles for Boltzmann equations. Of 
course, in each of the mentioned cases (fully multi-dimensional, semidiscrete, 
and kinetic equations), there remains the problem of carrying out a suitable 
low-frequency linearized analysis; nonetheless, it is a substantial reduction 
of the problem. 



63 



5 Low frequency analysis/completion of the proofs 



We complete our analysis by carrying out the remaining, low-frequency anal- 
ysis, establishing Propositions 14, 13*1 and |2~"T*"1 and Theorem 12. 181 Our main 
tools in this endeavor are the conjugation lemma, Lemma 13.131 by which 
we reduce the resolvent equation to a constant-coefficient equation with im- 
plicitly determined transmission condition, and the Evans function, through 
which we in effect obtain asymptotics for the resulting, unknown transmis- 
sion condition in the low- frequency limit (Proposition 12. 14*)) . 

5.1. Bounds on Si(t). Continuing in our backwards fashion, we begin 
by establishing the deferred low-frequency bounds cited in Proposition 14. 131 
of the previous section under appropriate resolvent bounds, thus reducing 
the remaining analysis to a detailed study of the resolvent and eigenvalue 
equations, to be carried out in the remainder of this section. The following 
low-frequency estimates were obtained in [Z.3] by pointwise estimates on 
the resolvent kernel; we sketch a different proof in Section 5.4, based on 
degenerate symmetrizer estimates as in [GMWZ.l, W]. 

Proposition 5.1 ([Z.3-4]). Assuming (Al)-(AS), (H0)-(H5), for p := 
|(£,A)| sufficiently small, and 

m = -0(|£| 2 + |9A| 2 ) 

for 9 sufficiently small, there holds 

(5.1) Ar^/U^o < C 7 i72P (1 - q)I/3| - 1 |/Il 1 (x,) 

for all 2 < p < oo, < \j3\ < 1, where 



(5.2) 7i(£, A) := 



1 in case of strong (uniform) dynamical stability, 
1 + J^jiP 1 I^A — iTj(£)\ + p]^ 1 otherwise, 



(5.3) 72 (l\) : =l + Y,[p- 1 \^-vt(0\+P}- t , 0<t<l, 

r}j(-) is as in (H5) and the zero-level set o/A(^,ir) for £, r real is given by 
Ujg&iTjd)), and 

I for Lax or overcompressive case, 

(5.4) a := < 

I 1 for undercompressive case. 
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(For a definition of overcompressive and undercompressive types, see, e.g., 
Section 1.2, [Z.4]-) More precisely, t := 1 — 1/K max , where K m&x := 
maxK^ 1 = maxs^ is the maximum among the orders of all branch sin- 
gularities rjj (•), and rjj defined as in (H5); in particular, t = 1/2 in the 
(generic) case that only square-root singularities occur. 

Remark 5.2. Rewriting the L 1 — > L 2 resolvent bound IB~T1 as 

we may recognize it as essentially a second-order correction of the corre- 
sponding 1? — > L 2 resolvent bound 

(5-6) |(^--A)-Vl^o< C|/ ^ (ai) 

of the inviscid theory, which in turn is approximately the condition for a 
bounded C° semigroup (Appendix A). The singular factor 72 appearing in 
the numerator of (|5.5|) . and the square root in the denominator, are new 
effects connected with the fact that the bound is taken between L 1 and I? 
rather than in L 2 alone; for further discussion, see [Z3]. 

Proof of Proposition \4-13[ For simplicity, restrict to the Lax-type, uniformly 
dynamically stable case; other cases are similar. By analyticity of the re- 
solvent on the resolvent set, we may deform the contour in Q4.13j) to obtain 

(5.7) 

Si(t)f(x) = 

[ I „ e*^ +A *(A - L~)~ 1 f{xi,l) d£d\, 

J\i\ 2 <6i+02 ./SA=,=£fc(|{[ 2 -|9Ap)>-0i 

where 8j > are as in (|4.37|) . and Corollary I4.1U1 with 9\ sufficiently small 
in relation to 62- Bounding 

( 5 - 8 ) I / U°° (si)) < \f\v(x 1 ,x) = l/|l 

using Hausdorff- Young's inequality, and appealing to the L 1 — ► L p resolvent 
estimates of Proposition 15.11 with 6 := s ^g 2 > taken suffiently small, we 
may thus bound 

(5-9) Wxi,£,\)\ L 2 ixi) < \f\ib(£,\), 



65 



for u := (Lg — A)/ and p := |£| + |A| sufficiently small, where 

(5.10) Ki\):=p- 1 l2 (i\). 
Denoting by r(£) the arc 



(5.11) 



KA 



(\C\ 2 + |QA| 2 ) > -e x 



h + o 2 



and using in turn Parseval's identity, Fubini's Theorem, the triangle inequal- 
ity, and our L 1 — > L 2 resolvent bounds, we may estimate 



(5.12) 

\ S l(t)f\ L 2^ Xl ^(t) 



e xt u(x 1 ,(,,\)d\ 



*er(Q 



2 „ 



e Ai «(aci,C,A)<ZA 



< 



1/2 



1/2 



L 2 (*i) 



w(a;i,f, \)\ L 2( Xl )d\ 



1/2 



<l/h 



e KAt 6(£, A)dA 



Asr(£) 



1/2 



from which we readily obtain the claimed bound on |Si(i)/|2 using (|5.1U|) 
and the definition of 72. Specifically , parametrizing r(£) by 

X(i,k) = ik-6 1 (k 2 + \i\ 2 ), k£R, 

and observing that in nonpolar coordinates 



9 72 



(5.13) 



< 



< 



(N+iiir'a+E 



l*-Ti(0K^-i 



where e := — - — - (0 < e < 1 chosen arbitrarily if there are no singularities), 
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we obtain a contribution bounded by 



(5.14) 



C\f\i 
<C\f\i 



-^+\l?)t {p) -i l2dk 



2 dt /2 



-20|£| 2 tj||-2e 



+00 



e -0\k\H lkl e-l dk 



2 .U/2 



-26\i\H^-2e 



+00 



-e\k\H 



\k-r 3 (Or l dk 



2 _\ 1/2 



<C|/| 



< ci/iit-^ 1 )/ 4 , 



-20|f| 2 ij£|-2e 



+00 



2 .U/2 



yielding the asserted bound for t > 1; for i < 1 on the other hand, we obtain 
the asserted uniform bound by local integrability together with boundedness 
of the region of integration. Derivative bounds, (3 = 1, follow similarly. □ 

Remark 5.3. In the proof of Proposition ^. l'A\ we have used in a fundamen- 
tal way the semigroup representation and the autonomy of the underlying 
equations, specifically in the parabolic deformation of contours revealing 
the stabilizing effect of diffusion. The resulting bounds are not available 
through Parseval's identity or Hausdorff- Young's inequality, correspond- 
ing to the choice of straight-line contours parallel to the imaginary axis. 
Such "parabolic," semigroup-type estimates do not immediately translate to 
the small-viscosity context, for which the linearized equations are variable- 
coefficient in time. Moreover, it is not clear in this context how one could 
obtain resolvent bounds between different norms, since standard pseudo- or 
paradifferential techniques are based on Fourier decompositions respecting 
L 2 but not necessarily other norms. The efficient accounting of diffusive ef- 
fects in this setting is an interesting and fundamental problem that appears 
to be of general mathematical interest [MeZ.l]. 

Remark 5.4. As noted in the acknowledgements, all estimates on the 
linearized solution operator, both on S\{t) and S2 (t), have been obtained 
through resolvent bounds alone. These could be obtained in principle by 
a variety of methods, an observation that may be useful in the more gen- 
eral situations (e.g., semidiscrete or Boltzmann shock profiles) described in 
Remark OI34. 
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5.2. Link to the inviscid case. It remains to study the resolvent 
(resp. eigenvalue) equation 



(5.15) 

where 
(5.16) 




L U 



(L i - X)U = (B 11 ^)' - {A l U)' -i^2A j CjU 

+ i B^U' + i Y,(B lk CkU)' - J2 B ik S&U - XU. 

j+\ k^l j,te£\ 

Up to this point our analysis has been rather general; from here on, we make 
extensive use of the property, convenient for the application of both Evans 
function and Kreiss symmetrizer techniques, that the (|5.15|) may be written 
as an ODE 



(5.17) W' -A(i,X,xi)W 



F, 
0. 



in the phase variable W := ^n^/j £ C n+r , where F ~ /. This is a 

straightforward consequence of block-diagonal structure, (Al), and invert- 
ibility of A\±, (HI); see [MZ.3, Z.4] for further details. (There is ample 
reason to relax this condition, however; see Remark 14.201 4.) Our first task, 
carried out in this subsection, is to make contact with the inviscid case. 

5.2.1. Normal modes. By the conjugation lemma, Lemma i'A.l'M we 
may reduce ()5.16|) by a change of coordinates W = P±Z±, f = P±f, with 
P± — > / as x\ — > ±oo, P analytic in (£, A), to a pair of constant-coefficient 
equations 



(5.18) Z' ± -A±&X)Z± 



f 
0. 



on the half-lines x\ ^ 0, coupled by the implicitly determined transmission 
conditions 

(5.19) P+Z + (0) - P_Z_(0) = 
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at the boundary x\ = 0, where 



(5.20) 



A±(£,A):= lim A(£,A,xi) 

X\— >±oo 



that is, a system that at least superficially resembles that arising in the 
inviscid theory. The following lemma generalizes a standard result of Hersch 
[H] in the inviscid theory; see Exercise 12.101 

Lemma 5.5 ([ZS, Z.3]). Assuming (Al)-(A2), (HI) (and, implicitly, ex- 
istence of a profile), the matrices A± have no center subspace on 



(5.21) 



- ' '^ lo ' 'cv\|2 J ' 



i + iei 2 + i^A| 



for 8 > sufficiently small; moreover, the dimensions of their stable and 
unstable subspaces agree, summing to full rank n + r. 

Proof. The fundamental modes of (|5.18|) are of form e^V, where (jl, V 
satisfy the characteristic equation 



(5.22) 



v = 0. 



The existence of a center manifold thus corresponds with existence of solu- 
tions fj, = V of (|5.22|) . £i real, i.e., solutions of the dispersion relation 



(5.23) 



(" B ± ~ 1 E ~ A/)v = 0. 

j,k j 

But, A G a(-B& - iA%) implies, by (|3~5Uj) . that 

KA<-^i^i 2 /(i + iei 2 ). 

Noting that for low frequencies |3A| = we thus obtain 

-e(||p + |3fA| 2 ) 



5ftA < 



1 + |£| 2 + |9A! 2 ' 
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in contradiction of 1)5.21(1 . 

Nonexistence of a center manifold, together with connectivity of the set 
described in (|5.21|) . implies that the dimensions of stable and unstable man- 
ifolds at +00/— 00 are constant on A. Taking A — > +00 along the real axis, 
with £ = 0, we find that these dimensions sum to the full dimension n + r 
as claimed. For, Fourier expansion about £1 = 00 of the one-dimensional 
(£ = 0) dispersion relation (see, e.g., Appendix A. 4 [Z.4]) yields n — r "hy- 
perbolic" modes 

= -itiiaj + • . • , j = 1, • • • , n - r, 

where a*j denote the eigenvalues of A\ := A\-^ — (b\ A\ 2 and r "parabolic 1 

modes 

\i~r-+j = + • • • , j = 1, • • • r, 

where bj denote the eigenvalues of b^ 1 ] here, we have suppressed the ± indices 
for readability. Inverting these relationships to solve for fj, := we find, for 
A — > 00, that there are n — r hyperbolic roots fij ~ — A/a*;, and r parabolic 
roots fJ, n _ r+ j ~ \J A/6j. By assumption (HI), &<diA\{x\) 7^ for all x\, and 
so the former yield a fixed number k/(n — r — k) of stable/unstable roots, 
independent of x±, and thus of ±. Likewise, (Hl(i) implies that the latter 
yields r stable, r unstable roots. (Note: we have here used the existence of a 
connecting profile: i.e., equality is not a consequence of algebraic structure 
alone.) Combining, we find the desired consistent splitting, with {k + r)/ 
(n — k) stable/unstable roots at both ±00. □ 

Corollary 5.6. Assuming (Al)-(A2) and (HI), on the set A defined in 
(|5.21|) . there exists an analytic choice of basis vectors Vj^ , j = 1, . . . k (resp. 
V- , j = k + 1, . . . ,n + r) spanning the stable (resp. unstable) subspace of 
A±. 

Proof. By spectral separation of stable and unstable subspaces, the associ- 
ated (group) eigenprojections are analytic. The existence of analytic bases 
then follows by a standard result of Kato; see [Kat], pp. 99-102. □ 

Definition 5.7. For x\ ^ 0, define normal modes for the variable-coefficient 
ODE (gmj) as 

(5.24) Wf ( Xl ) ■.= P±e A±xl V 3 ± . 

In particular, W^, . . . , span the manifold of solutions of the eigen- 
value equation decaying as x\ — > +00, and W^~ +1 , . . . , W^ +r span the mani- 
fold of solutions of the eigenvalue equation decaying as x\ — > -co. 
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5.2.2. Low-frequency asymptotics. For comparison with the invis- 
cid case, it is convenient to introduce polar coordinates 

(5.25) (£,A)=:(p£o,pAo), 

p G M 1 , (£o,A ) G M^ 1 x {5ft A > 0} \ {(0,0)}, and consider Wp, Vf as 
functions of (p, £q,\o). With this notation, the resolvent (resp. eigenvalue) 
equation (|5.16|) may be viewed as a singular perturbation as p — > of the 
corresponding inviscid equation. The following results quantify this obser- 
vation, separating normal modes into slow-decaying modes asymptotic to 
those of the inviscid theory (the "slow manifold" in the singular perturba- 
tion limit) and fast-decaying transient modes associated with the viscous 
regularization. 

Lemma 5.8 ([K, Me. 4]). Let there hold (Al), (H2), and (H4), or, more 
generally, o-(df^{u±)) real, semisimple, and of constant multiplicity for £ G 
R d \{0} anddetdf 1 ^ 0. Then, the vectors {rf ,•• • , r"^}, {rf + +1 , • • ■ ,r+} 

defined as in (j2.6|) - ()2.8|) . Section 2.2 (and thus also A(£,A)J may be chosen 
to be homogeneous degree zero (resp. one), analytic on £ G M^ -1 , 5ftA > 
and continuous at the boundary £ G M.^ 1 \ {0}, KA = 0. 

Proof. See Exercises 4.23-4.24 and Remark 4.25, Section 4.5.2 of [Z.3] for a 
proof in the general case (three alternative proofs, based respectively on [K], 
[CP], and [ZS]). In the case of main interest, when (H5) holds as well, this 
will be established through the explicit computations of Section 5.4. □ 

Lemma 5.9 ([ZS, Z.3, MeZ.2]). Under assumptions (A1)-(A2) and (H0)- 
(H4), the functions Vj may be chosen within groups ofr "fast", or "viscous" 
modes bounded away from the center subspace of coefficient A± , analytic in 
G°)£o) Ao) for p > 0, £o G IftAo > 0, and n "slow", or "inviscid" modes 

approaching the center subspace as p — > 0, analytic in (p,£,o, Xo) for p > 0, 
£o G !ftAo > and continuous at the boundary p = 0, with limits 

(5.26) ^(0,£o,Ao) = ( ( ^ rlr | (l °' Ao) )' 
rf defined as in (|2~fil) - (t2~Sj) . 

Proof. We carry out here the simpler case 5?Ao > 0. In the case of interest, 
that (H5) also holds, the case 3?A = follows by the detailed computations 
in Section 5.4; see Remark 15.171 For the general case (without (H5)), see 
[Me. 2]. 
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Substituting U = e /iXl v into the limiting eigenvalue equations written in 
polar coordinates, we obtain the polar characteristic equation, 



(5.27) 



- (i P & Si +p 2 Y1 BJk ^ + p xi )} v = °> 



where for notational convenience we have dropped subscripts from the fixed 
parameters £o> Ao- At p = 0, this simplifies to 

(p'-Bi 1 - pA l ± ) v = 0, 

which, by the analysis in Appendix A. 2 of the linearized traveling-wave 
ordinary differential equation [pB± — A±) v = 0, has n roots p = 0, and 
r roots $tp ^ 0. The latter, "fast" roots correspond to stable and unstable 
subspaces, which extend analytically as claimed by their spectral separation 
from other modes; thus, we need only focus on the bifurcation as p varies 
near zero of the n-dimensional center manifold associated with "slow" roots 
p = 0. 

Positing a first-order Taylor expansion 



(5.28) 



p = + p 1 p + o(p), 
v = v° + v 1 /? + o(p), 



and matching terms of order p in (|5.27|) , we obtain 
(5.29) {~^A l ± - i A j ij - AJ)v° = 0, 

or equivalently —p\ is an eigenvalue of (A 1 )^ 1 (X + iA^) with associated 
eigenvector v°. 

For 3?A > 0, (A 1 )~ 1 (X + iA*) has no center subspace. For, substituting 
p} = i£i in (|5.29|) . we obtain A G a(iA*), pure imaginary, a contradiction. 
Thus, the stable/unstable spectrum splits to first order, and we obtain the 
desired analytic extension by standard matrix perturbation theory, though 
not in fact the analyticity of individual eigenvalues p. □ 

Remark 5.10. The first-order approximation (|5.29f) is exactly the matrix 
perturbation problem arising in the inviscid theory [K, Me.l]. 
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Corollary 5.11. Denoting W = {U^b^u 1 + b^u ) as above, we may ar- 
range at p = that all Wj satisfy the linearized traveling-wave ODE 

(B u u r y - (A^y = o, 

with constant of integration 



(5.30) B n uf ~ Al Uf 



0, for fast modes, 
for slow modes, 



r • as above, with fast modes analytic at the p = boundary and independent 
of (£o, Ao) for p = 0, and slow modes continuous at p = 0. 

Proof. Immediate □ 

5.2.4. The Evans function and its low-frequency limit. We now 

complete the analogy with the inviscid case, introducing the Evans function 
D and establishing its relation to the Lopatinski determinant A in the limit 
as frequency goes to zero, Proposition 12.141 

Definition 5.12. Following the standard construction of, e.g., [E.l-5, AGJ, 
PW, GZ, ZS], we define on the set A an Evans function 



D(Z, A) := det [W+,..., W+, W£ +v . . . , W„ 

(5.31) 



\x=0,\ 

det (p + y+, . . . , P+V+, P-V k+1 , ...,P-V^ 



|x=0,A 



measuring the (solid) angle of intersection between the manifolds of solu- 
tions of the eigenvalue equation decaying as x\ — > +oo and x\ — > — oo, 
respectively, where P±, Vp, Wj are as in Definition 15.71 

Evidently, D vanishes if and only if there exists an exponentially decay- 
ing solution of the eigenvalue equation, i.e., A is an eigenvalue of L?, or, 
equivalently, there exists a solution of the boundary- value problem (|5.18|) 
satisfying boundary condition (|5.19|) . That is, D is the analog for the vis- 
cous problem of the Lopatinski determinant A for the inviscid one. 



Proof of Proposition \2.1J\ We carry out the proof in the Lax case only. The 
proofs in the under- and overcompressive cases are quite similar; see [ZS]. We 
are free to make any analytic choice of bases, and any nonsingular choice of 
coordinates, since these affect the Evans function only up to a nonvanishing 



73 



analytic multiplier which does not affect the result. Choose bases Wj as in 
Lemma EHTl W = (U,z' 2 ), z 2 = &}V + b\ l u n . Noting that L Q u' = 0, by 
translation invariance, we have that u 1 lies in both Span {U^, • • • , U^} and 
Span (Uft +1 , • • • , U~ +r ) for p = 0, hence without loss of generality 

(5.32) U+ = U~ +r = u', 

independent of £, A. (Here, as usual, ""' denotes d/dx\). 
More generally, we order the bases so that 

W+, ...,W+ and W~ +r _ k _ v W~ +r 

are fast modes (decaying for p = 0) and 

W+ +1 , . . . , W+ and W« +v . . . , W~ +M 

are slow modes (asymptotically constant for p = 0), fast modes analytic and 
slow modes continuous at p = f Corollary I5.11|) . 

Using the fact that and U~ +r are analytic, we may express 

C/ 1 + (p) = ^ 1 + (0) + t/ 1 > + O ( /) ), 

^n+r (P) = ^n+r (°) + ^n+r.pP + °(p) > 

Writing out the eigenvalue equation 

(fl u u/) = (AV) - *> J^fl^&V 

(5.34) - ipQ2 B lk i k w)' + ip^ A ^j w + P Xw 

k^i j^i 

in polar coordinates, we find that Y + := U^ p (0) and Y~ := U~ +r p (0) satisfy 
the variational equations 



(B u Y'y = {A X Y') - \iJ2B jl Cju' 



(5.35) J " > 
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with boundary conditions Y + (+oo) = Y (— oo) = 0. Integrating from +00, 
—00 respectively, we obtain therefore 

B U Y±' - A 1 Y ± ' = ift(u) - iB^(u)u' 

( 5 - 36 ) « r * 

-iB^(u)u' + Xu- ifH u ±) + Xu ± 

hence Y := (Y~ — Y + ) satisfies 
(5.37) B U Y' - A 1 ? = i[ft] + \[u}. 

/A 1 A 1 \ 

By (Al) together with (HI), ( ^11 £,11 J * s i nver tible, hence 

(5 3g) (u, 4) - te, -zi, 4 + (Ax - bf\A\ 2 - bf')u) 

= (z 2 ,B 11 U' -AU) 
is a nonsingular coordinate change, where \) '■= ( 7 11 ^ii 2 ) U. 



z 2 ) ■ v^i b 2 

Fixing £q, Aq, and using Wi_ + (0) = W~ +r (0), we have 



D{p) = det (W+(0) + pW+(0) + o(p), ■■■ , W+(0) + o(l), 

Wk +1 (0) + o(l), • • • , W- +r (0)pW- +rp (0) + o(p) 

= det (V+(0) + pW+{0) + o(p), ■■■ , W+(0) + o(l), 
W K+1 (0) + o(l),-- - ,pY(0)+o(p)) 

= det (V+(0) • • • , W+(Q) W£ +1 (0) • • • , pF(0)) + o(p) 

Applying now (|5.38j) and using (|5.3fl|) and (|5.37|) . we obtain 
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fast 



D(p) = Cdet 



2,1 

o, 



i++l> 



slow 



z. 



fast 

^^^^^^^ s 

2,n+r-k-V ' Z 2,n+r-l> 

0,--- ,0, 



i[/«(n)]+A[n] 



+ o(p) 



U 1= o 



= 7 A(e,A) + ( P ) 

as claimed, where 

C det fz^i, • 



7 := 



' Z 2,fc' Z 2,n+r-fc-l' ' 



' z 2,n+r-l ) 



Noting that {zj v ■ ■ ■ , zj k } and {% n+r _ fc _ 1 , • • • , 2 2 n+r } span the tangent 
manifolds at n(-) of the stable/unstable manifolds of traveling wave ODE 



(|3.2I) at U+/U-, respectively, with 



IV 



u m common, 



we see that 7 indeed measures transversality of their intersection; moreover, 
7 is constant, by Corollary 15.111 □ 

Remark 5.13. The proof of Proposition 12.141 may be recognized as a 
generalization of the basic Evans function calculation pioneered by Evans 
[E.4], relating behavior near the origin to geometry of the phase space of 
the traveling wave ODE and thus giving an explicit link between PDE and 
ODE dynamics. The corresponding one-dimensional result was established 
in [GZ]; for related calculations, see, e.g., [J, AG J, PW]. 

5.3. Spectral bounds and necessary conditions for stability. 

Using Proposition 12.141 we readily obtain the stated necessary conditions 
for stability, Theorem 12.181 Define the reduced Evans function as 



(5.39) 



A(£,A) := ]im p-'Dip^pX). 



By the results of the previous section, the limit A exists and is analytic, 
with 



(5.40) 



A = 7 A, 
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for shocks of pure type (indeed, such a limit exists for all types). Evidently, 
A(-, •) is homogeneous, degree I. 

Lemma 5.14 ([ZS]). Let A(0, 1) ^ 0. Then, near any root (£oi\>) of 
A(-, •), there exists a continuous branch A(£), homogeneous degree one, of 
solutions of 

(5.41) A(£,A(£)) = 

defined in a neighborhood V o/£o; with A(£o) = Ao- Likewise, there exists a 
continuous branch A*(£) of roots of 

(5.42) D(lK(i)) = 0, 

defined on a conical neighborhood V po := {£ = pfj : fj G V, < p < po}, 
Pq > sufficiently small, "tangent" to A(-) in the sense that 

(5-43) |A.(f)-A($|= (|f|) 
as ||| ^0, /or£e V 

Proof. Provided that A(£o, •) 0, the statement (|5.41() follows by Rouche's 
Theorem, since A(£o, •) are a continuous family of analytic functions. But, 
otherwise, restricting A to the positive real axis, we have by homogeneity 
that 

0= Urn A(£o,A)= lim A(f /A, 1) = A(0, 1), 

A^+oo A^+oo 

in contradiction with the hypothesis. Clearly we can further choose A(-) 
homogeneous degree one, by homogeneity of A. Similar considerations yield 
existence of a branch of roots A(£,p) of the family of analytic functions 

(5.44) g^(\):=p- l D(pip\), 

for p sufficiently small, since />° = A(f, •). Setting A*(£) := |£|A(£/|£|, |f|), 
we have 

O(|,A*(|)) = |||V 7lll ' lll (A <! ) = 0, 

as claimed. "Tangency," in the sense of (J5.43)) . follows by continuity of 
A at p = 0, the definition of A*, and the fact that A(£) = |£|A(£/|£|), by 
homogeneity of A. □ 



12 Here, and elsewhere, homogeneity is with respect to the positive reals, as in most 
cases should be clear from the context. Recall that A (and thus A) is only defined for 
real £, and 5RA > 0. 
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Proof of Theorem \2.18l Weak spectral stability is clearly necessary for vis- 
cous stability. For, unstable (L p ) spectrum of (necessarily point spec- 
trum, by TODO) corresponds to unstable (L p ) essential spectrum of the 
operator L for p < oo, by a standard limiting argument (see e.g. [He, Z.l]), 
and unstable point spectrum for p = oo. This precludes L p — > LP stabil- 
ity, by the generalized Hille-Yosida theorem, (Proposition lA.lTl Appendix 
A). Moreover, standard spectral continuity results [Ka, He, Z.2] yield that 
instability, if it occurs, occurs for a band of £ values, from which we may 
deduce by inverse Fourier transform the exponential instability of (|^.<S|) for 
test function initial data Uq € Cq°, with respect to any L p , 1 < p < oo. 

Thus, it is sufficient to establish that failure of weak refined dynamical 
stability implies failure of weak spectral stability, i.e., existence of a zero 
D(i,X) = for £ £ KA > 0. Failure of weak inviscid stability, or 

A(£, A) = for £ 6 R d_1 , KA > 0, implies immediately the existence of such 
a root, by tangency of the zero-sets of D and A at the origin, Lemma 15.141 
Therefore, it remains to consider the case that weak inviscid stability holds, 
but there exists a root -D(£, it) for £, r real, at which A is analytic, A^ ^ 0, 
and (3{£,,iT) < 0, where (5 is defined as in (|2.17j) . 

Recalling that D(p^,pX) vanishes to order i in p at p = 0, we find by 
L'Hopital's rule that 

(d/dp) e+1 D(pi,pX)\ p=0 , x=iT = (l/£\)(d/dp)g^(0) 

and 

(d/d\ )D(plpir)\ p=0 , x=iT = (l/£\)(8/dX)^(0), 
where g^ p {X) := p- e D(pi,pX) as in (|5.44j) . whence 

(d/dp)g^(0) 



(5.45) (3 



(d/dx ) 9 i>\o) 



for defined as in Definition l2~H)l with (d/dX )/' A (0) ^ 0. 

By the (analytic) Implicit Function Theorem, therefore, A(£,p) is ana- 
lytic in £, p at p = 0, with 

(5.46) (d/dp) A(£,0) =-(3, 

where A(£, p) as in the proof of Lemma f5.14l is defined implicitly by g^ ' (p) = 
0, A(£,0) := it. We thus have, to first order, 

(5.47) X(tp) = ir-(3p + 0(p 2 ). 
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Recalling the definition A*(£) := |£|A(£/|£|,|£|), we have then, to second 
order, the series expansion 

(5.48) \*(p0 = i P T-pp 2 + O(p 3 ), 

where A*(£) is the root of -D(£,A) = defined in Lemma 15.141 It follows 
that there exist unstable roots of D for small p > unless > 0. □ 

5.4. Low-frequency resolvent estimates. It remains to establish the 
low-frequency resolvent bounds of Proposition 15.11 Accordingly, we restrict 
attention to arcs 

(5.49) I* : 5?A = 0(|f| 2 + |3A| 2 ), < |(f,3A)| < 5, 

with 9 > and 5 taken sufficiently small. 

5.4.1. Second-order perturbation problem. We begin by deriving a 
second-order, viscous correction of the central matrix perturbation problem 
(|5.29j) underlying the inviscid stability analysis [K, Ma. 1-3, Me.l]. Introduce 
the curves 

(5.50) (i A)(p, Co, t ) := {pio, piT - 9 lP 2 ), 

where £o S and tq € M are restricted to the unit sphere S d : |£o| 2 + 

|ro| 2 = 1. Evidently, as (£,o,tq,p) range in the compact set S d x [0,8], (£, A) 
traces out the surface U?r^ of interest. 

Making as usual the Ansatz U =: e Ma:i v, and substituting A = ipTQ — 9\p 2 
into (|5.27j) . we obtain the characteristic equation 



(5.51) 



2 B£ + p{-A l ± + ip £ B£t Qj +i P Y, B lk £ 

m mi 



-(i P ]T a%, + P 2 Y, Bjk &M + (p iT ° ~ ^p 2 ) 1 ) 

&1 jk^l 



v = 0. 



Note that this agrees with Q5.27JI up to second order in p, hence the (first- 
order) matrix bifurcation analysis of Lemma 15 . 91 applies for any fixed 9. We 
focus on "slow" , or "inviscid" modes /i ~ p. 
Positing the Taylor expansion 



(5.52) 



H = + H 1 p + 
v = v° H 
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(or Puisieux expansion, in the case of a branch singularity) as before, and 
matching terms of order p in (|5.51|) . we obtain 

(5.53) " * E A '^i ~ iT I)v = 0, 
just as in (|5.29j) . or equivalently 

(5.54) [(A 1 )~ 1 (ir + iA & ) - a /]v = 0, 

with //i =: — Qo, which can be recognized as the equation occurring in the 
inviscid theory on the imaginary boundary A = itq. 

In the inviscid stability theory, solutions of (|5.53|) are sub categorized into 
"elliptic" modes, for which ao has nonzero real part, "hyperbolic" modes, 
for which ao is pure imaginary and locally analytic in (£,t), and "glancing" 
modes lying on the elliptic-hyperbolic boundary, for which ao is pure imag- 
inary with a branch singularity at (£o) r o)- Elliptic modes admit a straight- 
forward treatment, both in the inviscid and the viscous theory; however, 
hyperbolic and glancing modes require a more detailed matrix perturbation 
analysis. 

Accordingly, we now restrict to the case of a pure imaginary eigenvalue 
c^o =: i%0i] here, we must consider quadratic order terms in p, and the 
viscous and inviscid theory part ways. Using p = —ipS,0i + °{p)i we obtain 
at second order the modified equation: 

(5.55) [(A 1 )- 1 ^ + p(B^° - 0{) + iA^°) - fii]v = 0, 

where d is the next order correction to a ~ —p)p, and v the next order 
correction to v. (Note that this derivation remains valid near branch singu- 
larities, since we have only assumed continuity of p/p and not analyticity 
at p = 0). Here, as usual denotes £o 3 -£o fc ) where £o := (£oi>£o)- 

Equation (|5,55|) generalizes the perturbation equation 

(5.56) [(A 1 )- 1 ^ + 7 + iA^ ) - al]v = 0, 

7 := 5RA — > + , arising in the inviscid theory near the imaginary boundary 
A = it [K, Me.4]. 

Note that tq is an eigenvalue of A^°, as can be seen by substituting 
= i£oi i n <1'5.54|) . hence |tq | < C|£o| and therefore (since clearly also 
l£o|>|fo|) 

(5.57) |eo|>(l/C)(|(|o,r )| = l/C. 
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Thus, for B positive definite, and 9\ sufficiently small, perturbation p{B*°*° — 
0\ ) , roughly speaking, enters (|5,55|) with the same sign as does 7/ in (|5.56|) , 
and the same holds true for semidefinite B under the genuine coupling condi- 
tion (|1.20j) : see (Kl), Lemma l3.18l Indeed, for identity viscosity B^ := 8^1, 
(|5.55j) reduces for fixed (£o> r ) exactly to (|5.56|) for Q\ sufficiently small, 
by the rescaling p — > /V(|£o| 2 — #1). This motivates the improved viscous 
resolvent bounds of Proposition I5.lt see Remark 15.51 

5.4.2. Matrix bifurcation analysis. The matrix bifurcation analysis 
for 1)5.55(1 goes much as in the inviscid case, but with additional techni- 
cal complications due to the presence of the additional parameter p. The 
structural hypothesis (H5), however, allows us to reduce the calculations 
somewhat, at the same time imposing additional structure to be used in the 
final estimates of Section 5.4.3. 

For hyperbolic modes, (£oj r o) bounded away from the set of branch sin- 
gularities U(£,77j(£)), we may treat (|5.55(l as a continuous family of single- 
variable matrix perturbation problem in p, indexed by (£o> r o); the resulting 
continuous family of analytic perturbation series will then yield uniform 
bounds by compactness. For glancing modes, (£o> r o) near a branch singu- 
larity, on the other hand, we must vary both p and (£,o,To), in general a 
complicated multi-variable perturbation problem. Using homogeneity, how- 
ever, and the uniform structure assumed in (H5), this can be reduced to a 
two-variable perturbation problem that again yields uniform bounds. For, 
noting that r)j(£) = 0, we find that £0 must be bounded away from the origin 
at branch singularities; thus, we may treat the direction £o/|£o| as a fixed 
parameter and vary only p and the ratio |to|/|£o|- Alternatively, relaxing 
the restriction of (£o>7"o) to the unit sphere, we may fix £0 and vary p and 
ro, obtaining after some rearrangement the rescaled equation 

(5.58) [(A 1 )' 1 (iro + [ia + p(B^° - 0i(|f o | 2 + |r + a\ 2 )] +iA^) - al]v = 0, 
where a denotes variation in To. 

Proposition 5.15 ([Z.3— 4]). Under the hypotheses of Theorem \2.1!A let 
a o = i£oi be a pure imaginary root of the inviscid equation (|5.54j) for some 
given £0, tq, i.e. det(A^ +To) = 0. Then, associated with the corresponding 
root a in (|5.55[) . we have the following behavior, for some fixed e, 9 > 
independent o/(£o, tq): 

(i) For (£0, tq) bounded distance e away from any branch singularity 
(£,?7j(£)) involving a, r/j as defined in (H5), the root a(p) in ([5.55(1 such 
that a(0) := ao bifurcates smoothly into m roots a±, . . . , a m , where m is the 
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dimension ofkei(A^° +tq), satisfying 



(5.59) 



or $taj < 



moreover, there is an analytic choice of eigenvectors spanning the associated 
group eigenspace V, in which coordinates the restriction Ay of the limiting 
coefficient matrices for Q5.17|) satisfy 



(5.60) 



< KA± < Cp 2 or - Cp 2 < 5?A± < 



C > 0, in accordance with (|5.59j) . for < p < e. 

(ii) For (£o> To) ^2/^5 a ^ a branch singularity (£,%'(£)) involving a, the 
root a(p,a) in (|5.58|) smc/i i/ia£ a(0, 0) = ao bifurcates (nonsmoothly) into 
m groups of s roots each: 



(5.61) 



«i, ■ ■ • • • , |«i , • • • ,a s ), 



where m is the dimension of ker(A^° + tq) and s is some positive integer, 
such that, for < p < e and |o~| < e, 



(5.62) 



«i =«+^+o(|o-i + ipd i/s , 



and, in appropriately chosen (analytically varying) coordinate system, the 
associated group eigenspaces V J are spanned by 



( 



(5.63) 



\ 







e-k 

(W) 2 e k 




V o / 



+ o(\a\ + \p\) 



l/s 



where 
(5.64) 



7T 



J ._ 



e^pa - iq k p) 1/s , 
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e := \ l l s , for p(£o) real-valued and uniformly bounded both above and away 
from zero and denoting the eigenvalues of Q € <C mxm such that sgnpQ > 
9>0, 

(5.65) IF := eH(apI m - iQp) 1 ' 8 

(well-defined, by definiteness ofp, Q), ek denote the standard basis elements 
in C m , and, moreover, the restrictions of the limiting coefficient matri- 
ces for (|5.17|) to the invariant subspaces V± are of form IF + o(|<r| + I/?!) 1 / 5 , 
satisfying 

(5.66) e P m j < < c P m j , or - c P m j < KA± < -e P m j , 

■= i( a - i P y/ s , C >0, /or < p < e. 

Proof. We here carry out the proof in the much simpler strictly hyperbolic 
case, which permits a direct and relatively straightforward treatment. A 
proof of the general case is given in Appendix C. In this case, the dimension 
of ker(^° + ro) is one, hence m is simply one. Let /(£o> r o) an d K^Oi^o) 
denote left and right zero eigenvectors of (A^° +tq) = 1, spanning co- kernel 
and kernel, respectively; these are necessarily real, since (A^° + ro) is real. 
Clearly r is also a right (null) eigenvector of (A 1 )~ 1 (iTQ + iA*°), and I A 1 a 
left eigenvector. 

Branch singularities are signalled by the relation 

(5.67) \A x r = 0, 

which indicates the presence of a single Jordan chain of generalized eigen- 
vectors of (A 1 )~ l (iTQ + iA^°) extending up from the genuine eigenvector r; 
we denote the length of this chain by s. 

Observation 5.16 ([Z.3]). Bound (|3JtO|) implies that 

(5.68) lB^°r > 6 > 0, 
uniformly in ^. 

Proof of Observation. In our present notation, (|«-i.5()j) can be written as 

(5.69) $ta(-iAt° - P B^°) < -8 lP , 

for all p > 0, some 9\ > 0. (Recall: |£o| > 62 > 0, by previous discussion). 
By standard matrix perturbation theory [Kat], the simple eigenvalue 7 = iro 
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of — iA^° perturbs analytically as p is varied around p = 0, with perturbation 
series 

(5.70) 7(p) = ir - plB^°r + o(p). 
Thus, 

(5.71) 3ft 7 (p) = -pW^°r + o(p) < -0 lP , 

yielding the result. □ 

In case (i), d(0) = a is a simple eigenvalue of (^4 1 ) _1 (iro + iA^q) : and so 
perturbs analytically in (|5.55j) as p is varied around zero, with perturbation 
series 

(5.72) a(p) = a + pp 1 + o(p), 

where p 1 = l(A 1 )~ 1 r, /, r denoting left and right eigenvectors of (j4 1 ) -1 (zto+ 
A^q). Observing by direct calculation that f = r, I = I A 1 /lA 1 r, we find that 



(5.73) p 1 = W^°r/lA 1 r 

is real and bounded uniformly away from zero, by Observation ^. lffl yielding 
the result (|5.59|) for any fixed (£,o, T o)i on some interval < p < e, where 
e depends only on a lower bound for p 1 and the maximum of ^"(p) on the 
interval < p < e. By compactness, we can therefore make a uniform choice 
of e for which (|5.59|) is valid on the entire set of (£o> to) under consideration. 
As Ay are scalar for m = 1, (|5.6Uj) is in this case identical to (|5.59() . 

In case (ii), a(0,0) = a is an s-fold eigenvalue of (A 1 )~ 1 (iro + iA£o), 
corresponding to a single s x s Jordan block. By standard matrix perturba- 
tion theory, the corresponding s-dimensional invariant subspace (or "total 
eigenspace") varies analytically with p and <r, and admits an analytic choice 
of basis with arbitrary initialization at p, a = [Kat]. Thus, by restricting 
attention to this subspace we can reduce to an s-dimensional perturbation 
problem; moreover, up to linear order in p, a, the perturbation may be 
calculated with respect to the fixed, initial coordinization at p, a = 0. 

Choosing the initial basis as a real, Jordan chain reducing the restriction 
(to the subspace of interest) of {A 1 )^ 1 {itq + iA£o) to i times a standard 
Jordan block, we thus reduce (|5.58|) to the canonical problem 

(5.74) (iJ + iaM + pN- (d - a))v 7 = 0, 
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where 



(5.75) 



J := 



/O 1 ••• 
OOIO 
1 

\o ••• 



0/ 



v/ is the coordinate representation of v in the s-dimensional total eigenspace, 
and M and N are given by 



(5.76) 

and 

(5.77) 



M := L(A 1 )~ 1 R 



N 



L(A l )- l (B^° -9i)R, 



respectively, where R and L are the initializing (right) basis, and its corre- 
sponding (left) dual. 

Now, we have only to recall that, as may be readily seen by the defining 
relation 



(5.78) 



L(A 1 )- 1 (iT + iAZ )R = J, 



or equivalently (A 1 y 1 (iT +iAi )R = RJ and L(A 1 )- 1 (iT +iAi ) = JL, the 
first column of R and the last row of L are genuine left and right eigenvectors 
f and I of (^4 1 ) _1 (iro + iA^q), hence without loss of generality 



(5.79) 



1= pi A 1 



as in the previous (simple eigenvalue) case, where p is an appropriate nonzero 
real constant. Applying again Observation l5.1fa1 we thus find that the crucial 
s, 1 entries of the perturbations M, N, namely p and pl(B^° — 0\)r =: q, 
respectively, are real, nonzero and of the same sign. Recalling, by standard 
matrix perturbation theory, that this entry when nonzero is the only sig- 
nificant one, we have reduced finally (modulo o(\a\ + \p\)^^ s ) errors) to the 
computation of the eigenvalues / eigenvectors of 



(5.80) 









1 
1 




\pa — iqp 



0/ 
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from which results (|5.62|) — (|5.64|) follow by an elementary calculation, for 
any fixed (£o,to), and some choice of e > 0; as in the previous case, the 
corresponding global results then follow by compactness. Finally, bound 
(|5.66j) follows from ()5.62[) and (|5.64j) by direct calculation. (Note that the 
addition of further 0(\a\ + \p\) perturbation terms in entries other than the 
lower lefthand corner of (|5.8U|) does not affect the result. Note also that 
A V j are scalar in the strictly hyperbolic case m = 1, and V± are simply 
eigenvectors of A±.) This completes the proof in the strictly hyperbolic 
case. □ 

Remark 5.17. The detailed description of hyperbolic and glancing modes 
given in Proposition l5.15l readilv yield the result of Lemma f5.9l in the deferred 
case KAo = 0, under the additional hypothesis (H5) (exercise). 

5.4.3. Main estimates. Combining the Evans-function estimates of 
Proposition 12.141 (and, in the case that refined but not uniform dynamical 
stability holds, also those established in the course of the proof of Propo- 
sition I2.18|) with the matrix perturbation analysis of Proposition I5,15( we 
have all of the ingredients needed to carry out the basic L 1 —* I? resolvent 
estimates of Proposition 15.11 In particular, in the uniformly dynamically 
stable case, they may be obtained quite efficiently by Kreiss symmetrizer 
estimates generalizing those of the inviscid theory. We refer to [GMWZ.l] 
or the article of M. Williams [W] in this volume for a presentation of the ar- 
gument in the strictly hyperbolic, Laplacian viscosity case. With the results 
of Proposition 15.151 a block version of the same argument applies in the 
general case, substituting invariant subspaces V J for individual eigenvectors 
V. We omit the details, which are beyond the scope of this article. 

In the refined, but not uniformly dynamically stable case, the estimates 
may be obtained instead as in [Z.3-4] by detailed pointwise estimates on the 
resolvent kernel, using the second-order Evans function estimates carried out 
in the course of the proof of Proposition 12. 18l and the explicit representation 
formula for the resolvent kernel of an ordinary differential operator [MZ.3, 
Z.4]. We refer to [Z.3-4] for an account of these more complicated arguments. 

5.4.4. Derivative estimates. Improved derivative estimates, = 1, 
may now easily be obtained by a method introduced by Kreiss and Kreiss 
[KK] in the one-dimensional case; see [GMWZ.l] or the notes of M. Williams 
[W] in this volume. Specifically, given a differentiated source / = d xi g, in 
resolvent equation {L^ — X)U = /, consider first the auxiliary equation 

(5.81) L W = (B U W X1 ) X1 - {A l W) Xl = 8 xi g. 
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Using the conservative (i.e., divergence- form) structure of Lo, we may 
integrate (|5.81f) from — oo to x% to obtain a reduced ODE 

(5.82) z' 2 - a{x x )z 2 = y G , W = *(z 2 ) 

analogous to that obtained in [KK, GMWZ.l] for the strictly parabolic case, 
where z 2 ■= BjjW and G ~ g, W ~ z 2 . More precisely, the inhomogeneous 
version z' 2 — a{x\)z = is the linearization about u of the (integrated) 
traveling- wave ODE (|3.5j) . from which we obtain by transversality 7 / 
that the solution of (|5.82|) is unique modulo v! . That a solution exists fol- 
lows easily from the fact (|3.7|) relating the signatures of a± to those of A\. 
(a consequence of transversality, together with our assumptions on the pro- 
file; see Remark 13.5(1 . together with standard arguments for asymptotically 
constant ODE as in, e.g., [He, Co, CL]; see, in particular, the argument of 
section 10.2, [GMWZ.l] in the strictly parabolic case, for which a reduces to 
A . See also [Go. 2, KK], or the article [W] by M. Williams in this volume. 
Indeed, imposing an additional condition 

(5.83) (£, W) = 0, 

where £ is any constant vector satisfying (£, v!) = £ ■ [U] 7^ 0, we have [He, 
Co, CL, GMWZ.l, Go.2, KK, W] the bound \z 2 \ w i, P < C\G\ LP for any p, 
yielding in particular 

(5.84) \W\ L i + \B n W Xl \ L i <C\g\ L i. 

The reduction to form (|5.82|) goes similarly as the reduction of the nonlinear 
traveling- wave ODE in Section 3.1; we leave this as an exercise. 

Setting now U = W + Y, and substituting the auxiliary equation into 
the eigenvalue equation, we obtain equation 

{Li ~ A)y = 0(p)(\W\ L i + \B u W Xl \ L i) 

(5.85 ) 

= 0(p\g\ L i), 

for the residual Y, from which we obtain the desired bound from the basic 
L 1 — > I? estimate of Section 5.4.3 above. 

Alternatively, one may obtain the same bounds as in [Z.3-4] by direct 
computation on the original resolvent equation. Both methods are based 
ultimately on the fact that, for Lax-type shocks, the only L 1 time-invariants 
of solutions of the linearized equations are those determined by conservation 
of mass. This property is shared by over- but not undercompressive shocks, 
hence the degraded bounds in the latter case; for further discussion, see 
[LZ.2, Z.2]. This completes the proof of Proposition 15. 1| and the analysis. 
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A Appendix A. Semigroup facts 



Definition A.l. Given a Banach space X, and a linear operator L : T>(L) C 
X — > X, we say that L is densely defined in X if is dense in X. We 

say that L is closed if u n —> u and x n := Lu n — > x (with respect to | ■ |x) for 
u n 6 and x n 6 X implies that u £ D(L) and Lu = x. T>{L) is called 

the domain of L. We define associated domains V(L n ) by induction as the 
set of x € V(L n ~ l ) such that L^x G V(L). 

Exercises A. 2. 1. If L is closed, show that T>(L) is a Banach space under 
the canonical norm |«|x>(L) := \u\x + i-e., each Cauchy sequences 

u n € T>(L) has a limit u £ T>(L). (First note that u n and x n := Lu n are 
Cauchy with respect to \ ■ \x, hence have limits u, x in X.) With this choice 
of norm, L : T>(L) — > X is a bounded operator, hence (trivially) closed in 
the usual, Banach space sense. 

2. Show that L — A is closed if and only if is L. 

3. Show that \u\ < C\Lu\ for L closed implies that range(L) is closed. 

Definition A. 3. Given a Banach space X, and a closed, densely defined 
linear operator L : T>(L) C X — > X, the resolvent set p(L), written px(L) 
when we wish to identify the space, is defined as the set of A € C for which 
(A— L) has a bounded inverse (A — L)~ l : X — ► V(L). The operator (A — L)~ l 
is called the resolvent of L. The spectrum a(L) of L, written crx(L) when 
we wish to identify the space, is defined as the complement of the resolvent 
set, a(L) := p{L) c . 

Exercise A. 4. If L is densely defined and A £ p(L) for some X, show that 
V(L n ) = range(A — L)~ n is dense in X. (Show by induction that each 
V{L n+1 ) is dense inV(L n ).) 

Lemma A. 5. For a closed, densely defined operator L on Banach space X, 
the resolvent operator (A — L)~ l is analytic in A with respect to \ ■ \x for A 
in the resolvent set p(L). 

Proof. For A sufficiently near Ao € p(L), we may expand 

A)(Ao-L)- 1 )" 1 

A)(A -L)- 1 ) J , 

using | (Ao — A) (A — L)~ 1 \ < C\Xo — A| by the definition of resolvent set and 
the Neumann expansion (I — T) _1 = Yl'jLo'T-' for \T\ sufficiently small. □ 



(A - L) = (A - L)~ [I — (X 



(A.l) 



oo 

(Ao - L)~ 1 £ ((Ao 
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Exercise A. 6. 1. Let L(a) be a family of closed, densely defined operators 
on a single Banach space X, such that L is "relatively analytic" in a with 
respect to (L(ao) — Ao), Ao £ p{L{ao), in the sense that L(a)(Ao — L(ao)) 1 
is analytic in a with respect to \ ■ \x- Show that the family of resolvent 
operators (A — L(a))^ 1 is analytic in (a, A) in a neighborhood of (ao, Ao). 

2. Including in (|4.35j) formerly discarded beneficial w terms in the 
energy estimates from which it derives, we obtain the sharpened version 

(A.2) ^ X + B 1 )( y \W\ fll + \d 2 Xl w II \) <Ci(|/| Al + Ci|W|), 

revealing smoothing in w 11 for 3?A > —B\. Show that this implies analyticity 
in £ with respect to \ • \ of L^(\q — L^)" 1 (note: dgLg is a continuous- 
coefficient first-order differential operator, for which the derivative falls only 
on w 11 components), and conclude that (A — L^) _1 is analytic in (£, A) for 
A £ pjj 1 (L^) and 3?A > —B\ with B\ > sufficiently small. 

Lemma A. 7 (Resolvent identities). For X, [i in the resolvent set of a 
closed, densely defined operator L : T>(L) — > X on a Banach space X, 

(A.3) (A - L)- V " L)-' = (M " L) "'" (A " L) " 1 = 0* - L)-\X - L)- 1 

A — /x 

on X and L(A — L)^ 1 = A(A — L) _1 — 1= (A — L)~ 1 L onV(L): in particular, 



(A.4) (A-L)- X n = A~ 1 (u + (A-L)- 1 Lu) for u € T>{L). 

Proof. Rearranging (u — L)(u — L)~ l = I, we obtain (A — L)(/i — L) _1 = 
(A — jti)(/i — L)" 1 /, from which the first equality of (|A.3|) follows upon appli- 
cation of (A — L) _1 from the left, and the second by symmetry. Rearranging 
defining relation (A — L)(X — L)~ l = I = (A — L) _1 (A — L), we obtain the 
second assertion, whereupon ()A.4jl follows by further rearrangement after 
multiplication by A -1 . □ 

Exercise A. 8. Assuming that the resolvent set is open, recover the result of 
Lemma lA.fl directly from the definition of derivative, using resolvent identity 
(|A.3|) to establish differentiability, (d/dX)(X — L)^ 1 = —(A — L)~ 2 . 

Definition A. 9 ([Pa]). A C° semigroup on Banach space X is a family of 
bounded operators T(t) satisfying the properties (i) T(0) = /, (ii) T(t + s) = 
T(t)T(s) for every t, s > 0, and (iii) lim t _+ + T(t)x = x for all x € X. The 
generator L of the semigroup is defined as Lx = lim t _ >0 + (T(t)x — x)/t on the 
domain T>(L) C X for which the limit exists. We write T(t) = e Lt . Every 
C semigroup satisfies \e Lt \ < Ce 70 * for some 70, C; see [Pa], Theorem 2.2. 
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Remark A.10. For a C° semigroup, (d/dt)e Lt f = Le Lt f = e Lt Lf for all 
/ £ T>{L) and t > 0; see [Pa], Theorem 2.4(c). Thus, e Lt is the solution op- 
erator for initial-value problem ut = Lu, u(0) = /, justifying the exponential 
notation. 

Proposition A. 11 (Generalized Hille Yosida theorem). An operator 
L : T>(L) — > X is the generator of a C° semigroup \e Lt \ < Ce 70 * on X with 
domain T>(L) if and only if: (i) it is closed and densely defined, and (ii) 
A € p(L) and |(A — L)~ k \ < C\X — 7o|~ fc for sufficiently large real X, in which 
case also |(A - L)~ k \ < C|5RA - 7 |~ fc for all 3?A > 70. 

Proof. (=>) By the assumed exponential decay, the Laplace transform T := 
J °° e- Xs T(s) ds is well-defined for KA > 7o , with |T(A)| < C|KA - t- | _1 - 
By the properties of the solution operator, we have also TL = LT = LT on 
V(L), as well as LT = d t T = XT — I. Thus, (A - L)f = f (A - L) = I on 
V(L). But, also, 



(A.5) h-\T{h)-I)f = h-\e Xh -I)T -h- 1 , 

Jo 



-X(s-h) 



T(s)ds 



approaches AT - I as h -> 0+, yielding f X C V(L) and LT = XT — I 
on X, by definition. Thus, A € p(£) for 3ftA > 70, and T = (A — L) _1 , 
whence |(A — L)~ k < C|KA — 7o|~ fc for k = 1. The bound for general 

fc times 



follows from the computation f k = £ (T * T ■ ■ ■ * T) = £(t k T/k\), where £ 
denotes Laplace transform and "*" convolution, g*h(t) := g(t — s)h(s) ds, 
together with f °° e~ z z k /kl dz = 1 for all > (exercise). 

(■<=) Without loss of generality, take 70 0. For real numbers z, e z = 
{e^ z y 1 = lim„_ >0O (l — z/n)~ n gives a stable approximation for z < 0. This 
motivates the introduction of approximants 

(A.6) T n := (/ - Lt/n)- n , 

which, by resolvent bound (ii), satisfies \T n \ < C. Restricting to the dyadic 
approximants T 2 j, and using the elementary difference formula a n — b n = 
(a n_1 + a n ~ 2 b + • • • + b n ^ 1 )(a — b) for commuting operators a, b together 
with bound (ii), we find for x G T>(L 2 ) that 
(A.7) 

|(T 2J+ i - T 2j )x\ < VC 2 \ ((/ - Lt/y +1 )~ 2 -(I- Lt/V)- 1 ^ 



x 



2 3 C i 



(I - Lt/2 j+1 )- 2 {I - Lt/2 j 



t 2 \L 2 x\2- 2j - 2 



< CH 2 \L 2 x\2- j - 
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and thus the sequence converges geometrically to a limit Tx for all x € 
T>(L 2 ). Since T>(L 2 ) is dense in X fExercise IA.4(I . and the T n are uniformly 
bounded, this implies convergence for all x, defining a bounded operator 
T(t). 

Clearly T(0) = T n (0) = /. Also, T n x -> x as t -> for x £ V(L n ), 
whence we obtain Tx — > x for all a; by uniform convergence of the dyadic 
approximants and density of convergence of the approximants and density 
of each V(L n ) in X. Finally, for x E V(L), we have {d/dt)(I - Lt/n)x = 
—Lx/n, from which we obtain (d/dt)T n x = —n(I — Lt/n)~ n ~ 1 (—L/n)x = 
L n (t)T n x, where L n (t) := L(I—Lt/n)~ l is bounded for bounded t (exercise). 
Observing that \LT n (t)x\ = \T n Lx\ < C|Lx| for all t, we find for x E T>{L 2 ) 
that 

(A.8) e n (r, t) := u{t) - v(t) = T n {r + t)x - T n (t)T n (r)x, 

v! = L n (s)u, v' = L n (r + s)v, u(0) = v(0) := T n (r)x, satisfies 
(A.9) e' n = L n (s + r)e n + f n (r,s), e(0) = 0, 

f n (r, s) := (L n (s) - L n (s + r))T n {r + s)x, where 

\f n (r,s)\ = {((I-Ls/ny 1 -(I-L(s + r)/n)-^T n Lx\ 

( A -!0) =\(I ~ Ls/n)- 1 ^ - L{s + r)/n)- l T n {r /n)L 2 x\ 

< |L 2 x|r/rz 

goes to zero uniformly in n for bounded r. Expressing 

(A.ll) e n (r,t)= [ T n (t + r)T n {t + r-s)- 1 f n (r,s)ds 

Jo 

using Duhamel's principle/uniqueness of solutions of bounded-coefficient 
equations, and observing ( Exercise IA. 131 below) that \T n (t + r)T n (t + r — 
s )~ 1 \ < C, we thus have |e n (r, t)\ < Ct|L 2 x|r/n — > uniformly in n for 
bounded r, t, verifying semigroup property (ii) Definition IA.9I in the limit 
for x G T>(L 2 ), and thus, by continuity, for all □ 

Remark A. 12 ([Pa]). The approximants T n {t) correspond to the finite- 
difference approximation obtained by first-order implicit Euler's method 
with mesh t/n. In the language of numerical analysis, boundedness of \T n \ 
corresponds to "^4-stability" of the scheme, i.e., suitability for "stiff" ODE. 
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Exercise A. 13. Show by direct computation that 

(A.12) T n {r + ,)T n ( S )- 1 = ((-£-) (J - L(r + ^/n)" 1 + (^))" 

Assuming resolvent bound (ii), and thus \(I — L(r + s)/n)~ k \ < C, show 
that \T n (r + s)T n (s)' 1 | < C, using binomial expansion and the fact that 
resolvents commute. 

Exercise A. 14. Show by careful expansion of 

T n+1 -T n = ((I - Lt/{n + l))-(«+!) - (/ - Lt/(n + 1))"»>) 

(A. 13) 

+ ((/ - Lt/(n + l))" n - (/ - Lt/n))- n ) 

that \T n+ \x — T n x\ < C(\L 2 x\ + \Lx\ + |x|)/n 2 , so that the entire sequence 
{T n } is convergent for x G V(L 2 ). 

Exercise A. 15. Show directly, using the same Neumann expansion argu- 
ment used to prove analyticity of the resolvent that |(A — L)~ 1 \ < (A — 70)" 1 
for sufficiently large real A implies |(A — L)~ l \ < (JRX — jq)" 1 for all 3?A > 70. 

Definition A. 16. We say that a linear operator L is dissipative if it satisfies 
an a priori estimate 

(A.14) \\-lo\ k \U\x<C\(\-L) k U\ x 

for all real A > Ao- Condition (|A,14|) together with range(A — L) = X is 
equivalent to condition (ii) of Proposition IA.11I 

Corollary A. 17. An operator L : T>(L) — > X is the generator of a C° 
semigroup \e Lt \ < Ce 70 * on X with domain T>(L) if and only if: (i) it is 
closed and densely defined, and (ii 7 ) both L and L* are dissipative, where 
L* : X* — > V(L)* defined by {L*v,u} := {v,Lu) denotes the adjoint of L. 

Proof. See Corollary 4.4 of [Pa], or exercise lA~20l below. □ 

Corollary A. 18 (Generalized Lumer Phillips theorem). A densely 
defined operator L : T>(L) — > X generates a C° semigroup \e Lt \ < Ce 70 * on 
X with domain T>(L) if and only if L is dissipative and range(Ao — L) = X for 
some real Ao > 70 • In particular, for a densely defined dissipative operator L 
with constants C , 70, the real ray (70, +00) consists either entirely of spectra, 
or entirely of resolvent points. 
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Proof. See Theorem 4.3 of [Pa], or Exercise IA.21I below. 



□ 



Remark A. 19. In the contractive case, C = 1, dissipativity is equivalent 
to $t(v,u) < for some v € X* such that {v,u) = \u\\v\; see Theorem 4.2, 
[Pa]. If X is a Hilbert space, contractive dissipativity (C = 1) of both L and 
its adjoint L* reduce to the single condition 3l(u, Lu) < 0, which is therefore 
necessary and sufficient that L generate a C° contraction semigroup; see 
Exercise IA.22I This corresponds to (d/dt)(l/2)\u\ 2 < 0, motivating our 
terminology. 

Exercise A. 20. For a closed operator L, show that (ii) of Proposition VA . 1 1\ 
is equivalent to ( ii ') of Proposition \A.17\ The forward direction follows by 
the general facts that \A\ = \A*\ and (A" 1 )* = (A*)^ 1 . The reverse follows 
by closure of range(A — L) (Exercise \A.2l 3). the Hahn-Banach Theorem, 
and the observation that = (/, (A — L)u) = ((A — L)*f, u) for all u G T>{L) 
implies / = 0, which together yield range(A — L) = X . 

Exercise A.21. Let L be a dissipative operator with constants C, 70, such 
that range(Ao — L) = X for some real Ao > 70 • 

1. Using Exercise \A . 3L S. show that L is closed. 

2. If (Ai — L)u\ = (A2 — L)u2 for Xj > 70, show that \u± — U2I < 
C|Ai - 7or 1 |Ai - A 2 1 |ii 2 1 . 

3. Using the result of 2, show that range(A n — L) = X for X n — > A > 70 
implies range(A — L) = X. (Show that {uj} is Cauchy for (Xj — L)uj := x, 
then use closure of L.) Conclude that range(A — L) = X for all real X > 70, 
by the fact that p(L) is open. 

Exercise A. 22. For a densely defined linear operator L on a Hilbert space 
X, and u € T>(L), show by direct inner-product expansion that \(X-L)u\ 2 > 
I (A — 7o)^| 2 f or rea l A > 70 if and only if (u,Lu) < 7o|u| 2 , if and only if 
(u, L*u) < 7o['u| 2 . 

Remark A. 23. Proposition IA. Ill includes the rather deep stability esti- 
mate \e Lt \ < Ce 70 * converting global spectral information to a sharp rate of 
linearized time-exponential decay. For a Hilbert space, a useful alternative 
criterion for exponential decay \e LT \ < Ce 70 * has been given by Priiss [Pr]: 
namely, that, for some 7 < 70, {A : 5ftA > 7} C p(L), and [(A — L)^ 1 < M 
on 3?A = 7. A useful observation of Kapitula and Sandstede [KS, ProK] is 
that, in the case that there exist isolated, finite- multiplicity eigenvalues of 
L in {A : IRA > 7} C p(L), the same result may be used to show exponential 
convergence to the union of their associated eigenspaces. 
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However, these tools are not available in the here, that essential 

spectrum of L approaches the contour 5?A = 70 (70 = in our case). In this 
situation, one may take the alternative approach of direct estimation using 
the inverse Laplace transform formula, given just below. However, notice 
that we do not get from this "local" formula the Hille-Yosida bound, which 
is not implied by behavior on any single contour 5RA = 7 (recall that A — ► 00 
in the proof of Proposition lA.lTj) . Indeed, the spectral resolution formula 
holds under much weaker bounds than required for existence of a semigroup; 
see Exercise IA.26I A theme of this article is that we can nonetheless get 
somewhat weaker, time-averaged bounds by similarly simple criteria, and 
that, provided we have available a nonlinear smoothing estimate analogous 
to (|3.64j) . we can use these to close a nonlinear argument in which the 
deficiencies of our linearized bounds disappear. The result is effectively a 
"Priiss-type" bound on the high-frequency part of the solution operator, 
analogous to the bounds obtained by Kapitula and Sandstede by spectral 
decomposition in the case that slow- and fast-decaying modes are spectrally 
separated. 

Proposition A. 24. Let L : T>(L) — > X be a closed, densely defined operator 
on Banach space X, generating a C° semigroup \e Lt \ < Ce 70 *: equivalently, 
satisfying resolvent bound ()3.16j) . Then, for f £ T){V), 

/■7+ioo 

(A.15) e Lt f = P.X. e xt (X- L)' 1 f dX 

J 7— ioo 

for any 7 > 70, with convergence in L 2 ([0, 00); X). For f, Lf £ T>{V), 
()A.15|) converges pointwise for t > 0, with uniform convergence on compact 
sub-intervals t € [e, 1/e] at a rate depending only on the bound for \ f\x + 
\LF\ x + \L{Lf)\ x . 

Proof. By the bound \e Lt \ < Ce 70 *, the Laplace transform u of u(t) := e Lt f 
is well-defined for 9?A > 70, with 

(A.16) u(t) = P.V. / e xt u(X)dX 

J 7— ioo 

for 7 > 70, by convergence of Fourier integrals on L 2 , and the relation 
between Fourier and Laplace transform (see (14.581) ). Moreover, u(t) € 
C 1 ([0,oo);X) bv Remark lA.101 with Uj = Lu and \u t \, \Lu\ < Ce 70 *, whence 
ut = Xu — u(0) = Xu — f and also Lu = Lu. Thus, (A — L)u = f, giving 
(|A.15|) by (|A.16|) and invertibility of A - L on KA > 70. 
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Next, suppose that /, Lf G 'D(L), without loss of generality taking 
A ^ 0. 13 Expanding 

(A.17) (A - L)- 1 / = A- 2 ((A - L)- l L ■ Lf + Lf) + A" 1 / 

using (|A.4|1 . we may split the righthand side of <|A. 15|> into the sum of an 
integral 

(A.18) P.V. / e xt X~ 2 ((A - L)-' L L ■ Lf + Lf) dX 

J -/—too S 

that is absolutely convergent on bounded time-intervals t G [0, <5 _1 ] and 
the integral P.V. e A *A _1 dXf = f, uniformly convergent on compact 

intervals t E [6, 5 ]; see Exercise IA.27I below. Note that (|A.18|) converges 
to 

pry _|_200 

(A.19) P.V. / A~ 2 ((A - LY l L ■ Lf + Lf) dX = 

as t — > 0, directly verifying semigroup property (iii). For an alternative 
proof, see Corollary 7.5 of [Pa]. □ 

Proposition A. 25. Let L : T>{L) — > X be a closed, densely defined operator 
on Banach space X, generating a C° semigroup \e Lt \ < Ce 70 *: equivalently, 
satisfying resolvent bound (|oM6j) . Then, for f G L 1 ([0, T]; or f G 

L 2 ([0,T]-X), 

(A.20) / e L ^-^f{s)ds = PX. e xt (X- L)- 1 f T (X)dX 

JO J ry— ioo 

for any 7 > 70, < t < T, with convergence in L 2 ([0, T];A), where g 
denotes Laplace transform of g and f T (x,s) := f(x,s) for < s < T and 
zero otherwise. For f, Lf G L 1 ([0, T];V(L)) or f G L q ([0, T];V(L)), q > 1, 
the convergence is pointwise, and uniform for t G [0, T] . 

Proof. Equivalence and convergence in L 2 ([0, T];X) for / G ^ 1 ([0, T]; T>(L)) 
follows similarly as in the proof of Proposition IA.151 observing that u : = 
jt e L(t-s)f^ ds satisfies Ut - Lu = f T for <t <T with u(0) = 0, and u, 

f are well-defined, with (A — L)u = f T . Likewise, we may obtain pointwise 

13 The restriction 7 / may be removed by considering the shifted semigroup T{t) := 
e~ 70 'e Li . (Exercise: verify that T is generated by L — 70.) 
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convergence for /, Lf € L 1 ([0, T]; P(L)) using expansion l)A,17|) and the 
Hausdorff- Young inequality 



(A.21) \f(X)\ = | F e~ Xs g T {s)ds\ < (1 + e~^ T )\g\ Lh , , 

J o 



to obtain a uniformly absolutely convergent integral plus the uniformly con- 
vergent integral 
(A.22) 

/■7+ioo pT />7+«oo 

P.V./ \- 1 e xt f T (\)d\= /(s)dsP.V. / \- 1 e x(t ~ s) f T (X)d\ 

J 7— ioo JO J -/—too 

ft 



f(s)ds; 



see Exercise I A. 271 below. 

Convergence in 1? for / € L 2 ([0,T]; X) follows using Parseval's identity, 

(A.23) 

P.V./ |^(A)| 2 dA= / |e-^( S )| 2 d S <(l + |e-^|) 2 b|i 2[0iT] , 

J ■y—ioo JO 

together with the fact that |e A *(A — L)~ 1 \ < C(^) for 5?A = 7. We may 
then conclude equivalence by a limiting argument, using continuity with 
respect to / of the left-hand side and equivalence in the previous case 
/ £ L 1 ([0, T]; I?(L)). A similar estimate together with expansion (|A.17|) 
yields pointwise convergence for / € L q ([0,T];V(L)), q > 1, using Holder's 
inequality \X~ l g T \i J 2/x) ^ C\\~[q\g T \LP and the Hausdorff- Young inequality 
\g T \LP < C\g\ L 2 [QtT] , where 1/p + 1/q = 1. □ 

Exercise A. 26. Let L : T>(L) ^> X be a closed, densely defined operator on 
Banach space X, satisfying resolvent bound (|3.16|) (equivalently, generating 
a C° semigroup \e Lt \ < Ce 7oi J. 

1. For / E F>(L), show by judicious deformation of the contour that 
P.V. J2_fo£ e xt (\ — L)~ 1 fd\ = for any 7 > 70 and t < 0, w/zi/i uniform 
pointwise convergence on compact intervals [— e, — e _1 ], e > ; at a rate 
depending only on the bound for \f\x + \LF\x- Indeed, this holds under the 
resolvent growth bound |(A — L)~ 1 \ < C(l + [A[) r for -RA > 7, /or any r < 1. 

2. For f G T>{L), show without reference to semigroup theory that the 
righthand side of (jA.15|) converges in L 2 (x,t) to a weak solution in the 
sense of [Sm] of initial value problem ut = Lu, u(0) = f, using Parseval's 
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identity and resolvent identity (|A.4I) together with the fact that distributional 
derivatives and limits commute. (By moving 7 to u) > 0, show that the 
L 2 (t;X) norm of e-^\e Li f - f) = J e^^A^A - L^LfdX is order 
(u — 7o)~ 1 , yielding a trace at t = 0.) Indeed, this holds under the resolvent 
growth bound |(A - L)- l \ < C(l + |A|) r for = 7, for any r < 1/2. 

Exercise A. 27. Show that P.V. f 7+4 °° z~ x e z dz is uniformly bounded, in- 

J7 — ZOO J * ' 

dependent of^y, and converges uniformly to sgn(7) for 7 bounded away from 
zero. 

Definition A. 28. A C° evolutionary system on Banach space X is a family 
of bounded operators U(s,t), s < t, satisfying the properties (i) U(s,s) = 
I, (ii) U(s,t + r) = U(s,t)U(t,t + r) for every s < t, r > 0, and (iii) 
lim t ^ s + U(s, t)x = x for all x G X. The instantaneous generators L(s) 
of the system are defined as L(s)x = lim t ^ s +(U(s,t)x — x)/(t — s) on the 
domain D(L(s)) C X for which the limit exists. 

Remark A. 29. For a C° evolutionary system, (d/dt)U(s, s)f = L(s)U(s, s)f 
for all / € V(L(s)); see [Pa], p. 129. Thus, if there is a common dense sub- 
space Y C T>(L(t)) for all t that it is invariant under the flow U(s,t), then, 
for / G Y, U(s, t) is a solution operator for initial-value problem ut = L(t)u, 
u(s) = f, generalizing the notion of semigroup to the nonautonomous case. 

Proposition A. 30. Let L(s) : T>(L(s)) — > X be a family of closed, densely 
defined operators on Banach space X , each generating a contraction semi- 
group je^^'l < e 7o( - s )* . Suppose also that there exists a Banach space 
Y C f] s 'D(L(s)) contained in the common intersection of their domains, 
dense in X , invariant under the flow of each semigroup e L ^ t , on which e^)* 
generates a semigroup also in the Y -norm, and for which L(s) : Y — > X is 
continuous in s with respect to the X —* Y operator norm. Then, L(s) 
generates a C° evolutionary system on X with domains T>(L(s)). 

Sketch of proof. Substituting for L(t) the piecewise constant approxima- 
tions L n (t) := L(k/n) for t 6 [k/n, (k + l)/n), we obtain a sequence of 
approximate evolutionary systems U n (s,t) satisfying the uniform bound 
\U(s,t)\ < e 7 °(*- s ) for s < t. For x G Y C V(L(s)) for all s, the sequence 
U n (s,t)x may be shown to be Cauchy by a Duhamel argument similar to 
the one used in the proof of Proposition lA.lll See Theorem 3.1 of [Pa] for 
details. □ 

Remark A. 31. There exist more general versions applying to systems for 
which the generators are not contractive, but stability in this case is difficult 
to verify. See [Pa], Section 3, for further discussion. 
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B Appendix B. Proof of Proposition 11.211 

Exercise B.l ([Kaw]). (Optional) 1. Using (|3.5U|) . show that 

(B.l) |e p -(^| <Ce- e ^ 2t l^\ 2 \ 

where P_(»0 := -i£^ i&A* " Ej.fcfc^*)-- 

2. Denoting by S(t) the solution operator of the constant- coefficient prob- 
lem 

(B.2) U t + Y,*Lu xj = J2(Bi k U Xk ) Xj , 

j j,k 

we have S(t)Uo = e p ~^ l UQ, where' denotes Fourier transform in x. De- 
composing 

(B.3) S = 5i + 5 2 , 

where Sj(t)Ua = Xj(O eP ~^ t ^o> where xi(0 7S one f or |£| < an ^ 2:67,0 
otherwise, and \2 = 1 — Xl> show, using Parseval's identity \f\L 2 ( x ) = I/Il 2 (£ 
and the Hausdorff-Young inequality \f\L°°(s,) < 

(B.4) |Sa(*)<9°/| L2 < C(l + t)- d / 4 -l-l/ 2 ]/Ux, [a| >0, 

and 

(B.5) |S 2 (0/Il* <Ce- e *|/| L2 . 

Proof of Proposition \1.21\ Taylor expanding about U = U- , we may rewrite 
(Ol) as 

(b.6) r/ t + £ ^ic/ 2 , = YSP-U*!.)*! = d Mu, d x U), 

j j,k 

where \Q(U,d x U\ < C\U\\d x U\ and \d x Q(U,d x U\ < C(\U\\d%U\ + \d x U\ 2 ) so 
long as \U\ remains uniformly bounded. Using Duhamel's principle/ variation 
of constants, we may thus express the solution of ()1.2|) as 

(B.7) U(t) = S{t)U + f S(t - s)d x Q(U, dU x ){s) ds, 

Jo 

where £(•) = Si(-) + Sk(-) is the solution operator discussed in Exercisc lB.il 
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Defining 

(B.8) ((t):= sup \U(s)\ Hs (l + s) d ^, 

0<s<t 



and using ()B.4|) and (|B.5|) . we may thus bound 
(B.9) 

\U(t)\ I ?<\Si(t)Uo\ I s + \S 2 (iL)Uo\ I ?+ [ \S 1 (t-s)d x Q(U,dU x )(s)\ L 2ds 

Jo 

+ f \S 2 (t-s)d x Q(U,dU x )(s)\ L 2ds 
Jo 

<C{l + t)- d / A \Uv\ L i+Ce- et \U Q \ L 2 

+ C [ (l + t- 8 )- d ^ 2 \Q(U,dU x )(s)\ L i ds 
Jo 

+ C f e- e ^\d x Q(U,dU x )(s)\ L 2ds 
Jo 

<C((l + t)- rf / 4 |C/(0)| LlnL2 

+ CC(t) 2 [ (1 + t - s)- d/4 ~ 1/2 (l + s)- d ' 2 ds 
Jo 

<c(i + t)^/ 4 (|c/(o)| LlnL2 + c(t) 2 ) 

so long as |C/|j?« remains uniformly bounded. Applying now (|3.75|) . we 
obtain (exercise) 

\U{t)\ 2 Hs < Ce- 9t \U(0)\ 2 Hs +C f e- e ^\U(s)\ 2 L 2 ds 

Jo 

(R10) < Ce- et \U(0)\ 2 Hs + C(l + t)- d /\\U(0)\»niP + at?? 

<C(l + t)- d / 2 (\U(0)\ LlnH s+at) 2 ) 2 , 

and thus 

(B.n) at)<c(\u(o)\ LlnH s+at) 2 ). 

From (|B.11|I . it follows by continuous induction (exercise) that 

(B.12) C(*)<2C|l7(0)| £ i nH . 

for \U(0)\ii n jj s sufficiently small, and thus 

(B.13) P(t)\ H > < 2C(1 + t)- d '*\U(0)\v nH . 
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as claimed. Applying we obtain the same bound for \U(t)\L°°, and 

thus for \U(t)\LP, all 2 <p< oo, by the L p interpolation formula 

(B-14) \f\ LP * < 

for all 1 < pi < p* < P2, where j3 := p\(j>2 — P*)/p*(P2 — Pi) is determined 
by = @/pi + (1 — P)/P2 (here applied between I? and L°°, i.e., with 
px = 2, p* = p, p 2 = oo). □ 
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C Appendix C. Proof of Proposition 15.151 



In this appendix, we complete the proof of Proposition 15.151 in the general 
symmetrizable, constant-multiplicity case; here, we make essential use of 
recent results of Metivier [Me. 4] concerning the spectral structure of matrix 
{A v )~ l (iT + L4>). Without loss of generality, take symmetric; this may 
be achieved by the change of coordinates — > A^ 2 A^A 1 ^ 2 . 

With these assumptions, the kernel and co-kernel of (A^° +to) are of fixed 
dimension m, not necessarily equal to one, and are spanned by a common set 
of zero-eigenvectors v\ , . . . , r m . Vectors v\ , . . . , r m are necessarily right zero- 
eigenvectors of (A 1 )~ l (iTo+iA^°) as well. Branch singularities correspond to 
the existence of one or more Jordan chains of generalized zero-eigenvectors 
extending up from genuine eigenvectors in their span, which by the argument 
of Lemma 3.4 is equivalent to 

(C.l) det(r t j A 1 r k ) = 0. 

In fact, as pointed out by Metivier [Me. 4], the assumption of constant mul- 
tiplicity implies considerable additional structure. 

Observation C.l ([Me.2]). Let (£o> r o) lie at a branch singularity involv- 
ing root ao = i£oi i n (|5,54|) . with tq an m-fold eigenvalue of ^4^°. Then, 
for (£, t) in the vicinity of (£o> r o)) the roots a bifurcating from ao i n <|5.54|) 
consist of m copies of s roots ai,...,a s , where s is some fixed positive 
integer. 

Proof of Observation. Let a(£, a) denote the unique eigenvalue of A^ lying 
near —to, where, as usual, — i£i := a; by the constant multiplicity assump- 
tion, a(-, •) is an analytic function of its arguments. Observing that 



(C.2) 



detail 1 )- 1 (it + iA*) - a] = det ^A 1 )' 1 det(r + A^) 

= e(i,r,a)(r + a(i,a)r, 



where e(-, •, •) does not vanish for (£, t, a) sufficiently close to (£o, To, ao), we 
see that the roots in question occur as m-fold copies of the roots of 

(C.3) T + a(i,a) = 0. 

But, the lefthand side of (|Q3|) is a family of analytic functions in a, con- 
tinuously varying in the parameters (£,t), whence the number of zeroes is 
constant near (£o> r o)- D 
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Observation C.2 ([Z.3]). The matrix (rjA r k ), j, k = 1, . . . , m is a real 
multiple of the identity, 

(C.4) {r)A x r k ) = (da/d^)I m , 

where a(£) denotes the (unique, analytic) m-fold eigenvalue of perturbing 
from —To. 

More generally, if 

(C.5) (da/dCi) = ■■■ = = °> (^a/dd) + 

at £o, then, letting ri(£), . . . ,r m (£) denote an analytic choice of basis for 
the eigenspace corresponding to a(£), orthonormal at (£o> r o)> we have the 
relations 

(C.6) (A 1 r 1 (T + ^ )r i)P = r i)P _ 1 , 

for 1 < p < s — 1, and 

(C7) (r^AVfe,^!) = p! {&>a/d$)I m , 

for 1 < p < s, where 

(C.8) r iiP := {-iy>p(jPrj/de x ). 

In particular, 

(C.9) {r ii0 ,...,r i>s _i}, j = 1, . ..,771 

is a right Jordan basis for the total zero eigenspace of (A 1 ) _1 (ro + A* ), for 
which the genuine zero-eigenvectors lj of the dual, left basis are given by 

(CIO) lj = {l/s\{d s a/di{))A 1 r j . 

Proof of Observation. Considering A^ as a matrix perturbation in £i, we 
find by standard spectral perturbation theory that the bifurcation of the in- 
fold eigenvalue To as £i is varied is governed to first order by the spectrum 
of (r t jA 1 rf : ). Since these eigenvalues in fact do not split, it follows that 
(rjA 1 rf t ) has a single eigenvalue. But, also, (rjA rf.) is symmetric, hence 
diagonalizable, whence we obtain result (jC,4|) , 

Result (jC.7|l may be obtained by a more systematic version of the same 
argument. Let denote the matrix of right eigenvectors 

(Cll) #(6):=(ri,...,r m )(6). 
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Denoting by 

(C.12) o(Ci + h) =: a + a l h + • • • + a p h p + ... 
and 

(C.13) + h) =: R° + R l h + • • • + RPhP + . . . 

the Taylor expansions of functions a(-) and R(-) around £o as £1 is varied, 
and recalling that 

(C.14) At = A*° + hA 1 , 

we obtain in the usual way, matching terms of common order in the ex- 
pansion of the defining relation (A — a)R = 0, the heirarchy of relations: 



(A& - a°)R° = 0, 
(A^-a^R 1 = -{A 1 - a 1 )R°, 
(C.15) 04 fo " a°)R 2 = -{A 1 - a^R 1 + a 2 R°, 



(A* - a°)R p = -{A 1 - a l )R p - 1 + a 2 R p ~ 2 + ■■■ + a p R° 



Using a = • • • = a s_1 = 0, we obtain l)C.6|) immediately, from equations 
p = 1, . . . , s — 1, and RP = (l/p\){r\ jP , . . . , r mjP ). Likewise, (|C.7|I . follows 
from equations p = l,...,s, upon left multiplication by L° := (R )^ 1 = 
(R°Y, using relations L (A& - a ) = and a p = (d p a/d^)/pl. 

From (|(J.6|) . we have the claimed right Jordan basis. But, defining lj as 
in IjC.lQI) . we can rewrite ()C.7|) as 



(C.16) {l)r k 



j'k,p-l) 



1 <p < S- 1, 



these ms criteria uniquely define lj (within the ms-dimensional total left 
eigenspace) as the genuine left eigenvectors dual to the right basis formed 
by vectors r^ v (see also exercise just below). □ 

Observation l( ) . 2l implies in particular that Jordan chains extend from all 
or none of the genuine eigenvectors r%, ... ,r m , with common height s. As 
suggested by Observation IC.ll (but not directly shown here) , this uniform 
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structure in fact persists under variations in £, r, see [Me. 4]. Observation 
EH is a slightly more concrete version of Lemma 2.5 in [Me. 4]; note the close 
similarity between the argument used here, based on successive variations in 
basis rj, and the argument of [Me. 4], based on variations in the associated 
total projection. 

With these preparations, the result goes through essentially as in the 
strictly hyperbolic case. Set 



(C.17) 
and define 
(C.18) 



p := l/{s\(d s a/dt{)) 



Note, as claimed, that p ^ by assumption (<9 s a/<9£f) ^ in Observation 
RT2l and sgn(p)Q > by (Kl), Proposition ITT81 

Thus, working in the Jordan basis defined in Observation IC.21 we find 
similarly as in the strictly hyperbolic case that the matrix perturbation 
problem 1)5.58(1 reduces to an ms x ms block-version 

(C.19) (U + icrM + pN -(a-aflvij = 

of (|A.2|) in the strictly hyperbolic case, where 







(° 


Im 


■•• 


o\ 












Im 




(C.20) 


J : = 








I m 













••• 


0/ 



denotes the standard block-Jordan block, and the lower-lefthand block of 
iaM+pN is apI m — ipQ ~ |o"| + [p[. To lowest order 0(|<7| + |p|) 1//s , therefore, 
the problem reduces to the computation of eigenvectors and eigenvalues of 
the perturbed block-Jordan block 



(C.21) 



diag{ r t 



Im 

I m 





\aplm - ipQ 



0/ 
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from which results (|5.62|) - (|5,66[) follow as before by standard matrix per- 
turbation theory; see, e.g., Section 2.2.4, Splitting of a block-Jordan block of 
[Z.4]. (Note that the simple eigenvalue case s = 1 follows as a special case 
of the block-Jordan block computation, with a = 0.) This completes the 
proof of Proposition 15.151 in the general case. 
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